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EDITORIAL 


Una objeción recurrente en algunos medios en lo que respecta a la Investigación científica tiene que 
ver con los temas que abordan los investigadores en sus estudios. Los proponentes de un lado de 
esta confrontación, manejan que no importa lo que se investigue si se cumple el objetivo fundamental 
de la investigación que es avanzar las fronteras de la ciencia. Esta posición no está exenta de 
méritos pues, por la esencia misma de la investigación, no se puede saber de antemano cuáles 
serán sus resultados. De hecho, muchos de los grandes descubrimientos científicos de la humanidad 
se han alcanzado por accidente (por ejemplo la penicilina), así que investigar por investigar tiene sus 
virtudes, sobre todo si los investigadores son talentosos. Esta estrategia es seguida por algunas de 
las más grandes empresas a nivel mundial, sobre todo por las del área de informática y de teleco- 
municaciones. 

Los defensores del otro punto de vista aducen que países con recursos limitados como México, no 
se pueden dar el lujo de desperdiciar sus escasos recursos en investigaciones sin objetivos claros 
(desde el punto de vista utilitario) o que no sean de aplicación inmediata. Esta postura también tiene 
visos de legitimidad, aunque tiene el defecto de dejar de lado buena parte de la investigación básica 
que no tiene aplicaciones con impacto económico inmediato. Las discusiones filosóficas sobre el 
propósito de la ciencia son en sí una rica veta de análisis y reflexión y no es ocioso involucrarse en 
ellas pues el hecho mismo de meditar y argumentar a favor o en contra de alguna de las posiciones 
citadas (o de alguna otra) enriquece el debate científico y contribuye a ensanchar sus miras y elevar 
sus metas. 

Este número es un interesante ejemplo de la diversidad en la investigación sobre temas forestales 
que actualmente se abordan en América Latina y su relación con los acuciantes problemas ambien- 
tales que aquejan a nuestros países. Los siete trabajos incluidos abarcan una amplia gama de inte- 
reses y en general están muy relacionados con los temas de actualidad en el ámbito científico 
forestal: la morfología de las fibras en madera de rápido crecimiento de plantaciones en un país 
centro americano, la distribución de una especie forestal de gran importancia en escenarios de 
cambio climático, el historial de incendios y la estructura forestal en la Sierra Madre Oriental de 
México, el ataque de plagas en frutos de especies forestales, el monitoreo de ecosistemas con apli- 
caciones de gran escala, el desarrollo de raíces y la supervivencia de propágulos en manglares y las 
características anatómicas y mecánicas de la madera. 

Con esta amplitud de temas, pareciera que todos los puntos de vista o corrientes de pensamiento 
tienen cabida en los temas de investigación vigentes. De lo que sí estamos seguros es que se trata 
de contribuciones relevantes y de gran calidad y rigor científico y que con toda certeza ayudarán a 
ir aclarando algunas de las dudas y falta de conocimiento que actualmente se tienen sobre estos 
temas. 


Raymundo Dávalos Sotelo 

Editor 
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ARTÍCULO DE INVESTIGACIÓN 


Fiber morphology in fast growth 
Gmelina arbórea plantations 

Morfología de la fibra en plantaciones de rápido crecimiento 

de Gmelina arbórea 

Róger Moya Roque 1 *, Mario Tomazelo Fo 2 
y Edwin Canessa Amador 3 

ABSTRACT 

Gmelina arbórea is planted ¡n large areas of forest with the objective of producing solid wood using 
well known silvicultura! techniques and taking advantage of the properties of the wood quality of fast-growing 
tree species managed in short rotation Systems. The aim of this study was to analyze the morphology and 
dimensión of fibers from the pith to the bark in trees from fast growth plantations in northern Costa Rica. 

The results indícate that fiber morphology is irregular in both diameter and shape; with 1 to 4 
septa, abundant crystals deposited in fiber lumina and minutely bordered pits with oblique and non- 
vestured apertures. Fiber length, width and cell wall thickness increased with tree age in the early 
stages. Lumen diameter was not correlated with tree age. 

KEY WORDS: 

Fiber width, lumen diameter, cell wall thickness, fiber length. 


RESUMEN 

Gmelina arbórea se planta en grandes áreas de bosques con el objetivo de producir madera 
sólida mediante el uso de técnicas silvícolas bien conocidas, así como para aprovechar las ventajas 
de las propiedades de la calidad de la madera de especies arbóreas de rápido crecimiento y mane- 
jadas por sistemas de rotación corta. El objetivo de este estudio fue analizar la morfología y dimen- 
sión de las fibras, desde la médula hasta la corteza, en los árboles de rápido crecimiento de las plan- 
taciones del norte de Costa Rica. Los resultados indican que la morfología de la fibra es irregular, tanto 
en diámetro como en forma; con 1 a 4 septa, cristales abundantes que se depositan en los lúmenes 
de la fibra y con hendiduras minuciosamente bordeadas con orificios no cubiertos y oblicuos. La 
longitud y ancho de la fibra y el grosor de la pared celular se incrementaron con la edad del árbol en 
las etapas tempranas. El diámetro del lumen no se correlacionó con la edad del árbol. 

PALABRAS CLAVE: 

Anchura de la fibra, diámetro del lumen, grosor de la pared celular, longitud de la fibra. 


INTRODUCTION 

Gmelina arbórea was introduced in species was chosen for its resistance to 

the Americas between 1970 and 1975, pests and diseases, its fast growth and 

mainly in Brazil (Kailsh, 1975). This the suitability of its wood for pulp and raw 

1* Centro de Investigación en Integración Bosque Industria, Instituto Tecnológico de Costa Rica, Apartado 
159-7050, Cartago-Costa Rica, rmoya@itcr.ac.cr. Autor para correspondencia. 

2 Departamento de Ciencias Florestais, ESALQ/Universidade de Sao Paulo, Av. Pádua Dias, 11-Caixa Postal 
09, CEP 13418-900-Piracicaba, Sao Paulo Brazil, mtomazel@esalq.usp.br 

3 Centro de Investigación en Integración Bosque Industria, Instituto Tecnológico de Costa Rica, Apartado 
159-7050, Cartago-Costa Rica, ecanessa@itcr.ac.ac 
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material for solid producís. Very soon after 
it was planted in the north of the Brazilian 
Amazon, it became evident that the trees 
grew very poorly at some sites (Hornick et 
al., 1984). In other countries, G. arbórea 
has been planted extensively to produce 
wood for solid products. In Costa Rica, G. 
arbórea is the most important species for 
timber production, and approximately 
65000 hectares have been planted in 
different ecological zones, from wet to dry 
(Moya, 2004) and under a variety of silvi- 
cultura! management regimes (Meza, 
1999). 

In terms of the wood’s anatomical 
composition, fibers are the predominant 
cell type and occupy 57-78% of the tissue 
volume (average: 76,6%). Several studies 
have shown that this percentage is not 
influenced by site, tree height, early and 
latewood, or radial distance from the pith 
(Hughes, 1968; Sosanwo and Lindberg, 
1975; Akachuku and Burley, 1979; Frim- 
pong-Mensah, 1992; Nobuchi et al., 
1997). The ¡mportance of fiber dimensión 
in making paper is well documented: the 
Runkel coefficient indicates this rela- 
tionship, but has the limitation of only 
taking cell wall thickness and lumen 
diameter into consideration (Sosanwo 
and Lindberg, 1975). Variations in other 
fiber dimensions, i.e., fiber length, wall 
area, cell wall thickness and lumen peri- 
meter and their effect on the properties of 
paper are mentioned in the literature as 
well (Corsan, 2002). Fiber dimensions are 
also important in solid wood products. For 
example, in Chínese poplar ( Popular 
trichocarpa) Bao and Liu (2001) detailed 
the relationship between fiber width, 
vessel percentage and other wood 
characteristics with the quality of veneer 
and plywood. 

The distribution, shape and dimen- 
sión of the fibers in the wood of Gmelina 
arbórea have been described in detail by 
Pearson and Brown (1932) based on a 


macroscopic examination of two Indian 
samples. They observed that the fibers 
were non-libriform, septate, with long 
tapering and numerous simple or 
bordered pits. Other studies report the 
anatomical features; mainly variation in 
fiber length (Hughes and Esan, 1969; 
Akachuku and Burley, 1979; Ohbayashi 
and Shiokura, 1989a; Frimpong-Mensah, 
1992). The cell and lumen diameter, and 
cell wall have been studied as well, by 
Frimpong-Mensah (1992) for example, 
and they found that cell wall thickness 
was significantly correlated with cambial 
age. Hughes and Esan (1969) found 
strong correlations between fiber length 
and tree age with distance from the pith in 
9-year-old trees in Nigeria. Also in 
Nigeria, for 7-year-old trees, it was found 
that fiber length was different at four sites 
(Akachuku and Burley, 1979). In contrast, 
Frimpong-Mensah (1992) found no varia- 
tion in fiber length with cambial age in 20- 
year-old Gmelina arbórea trees in Ghana. 
Growth rate affects fiber dimensions too. 
Ohbayashy and Shiokura (1989b) carried 
out a study on fiber length in 15-year-old 
trees and found that a high growth rate 
was strongly correlated with short fiber 
length. 

OBJECTIVE 

The objective of the present research was 
to study the variation in shape and size of 
wood fibers in Gmelina arbórea trees 
under fast growth conditions in Costa 
Rica. 

MATERIALS AND METHODS 

Sampled areas: G. arbórea plantations 
were selected in different regions of Costa 
Rica (Fig. 1), where mean annual precipi- 
tation varíes between 1700 and 5000 mm, 
temperatures range from 20 to 28 °C, and 
there is a modérate to short dry season (3 
to 1 dry months). 
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Figure 1. Location of the plantations sampled in northern Costa Rica 


Plot establishment and tree selec- 
tion: A total of 30 plots (400 m 2 each) were 
established in 9- to 12-year-old planta- 
tions. The coordinates of each plot were 
recorded with a Global Positioning 
System (gps). On each plot, the diameter 
at breast height (dbh), total and crown 
height were measured. One tree with an 
average dbh, a straight stem, and no 
disease or pest symptoms was selected 
on each plot (Table 1). 

Wood samples: Each sampled tree 
was marked on the side facing north and 
a stem section (3,0 cm thickness) was cut 


at dbh. A 1 ,0 mm wide, 2,0 mm thick slice 
was cut from pith to bark. Each growth 
ring was separated using a razor blade. 
We decided to only use the beginning part 
of the annual ring - which is similar to 
earlywood in températe species - given 
that the profile patterns of fiber dimen- 
sions are more or less stable from the 
beginning of the ring to 50% (Moya, 
2005). This part of the ring was macerated 
using Franklin’s method (Ruzin, 1999). 

Fiber morphology: Fiber shape and 
inclusions were identified under a light 
microscope. In addition, 6 block samples 
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Table 1. Description of plantation samples 


Plantation 

sampled 

Stand 
density 
(n ha' 1 ) 

Diameter 
at breast 
height (cm) 

Average 
growth rate 
(cm year' 1 ) 

Total 

height 

(m) 

Crown 

height 

(m) 

Age 

(years) 

1 

159 

38,0 

4,22 

15,5 

11,0 

9 

2 

300 

31,5 

3,94 

22,0 

11,0 

8 

3 

191 

33,2 

3,32 

18,5 

11,2 

10 

4 

223 

32,0 

4,00 

23,5 

14,5 

8 

5 

127 

30,7 

3,41 

18,0 

13,0 

9 

6 

659 

25,0 

3,13 

19,0 

6,0 

8 

7 

477 

28,5 

3,56 

22,5 

14,0 

8 

8 

350 

25,2 

3,15 

21,4 

14,0 

8 

9 

446 

30,0 

3,75 

23,2 

16,0 

8 

10 

505 

28,3 

2,83 

24,0 

18,0 

10 

11 

732 

24,0 

2,67 

21,8 

15,0 

9 

12 

1032 

21,5 

2,69 

21,5 

17,5 

8 

13 

764 

21,0 

1,91 

19,1 

11,6 

11 

14 

891 

24,1 

2,41 

15,0 

12,0 

10 

15 

1496 

24,3 

2,03 

20,0 

13,8 

12 

16 

318 

31,9 

3,19 

20,8 

11,0 

10 

17 

350 

29,8 

2,71 

19,5 

12,5 

11 

18 

223 

32,5 

2,71 

24,7 

13,5 

12 

19 

344 

30,5 

3,05 

24,0 

11,0 

10 

20 

250 

33,2 

2,77 

27,1 

17,1 

12 

21 

509 

19,7 

1,79 

19,5 

13,2 

11 

22 

477 

22,5 

2,50 

23,3 

13,4 

9 

23 

477 

23,4 

2,34 

17,8 

13,5 

10 

24 

605 

23,5 

2,61 

18,5 

7,7 

9 

25 

509 

23,9 

2,99 

22,1 

15,1 

8 

26 

827 

20,5 

2,05 

18,7 

13,0 

10 

27 

732 

18,8 

2,09 

19,0 

9,5 

9 

28 

836 

20,5 

2,05 

22,1 

14,2 

10 

29 

796 

20,7 

2,07 

18,3 

9,1 

10 

30 

1025 

20,2 

2,24 

22,0 

13,0 

9 


(5x5x5 mm) were cut from six different 
trees. They were dried at 80°C and 
coated with platinum for later observation 
with the scanning electronic microscope 
(sem) at the Electronic Microscopy Unit of 
the University of Costa Rica. 

Fiber dimensions: Initially, a sample of 
30 fibers was evaluated to determine the 
number of replicates it would be necessary 
to measure to obtain fiber length, width and 
lumen diameter. Thirty coefficients of varia- 
tion (cv) were calculated for the first fiber 
measured to the last one. A scatter plot was 
then graphed with the cv valúes and the 
measurements. The number of fibers to 
measure was determined as the point where 
the coefficient stabilized. It was necessary to 


measure 25 fibers for fiber width and lumen 
diameter, and 24 fibers for length. Each 
annual ring was macerated and three glass 
slides of the macérate prepared for observa- 
tion. Eight fibers were measured for length 
and five fibers for diameter by observing 
them under magnifying lenses at 250x and 
lOOOx, respectively, on each glass slide. 
Fiber wall thickness was calculated by subs- 
tracting lumen diameter from fiber width. 

Data analysis: Different models were 
tested for the analysis of variation in fiber 
dimensión from pith to bark, setting tree age 
as the independent variable and fiber 
dimensión as the dependent variable. An 
analysis of variance (anova) was carried out 
using the sas program (sas Institute, 1997). 
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RESULTS 

A total of 1 to 4 septa per fiber were found 
in the wood samples (Figs. 2a and 2b). 
Non-septate fibers were observed as well, 
though there were few of these. Crystals 
were found in several sections of the 
lumen or in each septum of the septa, 
when present (Figs. 2a and 2b). The 
fibers were irregular in both diameter and 
shape (Fig. 2b), with one- or two-sided, 
long tapering and circular forked ends 
(Figs. 2c and 2d). Few pits were found 
along the fiber walls and there was a 


mean of 20 mm between pits (Fig. 3a). 
Minutely bordered pits, with borders less 
than 2 mm in diameter (Fig. 3b), oblique, 
longitudinal apertures (Fig. 3c) and 
lacking vesturing (Fig. 3d) were observed 
in the wood samples. 

Dimensión fibers: The mean fiber 
length of G. arbórea wood samples was 
1,29 mm, fiber width was 30,67 mm, cell 
wall thickness was 4,01 mm, and lumen 
diameter was 22,64 mm (Table 2). Fiber 
length tended to increase with tree age 
(Fig. 4a), varying from 0,75 to 1,75 mm 



Figure 2. Gmelina arbórea fibers. a) Crystals both in the lumen and between septa, b) Cell 
with an irregular diameter (arrow indicates septum), c) One-sided forked ends (delineated 
by white dots), d) Two-sided forked ends (dashed line showing indicates end of fiber) 
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Figure 3. Gmelina arbórea pits in the cell wall. a) Pit (arrow) arrangement, 
b) Obligue aperture in longitudinal direction. Fig. 3c-d. Bordered pit. 


Scale bar 3a = 30 pin; 3b = 


(Table 2), but increased rapidly with tree 
age up to the 6 th year. After this, fiber 
dimensión decreased drastically and 
reached a plateau. The model tested for 
predicting fiber length from tree age had a 
high correlation coefficient (Table 3). Fiber 
wall thickness increased with tree age (Fig. 


10 pm; 3c and 3d = 1,5 pm. 


4b), varying from 1,74 to 5,83 mm (Table 
2), with an increase of 45% between the 1 st 
and 12 th year (Table 3). Fiber width incre- 
ased with tree age from 24,89 to 37,27 
mm; however the correlation coefficient 
decreased and was more variable within 
each growing year. The greatest increase 



Madera y Bosques 13(2), 2007:3-13 


9 


Table 2. Fiber dimensión of Gmelina arbórea from fast growth plantations in Costa Rica 


Fiber dimensión 

Minimum 

Máximum 

Mean 

Number 

of 

Samples 

Standard 

deviation 

Coefficient 
of variation 

Fiber length (mm) 

0,75 

1,75 

1,29 

7125 

0,210 

16,26 

Cell wall thickness 

1,74 

5,83 

4,02 

7125 

0,636 

15,84 

(nm) 







Fiber width (pm) 

24,89 

37,27 

30,67 

7125 

2,400 

7,72 

Lumen diameter (|um) 

17,57 

28,38 

22,64 

7125 

2,159 

9,54 


occurred between the 1 st and 6 th year of 
growth and after the 6 th year fiber width 
rate of change decreased or remained 
constant (Fig. 4c). Average lumen diameter 
was not strongly correlated with tree age, 
and variation within each growing year was 
high (Fig. 4d, Table 3). 

DISCUSSION 

Septate fibers have been previously 
reported for Gmelina arbórea by Pearson 
and Brown (1932), and Flughes and Esan 
(1969). These are fibers that contain a 


septum or a thin transverse wall across 
the lumen, which is present in dicotyledo- 
nous plants and is formed during cell wall 
maturation, in which cells lose their cyto- 
plasmic contents and lignification 
processes are absent (Itabashi et al., 
1999). Septa were not observed in all 
fibers and their variable occurrence (from 
1 to 4 septa per fiber) is suspected to be 
influenced by anatomical features such as 
tree age or the proximity of other woody 
cells (Itabashi et al., 1999); possible 
causes that are beyond the scope of the 
present research. 


Table 3. Parameter valúes and regression statistics for the best allometric models 
developed for fiber dimensions of G. arbórea from fast growth plantations in Costa Rica 


Fiber dimensión 

a 

b 

c 

Coefficient 

correlation 

RMSE 

P Valué 

Fiber length 1 

0,8306 

0,1327 

-0,0068 

0,6801 

0,119 

<0,0001 

Cell wall 


3,0136 

0,2416 

-0,0079 

0,4976 

0,4524 

<0,0001 

Thickness 1 

Fiber width 1 

27,3030 

0,7984 

-0,0505 

0,2772 

2,0483 

<0,0001 

Lumen diameter 2 

22,0810 

0,1043 

- 

0,0160 

2,1415 

00185 


1 Quadratic model: Y = a + b*age + c*(age ) 2 

2 Linear model: Y = a + b*age 
RMSE: Root Mean Square Error 
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Figure 4. Variation in the fiber dimensión of Gmelina arbórea from a fast growth 

plantation in Costa Rica 


Crystal deposition within the fiber 
lumen on secondary xylem has been 
reported for G. arbórea. Deshpande and 
Vishwakarma (1992), and Rao and Dave 
(1984) carried out a detailed analysis of 
the occurrence of crystals at fusiform and 
at radial/axial parenchyma cells in the 
secondary phloem and primary xylem of 
G. arbórea. According to those authors, 
the crystals are formed by deposits of 
calcium oxalate, which result from the 
high metabolic rate generated by a 
physiologically active vascular cambium. 
According to our results, it can be 
concluded that the crystals in the secon- 
dary xylem correspond to those deposited 
in the fusiform initial cells of the vascular 
cambium. 


The irregular diameter and shape of 
fibers in G. arbórea could have originated 
in the cambial zone. Rao and Dave (1985), 
and Dove and Rao (1982) have studied the 
causes of the wide variation in size, shape, 
and arrangement of the fusiform cells in the 
vascular cambium and primary/secondary 
xylem of G. arbórea. Vascular cambium 
activity in primary xylem was intense; the 
fusiform initials were long, radially flattened 
and tangentially tapering with overlapping 
ends. Múltiple división types (periclinal, 
anticlinal, pseudotransverse) were 
observed within radial and circumferential 
increments. 

According to our results, trees from 
fast growth plantations are expected to 
have the greatest variation in shape and 
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fiber diameter. On the other hand, the 
presence of circular forked ends usually 
occurs when those cell tips endose a ray, 
and forking appears during the elongation 
period of the fusiform cell in the cambial 
zone (Rao and Dave, 1985). This type of 
cell end is considered abnormal because 
once arrested at vertical ends, cell growth 
may take place at the lateral walls. 

Pearson and Brown (1932) described 
numerous and simple or bordered pits in 
G. arbórea from India, whereas in our 
study only bordered pits were found. 

The average and variation range for 
the fiber dimensions obtained in the 
present study coincide with those of 
previous studies of Nigerian wood in India 
by Pearson and Brown (1932). Akachuku 
and Burley (1979) recorded fiber lengths 
of 0,80 to 1,47 mm in Ghana; Ohbayashi 
and Shiokura (1989a) found that fiber 
length was 0,5 mm near the pith and 1,4 
to 1,6 near the bark in Malaysia, and 
Nobuchi et al. (1997) have reported less 
variation in fiber width, as supported by 
this research. However, Frimpong- 
Mensah (1992) reported fiber lengths of 
0,7 to 1 ,09 mm and wall thickness ranging 
from 2,70 to 3,38 mm for 16-year-old 
trees in Japan; all of which are lower than 
our valúes. 

G. arbórea trees increase fiber 
length, cell wall thickness and fiber width 
as they get older. This is a result of the 
many molecular and physiological 
changes that occur in the vascular 
cambium during the aging process 
(Plomion et al., 2001). The cells produced 
in the primary xylem divide less frequently, 
thus allowing more time for the fusiform 
initial section to elongate longitudinally 
and transversally (Horacek et al., 1999). 

At 6 years of age an inflection in the 
allometric model for length, width and 
fiber wall thickness suggests the begin- 


ning of a transition between juvenile and 
mature wood. In Gmelina arbórea this 
transition was previously identified at 4,0- 
9,0 cm from the pith by Ohbayashi and 
Shiokura (1989a). If the annual increment 
(caí) in dbh G. arbórea trees is 2,0 cm 
year 1 in Costa Rica (Moya, 2004), the 
juvenile period of this species can be set 
at 5-6 years oíd, when the inflection 
occurs in the fiber dimensión curve. It is 
necessary, however, to confirm this result 
for older trees. Hughes and Esan (1969) 
found strong correlations between fiber 
length, tree age and distance from the pith 
at 1,0 mm intervals from pith to bark in 
9-year-old trees in Nigeria. Akachuku and 
Burley (1979) found that fiber length 
increased until 6, 0-8,0 cm from pith to 
bark in 51-year-old Nigerian trees. In 
contrast, Frimpong-Mensah (1992) found 
no variation in fiber length with cambial 
age in 20-year-old G. arbórea trees in 
Ghana. 

The average fiber diameter incre- 
ased with tree age (Fig. 4c) but there was 
no variation in lumen diameter with tree 
age (Fig. 4d), and according to this result, 
the increase in fiber diameter was mainly 
the result of the increase in cell wall thick- 
ness (Fig. 4b) and not that of lumen 
diameter (Fig. 4b). 

CONCLUSIONS 

Characteristics were found in G. arbórea 
fibers from fast growth plantations, which 
have not been reported in the literature to 
date, i.e. the presence of minutely 
bordered pits and the lack of regular 
morphology related to its width and tips. 
Fiber length and width increased with 
increasing tree age. However, these fiber 
dimensions increased rapidly with tree 
age up to the 6th year, after which the rate 
of increase slowed down (and decreased) 
drastically until it reached a plateau. Wall 
thickness was associated with tree age 
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but no inflection in the rate of increase for 
the length or width of the fiber. The best 
allometric models were the linear model 
for lumen diameter and the quadratic 
model for width, length, and fiber wall 
thickness. 
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ARTÍCULO DE INVESTIGACIÓN 


Efecto de factores abióticos 
en el desarrollo de raíces primarias, 
crecimiento y supervivencia de propágulos 
en Rhizophora mangle L. 

The abiotic factors effects on the development 
of primary roots, growth and survival of Rhizophora mangle L. 

propagules 

José Luis Febles-Patrón*, Jorge Novelo López 
y Eduardo Batllori Sampedro 

RESUMEN 

Se realizó un experimento para conocer el efecto de los factores abióticos en el desarrollo de 
raíces primarias, crecimiento y supervivencia de propágulos de mangle rojo ( Rhizophora mangle L.). 

Se colectaron 180 propágulos en la ciénaga de Santa Clara, en la costa norte de Yucatán, México, y 
se formaron 6 grupos de 30 individuos. Los propágulos se colocaron en orificios de una lámina de 
poliestireno expandido y se mantuvieron durante 48 días en recipientes de plástico, en posición 
vertical y con el hipocótilo inmerso en agua. Se establecieron tres salinidades: dulce (5-7 ups), marina 
(30-40 ups) e hipersalina (70-75 ups), bajo dos condiciones de iluminación (sombra y radiación solar 
directa) con un máximo albedo. Las plantas que permanecieron bajo condiciones de sombra (SD, SM 
y SH) presentaron una sobrevivencia de 100%, mientras que en los tratamientos expuestos directa- 
mente a la radiación solar directa (LD, LM y LH) presentaron una sobrevivencia inversa a la salinidad. 

El desarrollo de las raíces fue gradualmente menor al aumentar la salinidad, disminuyendo cuando los 
propágulos se encontraron expuestos a la insolación directa y a una mayor temperatura del agua. Los 
propágulos sometidos al tratamiento sombra-dulce (SD), presentaron una mayor ganancia de peso 
por el desarrollo de raíces y hojas (sólo en este tratamiento se produjeron hojas). La salinidad es el 
factor que controla la producción de raíces en los propágulos de R. mangle en la etapa de estableci- 
miento. La insolación y la temperatura provocan un incremento del estrés fisiológico en los propá- 
gulos, disminuyendo la producción de raíces y aumentando la mortalidad. 

PALABRAS CLAVE: 

Fase de establecimiento, insolación, mortalidad, Rhizophora mangle L., salinidad del agua, Yucatán. 


ABSTRACT 

An experiment was made to know the abiotic factors effect on the development of primary roots, 
growth and survival of red mangrove propagules ( Rhizophora mangle L.). 180 propagules were 
collected at the Santa Clara swamp, in the North coast of Yucatán, México, forming 6 groups of 30 indi- 
viduáis. Propagules were placed in orífices of a white exponded polyestirene sheet, maintaining them 
during 48 days in plástic containers, with a vertical position and with hipocotil immersed in water. There 
were considered three salinities: freshwater (5-7 ups), sea water (30-40 ups) and hypersaline water 
(70-75 ups), and two conditions of illumination (shade and direct solar radiation) with a máximum 
Albedo, representing a total of six treatments. Plañís under shade conditions (SD, SM and SH) 
presented a survival of 100%, whereas the treatments exposed directly to solar radiation (LD, LM and 
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nico Nacional (CINVESTAV-IPN), Unidad Mérida. Km 6 antigua carretera a Progreso, Mérida, Yucatán, 
México. C.P. 97310. febles@mda.cinvestav.mx 
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LH) presented an ¡nverse survival to the 
water salinity. Roots development was 
gradually less when increasing the sali- 
nity, intensifying when the propagules 
were exposed to direct radiation and a 
greater temperature of the water. The 
propagules at the shade-freshwater treat- 
ment (SD), presented a greater gain of 
weight due to the development of roots 
and leaves (only in this treatment leaves 
were developed). Salinity is the factor that 
Controls the roots production in the R. 
mangle propagules at the establishment 
stage. Insolation and temperature cause 
an increase of physiological stress in the 
propagules, diminishing the production of 
roots and increasing mortality. 

KEY WORDS: 

Establishment phase, insolation, morta- 
lity, Rhizophora mangle L., water salinity, 
Yucatán. 

INTRODUCCIÓN 

Una de las actividades de restauración 
ambiental en los humedales costeros de 
México ha sido la reforestación del 
manglar. En este sentido, una de las 
técnicas propuestas es la siembra directa 
de propágulos de mangle rojo ( Rhizo- 
phora mangle L.) (Elster, 2000; Reyes y 
Tovilla, 2002). No obstante, hace falta 
información experimental que permita 
demostrar el efecto de factores como la 
salinidad, temperatura, radiación solar, 
inundación, nutrimentos, etc., en el esta- 
blecimiento, crecimiento y supervivencia 
de los propágulos. Es importante consi- 
derar estos factores en las pruebas de 
reforestación con propágulos y plántulas 
de R. mangle, debido a que la mortalidad 
es alta durante la fase de establecimiento, 
justo cuando las raíces están desarrollán- 
dose (Elster et al. 1999; Elster, 2000). 

Los manglares son dependientes 
del reclutamiento de plántulas para su 
mantenimiento (Craighead, 1971; 
Tomlinson, 1994); por lo tanto, el estable- 
cimiento y sobrevivencia de las plántulas 
afecta directamente los patrones de 
distribución y abundancia (McKee, 1995 
a y b). Existe una amplia tolerancia de los 


manglares a diferentes condiciones 
ambientales que se expresan en 
respuestas funcionales y estructurales 
(Jiménez, 1994). La salinidad del agua 
superficial e intersticial son factores 
importantes que pueden afectar el esta- 
blecimiento, la zonación, estructura y la 
supervivencia (Elster et al., 1999). 

Las especies de mangle presentan 
diversas adaptaciones fisiológicas para 
tolerar la salinidad del suelo. Entre éstas, 
podemos mencionar la exclusión y excre- 
ción de sales que son muy importantes en 
el balance osmótico de la planta (Jiménez, 
1 994). Las especies exclusoras son aque- 
llas que limitan el ingreso de sales a las 
plantas hasta en un 99%, mientras que las 
especies excretoras permiten un mayor 
ingreso de sales al sistema circulatorio, 
secretándola posteriormente a través de 
glándulas ubicadas en la superficie foliar. 
R. mangle L. es una especie exclusora y 
es menos tolerante a la salinidad del suelo 
(máximo 65 ups), comparada con especies 
como Avicennia germinans L. que puede 
tolerar hasta 100 ups, debido a su adapta- 
ción secretora (Jiménez, 1994). 

Elster et al. (1999), reportan para la 
Ciénaga Grande de Santa Marta, 
Colombia, que en áreas de bosque el creci- 
miento de plántulas está controlado princi- 
palmente por los cambios en la salinidad, 
mientras que la inundación estacional 
regula la mortalidad. En áreas desprovistas 
de vegetación, la regeneración es afectada 
por la salinidad extremadamente alta, altas 
temperaturas y desecación de la superficie 
del suelo. En este sentido, las condiciones 
físicas extremas pueden ser críticas para la 
sobrevivencia de plántulas de manglar 
(Rabinowitz, 1978). 

Así, en este estudio se realizó un 
experimento para conocer el efecto de la 
salinidad y la insolación, en el desarrollo 
de raíces primarias, crecimiento y super- 
vivencia en propágulos de mangle rojo. 
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ÁREA DE TRABAJO 

El trabajo experimental se realizó en la 
ciudad de Mérida, Yucatán, ubicada al 
sureste de México. El clima en la zona 
corresponde al subtipo Awo(i’)g (García, 
1988), con una marcada estacionalidad 
pluvial concentrada durante el período 
mayo-octubre. Por lo regular septiembre 
es el mes más lluvioso (186,9 mm 
promedio) y abril es el más seco (15,8 
mm promedio). La precipitación media 
anual varía de 838 a 1128 mm, con un 
porcentaje de lluvia invernal entre 6 y 
10,2% (Duch, 1988). La temperatura 
media anual varía entre 24,5 y 27°C, con 
un máximo promedio de 29,5°C en el mes 
de mayo y mínima promedio de 23,4°C 
en enero. No obstante, puede presentar 
una temperatura máxima extrema de 
39,9°C en mayo y una mínima extrema 
de 12,6°C en enero (Duch, 1988). 

METODOLOGÍA 

Durante el mes de agosto de 2004 se 
colectaron 180 propágulos de mangle rojo 
(. R . mangle) en la ciénaga de Santa Clara, 
en la costa norte de Yucatán. Cada propá- 
gulo fue pesado en una balanza analítica 
electrónica, modelo Ohaus Navigator, y 
se obtuvo el Peso Inicial (P¡). Con estos 
propágulos se formaron seis grupos de 
30 individuos para someterlos a trata- 
mientos con diferentes condiciones de luz 
y salinidad. Para controlar la intensidad 
de luz sobre los propágulos se instaló una 
tela de mallasombra (número 80 color 
verde) en un área de 6 m 2 . Para controlar 
la salinidad del agua en los tratamientos, 
se colectó agua dulce de manantial (5 a 7 
ups), agua marina (35 a 40 ups) y agua 
hipersalina (70 a 75 ups) en la costa 
yucateca. La salinidad fue medida usando 
una sonda YSI Modelo 85, con una 
amplitud de 0 a 80 ups y una exactitud de 
± 0,1 ups (Unidades potenciales de sali- 
nidad, equivale a partes por mil o gramos 
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por litro de agua). Se utilizaron seis reci- 
pientes de plástico de 1 m de largo por 40 
cm de ancho y 20 cm de alto. Posterior- 
mente se vertió agua con diferentes 
concentraciones de sal en cada uno de 
los recipientes hasta lograr una columna 
de agua de 7 cm de profundidad. En cada 
recipiente, se colocó sobre el agua una 
lámina de poliestireno expandido de 1 cm 
de grosor, a la cual se le realizaron treinta 
orificios de 1,5 cm de diámetro. Los 
propágulos se colocaron en cada orificio 
de la lámina y se mantuvieron en posición 
vertical y con el hipocótilo inmerso en 
agua durante 48 días. Para evitarles 
daños por el roce con la lámina se colocó 
alrededor de cada propágulo, a la altura 
de la lámina, una esponja sintética de 1 
cm de ancho. 

Tres recipientes se colocaron en el 
área acondicionada con la mallasombra y 
tres permanecieron expuestos a la luz solar 
directa. De esta manera quedaron estable- 
cidos los siguientes seis tratamientos: (1) 
sombra - agua dulce (sd), (2) sombra - 
agua marina (sm), (3) sombra - agua hiper- 
salina (sh), (4) luz - agua dulce (ld), (5) luz 
- agua marina (lm) y (6) luz - agua hipersa- 
lina (lh). 

Los propágulos se colocaron en los 
recipientes el 2 de agosto de 2004 (Tiempo 
1) y para mantener el nivel de salinidad del 
agua en cada tratamiento, se realizaron 
cambios de agua constantes durante 48 
días (Tiempo 2), fecha en que se realizó la 
evaluación de la respuesta de los propá- 
gulos sometidos a los seis tratamientos. 
Los parámetros que se utilizaron para la 
evaluación fueron: sobrevivencia, número 
de raíces, largo de la raíz, número de 
hojas, peso inicial (P¡) vs. peso final (P f ) y 
balance de biomasa neto. 

Para estimar la variación diurna de la 
temperatura del agua en cada tratamiento, 
se realizó un ciclo de 24 horas el 9 de 
septiembre, en el que se realizaron medi- 


18 


Efecto de factores abióticos en el desarrollo de raíces primarlas... 


ciones cada cuatro horas utilizando la 
sonda YSI 85 (rango de -5 a 65°C ±0,1). 

Se realizó una prueba de Kruskal- 
Wallis para muestras independientes del 
número y largo de las raíces, además de la 
biomasa neta obtenida entre los trata- 
mientos. También se realizó una prueba de 
Mann-Whitney para la comparación del 
peso inicial y final de cada tratamiento. 
Para comparar la temperatura obtenida en 
el ciclo de 24 h, se realizó otra prueba de 
Kruskal-Wallis, además de una correlación 
no paramétrica de Spearman entre los 
valores promedio de la temperatura, la 
sobrevivencia, el número y largo de raíces, 
así como la biomasa neta obtenida para 
cada tratamiento. 


RESULTADOS 

Sobrevivencia 

Los propágulos que permanecieron bajo 
condiciones de sombra (sd, sm y sh) 
presentaron una sobrevivencia de 100%, 


mientras que en los tratamientos expuestos 
directamente a la insolación, se observó un 
gradiente de sobrevivencia inverso a la 
concentración salina; en el tratamiento ld el 
valor de sobrevivencia fue de 53,33%, en 
lm de 40% y en LH de 0% (Figura 1). 


Número de raíces 

Al aumentar la salinidad se incrementó la 
inhibición en el desarrollo de las raíces 
en los propágulos, aumentando cuando 
los propágulos se encontraron 
expuestos a la insolación directa y a una 
mayor temperatura del agua. Se regis- 
traron diferencias significativas entre los 
tratamientos, debido principalmente a la 
mayor producción de raíces en sd y sm 
(Figura 2). 

Todos los propágulos del trata- 
miento sd presentaron raíces. La 
mayoría entre tres y cuatro raíces (40 y 
46,7%, respectivamente), una proporción 
relativamente baja presentó cinco raíces 
(6,7%) y los demás presentaron entre 


100 100 100 



Sombra Luz dulce Sombra Luz marina Sombra Luz 

dulce marina hipersalina hipersalina 


Figura 1. Sobrevivencia obtenida en los tratamientos del experimento 
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Tratamientos 


Figura 2. Prueba de Kruskal-Wallis del número de raíces obtenido en el experimento 


una y dos raíces (6,6%). En el trata- 
miento sm no todos los propágulos 
presentaron raíces. El 66,7% desarrolló 
dos o tres raíces, el 20% presentó cuatro, 
el 3,3% desarrolló una, mientras que el 
3,3% no presentó ninguna raíz. En el 
tratamiento SH el 23,3% presentó única- 
mente una raíz, pero el 30% no presentó 
ninguna. De los propágulos sometidos al 
tratamiento ld, el 46,6% no presentó 
desarrollo de raíces. En el tratamiento lm 
la proporción de propágulos sin raíces 
aumentó al 60%, mientras que en el 
tratamiento LH ninguno desarrolló raíces 
(Figura 3). 

Largo de las raíces 

Se encontraron diferencias significativas 
en el largo de las raíces entre los trata- 
mientos. El largo promedio de la raíz fue 


menor en los tratamientos sometidos a la 
insolación directa, para cada condición 
de salinidad. El agua marina (con salini- 
dades de 35 ups) estimuló el mayor 
crecimiento longitudinal de las raíces de 
los propágulos, en ambas condiciones 
de luz y sombra (Figura 4). Para el trata- 
miento SM el valor promedio fue de 5,49 
cm (EE=±0,37) y para lm de 2,17 cm 
(EE=±0,51). Los tratamientos con agua 
dulce ocuparon un lugar intermedio entre 
los tratamientos marinos e hipersalinos; 
para SD el valor promedio fue de 4,16 
cm (EE=±0,20) y 1 ,46 cm (EE=±0,27) en 
ld. Los tratamientos con aguas hipersa- 
linas, ocasionaron que los propágulos 
presenten poco o nulo desarrollo de 
raíces; para el tratamiento SH el valor 
promedio del largo de la raíz fue de 0,61 
cm (EE=±0,09) y para LH el valor fue de 
0 cm. 
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□ 5 q4 o3 q2 □ 1 DÜ 


Figura 3. Proporción del número de raíces que presentaron los propágulos en cada 
tratamiento: sombra - agua dulce (SD), sombra - agua marina (SM), 
sombra - agua hipersalina (SH), luz - agua dulce (LD), luz - agua marina (LM) 

y luz - agua hipersalina (LH) 



Tratamientos 


Figura 4. Prueba de Kruskal-Wallis del largo de las raíces obtenido en el experimento 
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Peso inicial (P¡), Peso final (P f ) y 
balance de biomasa neto 

Tanto los tratamientos de luz como los 
de sombra presentaron un mayor 
balance de biomasa neto a menor sali- 
nidad del agua. Los propágulos del trata- 
miento SD presentaron un P¡ promedio 
de 11,17 g (EE=±0,56) y un P f de 13,04 
g (EE=±0,67), presentando un incre- 
mento promedio neto de 1,87 g 
(EE=±0,13), correspondiente al 16,73%. 
Para el tratamiento SM el valor promedio 
de P¡ fue de 9,96 g (EE=±0,41) y P f de 
10,25 g (EE=±0,46), obteniendo una 
ganancia promedio en peso de 0,29 g 
(EE=±0,41), correspondiente al 2,94 %. 
Los propágulos del tratamiento SH 
presentaron un P¡ de 13,47 g (EE=±0,67) 
y un P f de 10,40 g (EE=±0,59), regis- 
trando una pérdida de 3,06 g 
(EE=±0,20), equivalente al 22,7%. El P¡ 
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de los propágulos del tratamiento ld fue 
de 12,54 g (EE=±0,50) y el P f de 13,17 g 
(EE=±0,55), con una ganancia promedio 
de peso de 0,63 g (EE=± 0,21) (5,03%). 
En el tratamiento LM los propágulos 
obtuvieron un P¡ de 11 ,55 g (EE=±0,40) y 
un P f de 11,04 g (EE=±0,57), represen- 
tando una pérdida de peso de 0,50 g 
(EE=±0,31) (4,38%). Para el tratamiento 
en condiciones extremas LH, el P¡ 
promedio fue de 11 ,87 g (EE=±0,43) y el 
P f de 7,57 g (EE=±0,32), teniendo una 
pérdida considerable de 4,30 g 
(EE=±0,15), equivalente al 36,21% 
(Figura 5). 

La prueba de Mann-Whitney entre 
los pesos iniciales (P¡) y finales (P f ) de 
los propágulos, indicó diferencias signi- 
ficativas en tres tratamientos: sd, sh y lh 
(Tabla 1). Para el caso del tratamiento 
sd, la diferencia significativa ( p<0,03 ) se 



Tratamientos 


Figura 5. Peso inicial (1) y final (2) para los diferentes tratamientos del experimento 
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Tabla 1. Prueba de Mann-Whitney para los pesos inicial y final de cada tratamiento 

utilizado en el experimento 



SD 

LD 

SM 

LM 

SH 

LH 

U de Mann-Whitney 

303,000 

394,000 

398,000 

396,000 

237,000 

54,000 

W de Wilcoxon 

768,000 

859,000 

863,000 

861,000 

702,000 

519,000 

Z 

-2,173 

-,828 

-,769 

-,798 

-3,149 

-5,855 

Sig. asintót. (bilateral) 

,030 

,408 

,442 

,425 

,002 

,000 


a Variable de agrupación: peso 1 y 2. 


debió a una considerable ganancia de 
peso por el desarrollo de raíces y hojas 
(sólo en SD se produjeron hojas), por 
encontrarse los propágulos en condi- 
ciones de bajo estrés salino y térmico; 
en el caso contrario, los tratamientos sh 
y lh presentaron diferencias significa- 
tivas ( p<0,002 y p<0,001) debido a la 
pérdida de peso de los propágulos por 
deshidratación al encontrarse en condi- 
ciones ambientales extremas. De 
manera general, los resultados indican 
que los propágulos responden al 
gradiente que va desde el tratamiento 
sd hasta lh y que la diferencia entre 
esta amplitud, determina la sobrevi- 
vencia y la ganancia o pérdida de 
biomasa (Figura 6). 


Temperatura del agua 

A través del ciclo de 24 h se observó un 
gradiente de menor a mayor tempera- 
tura, salinidad e insolación. Para los 
tratamientos que se encontraban en 
sombra los valores máximos fueron los 
siguientes: SD= 31,70°C, SM= 31,90°C 
y SH= 33,20°C (entre 16 y 20 h). En los 
tratamientos que se encontraban direc- 


tamente a la luz, la temperatura 
aumentó notoriamente (3,4 °C) con 
relación a los que se encontraban en la 
sombra; el tratamiento ld presentó un 
valor máximo de 35,90°C, para lm fue 
de 35,20°C y para la condición extrema 
de luz y salinidad (LH) el valor fue de 
36,60°C (entre 16 y 20 h) (Figura 7). 

A través de la prueba de Kruskal- 
Wallis no se encontraron diferencias 
significativas entre sitios. No obstante, 
es evidente la menor temperatura en 
los tratamientos a la sombra, contra los 
valores máximos registrados bajo 
condiciones de insolación directa 
(Figura 8). 

En la tabla 2 se observa de 
manera detallada que existe una corre- 
lación inversa y altamente significativa 
entre los valores promedio de la tempe- 
ratura del agua con la sobrevivencia y 
con el número de raíces encontrado en 
los tratamientos. Por otra parte, es inte- 
resante notar que la sobrevivencia se 
relaciona de forma directa y significa- 
tiva con el número promedio de raíces 
en los tratamientos. 
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Figura 6. Prueba de Kruskal-Wallis del balance neto de biomasa en los tratamientos del 

experimento 
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Figura 7. Fluctuación de la temperatura del agua ( C) a través de un ciclo de 24 h 

(9 de septiembre de 2004) 
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Figura 8. Prueba de Kruskal-Wallis para la temperatura del agua en los tratamientos del 

experimento 


Tabla 2. Correlación de Spearman entre los valores promedio de los parámetros 
obtenidos en los tratamientos del experimento 




Temperatura 

Sobrevivencia 

Número de raíces 

Largo de raíces 

Biomasa 

Temperatura 

Coeficiente de correlación 

1 





Sig. (unilateral) 

** 

** 



N 

6 





Sobrevivencia 

Coeficiente de correlación 

-0,94 

1 




Sig. (unilateral) 

0,003 


** 



N 

6 

6 




Número de raíces 

Coeficiente de correlación 

-1 

0,94 

1 



Sig. (unilateral) 

0,000 

0,003 



N 

6 

6 

6 



Largo de raíces 

Coeficiente de correlación 

-0,71 

0,58 

0,71 

1 



Sig. (unilateral) 

0,055 

0,115 

0,055 




N 

6 

6 

6 

6 


Blomasa 

Coeficiente de correlación 

-0,71 

0,52 

0,71 

0,71 

1 


Sig. (unilateral) 

0,055 

0,147 

0,055 

0,055 



N 

6 

6 

6 

6 

6 



** La correlación es significativa al nivel 0,01 (unilateral). 
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DISCUSIÓN 

Los manglares pueden establecerse en 
condiciones variables de salinidad. No 
obstante, en la región del Caribe y Golfo 
de México, R. mangle se desarrolla mejor 
en salinidades cercanas al agua de mar o 
menores (Lugo y Snedaker, 1974; Cintrón 
et al., 1978; Elster, 2000). En general, 
varios estudios muestran que el óptimo 
desarrollo de R. mangle se alcanza a sali- 
nidades más bajas que el agua de mar 
(< 35 ups) (Medina, 1998; Marchand et 
al., 2004). Davis (1940) describió que los 
propágulos de R. mangle producen raíces 
en menor tiempo en agua dulce que en 
agua marina. Pannier (1959) registró una 
máxima producción de hojas y raíces en 
plántulas de R. mangle de seis meses de 
edad con una salinidad del agua de 9 ups. 
Por su parte, Werner y Stelzer (1990) 
confirmaron un incremento en el creci- 
miento de raíces de R. mangle con salini- 
dades de 14 ups. 

Nuestros resultados indican que el 
agua hipersalina (70 ups) y la exposición 
directa a la insolación, inhiben el 
desarrollo de las raíces en los propá- 
gulos. Por el contrario, en el tratamiento 
de sombra-dulce (sd) y luz-dulce (ld) se 
observó el desarrollo de raíces vigorosas. 
El vigor de las raíces del tratamiento sd 
favoreció el desarrollo de hojas; después 
de “30 días” iniciaron la formación de un 
par de hojas y el 86,6% de los propágulos 
concluyeron su desarrollo para el final del 
experimento (48 días). 

De acuerdo con Reyes y Tovilla 
(2002), los factores ambientales como la 
temperatura, insolación y evaporación 
contribuyen a elevar la salinidad, dismi- 
nuir la disponibilidad de agua e incre- 
mentar la mortalidad de las plantas. Por 
su parte, Lugo y Snedaker (1974) y 
Clough et al. (1982), mencionan que las 
temperaturas óptimas para el desarrollo 
adecuado de los procesos fotosintéticos 
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de los manglares se registran alrededor 
de los 35°C. En este estudio encontramos 
una correlación inversa significativa entre 
la temperatura y el número promedio de 
raíces en los tratamientos. 

Con relación al efecto de la insola- 
ción, Tovilla y Orihuela (2002) observaron 
el desarrollo de propágulos de R. mangle 
en condiciones de sombra (50% de ilumi- 
nación) e iluminación directa (100%) 
encontrando durante el primer año que la 
mayor mortalidad se registró al sol, ya 
que el 49,5% perecieron; en comparación 
con el 38,5% de aquellas desarrolladas a 
la sombra. Al respecto, McKee (1995c) 
indica que el proceso de arraigo 
temprano de los propágulos está más 
relacionado con la disponibilidad de 
nutrientes en el medio que con la luz. La 
mortalidad al sol, está relacionada 
también con el flujo de mareas a lo largo 
del año (Tovilla y Orihuela, 2002). Ellison 
y Farnsworth (1993) en Belice, obser- 
varon que R. mangle presentó mayor 
sobrevivencia relacionada con una mayor 
influencia de las mareas presentes en 
cada sitio, independientemente de si las 
plántulas estuvieran bajo el dosel del 
bosque o al sol. 

Los resultados de Smith y Lee 
(1999), quienes observaron que la 
cantidad de luz afecta la producción de 
raíces en R. mangle, indican que puede 
existir un mejor desarrollo a la sombra en 
la fase inicial de desarrollo. En este 
sentido, es posible que existan asocia- 
ciones positivas importantes entre indivi- 
duos o nodricismo, protegiendo a los 
propágulos de la insolación y tempera- 
turas extremas (Lema etal., 2003), lo cual 
sugiere que el establecimiento podría ser 
un factor clave en la formación de 
patrones agregados o parches de vegeta- 
ción, tan comúnmente observados en las 
ciénagas y lagunas costeras de México. 
Así, cualquier proceso de restauración 
ambiental en ecosistemas de manglar 
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debe considerar la sombra como un 
factor críticamente importante, principal- 
mente durante la fase de establecimiento. 

CONCLUSIONES 

Los propágulos sometidos al tratamiento 
sombra-dulce (sd), presentaron una 
mayor ganancia de peso por el desarrollo 
de raíces y hojas (sólo en este trata- 
miento se produjeron hojas), por encon- 
trarse los propágulos en condiciones de 
bajo estrés salino. En el caso extremo, los 
tratamientos sombra-hipersalino (sh) y 
luz-hipersalino (lh), presentaron diferen- 
cias significativas en comparación con los 
otros tratamientos debido a la pérdida de 
peso de los propágulos por deshidrata- 
ción, al encontrarse en condiciones de 
salinidad extrema. 

Con base en las observaciones 
realizadas en este experimento, podemos 
corroborar que la salinidad es el factor 
que controla la producción de raíces en 
los propágulos de R. mangle en la etapa 
de establecimiento. La insolación provoca 
un incremento del estrés fisiológico en los 
propágulos, disminuyendo la producción 
de raíces y aumentando la mortalidad. 
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ARTÍCULO DE INVESTIGACIÓN 


Distribución del cedro rojo 
(Cedrela odorata L.J en el estado de Hidalgo, 
bajo condiciones actuales y escenarios 

de cambio climático 

Distribution of red cedar (Cedrela odorata L.) 
in the Hidalgo State, under current conditions 
and scenarios of climate change 

Jesús David Gómez Díaz 1 , Alejandro Ismael Monterroso Rivas 1 

y Juan Ángel Tinoco Rueda 1 

RESUMEN 

Se obtuvieron las áreas climáticas potenciales de distribución del cedro rojo ( Cedrela odorata 
L.) en el estado de Hidalgo bajo condiciones actuales, con información meteorológica de 1961-1990 
como escenario base, y bajo escenarios de cambio climático obtenidos con los modelos GFDL-R30 y 
el HadCM3 para los años 2020 y 2050. Se aplicaron las tasas de cambio en temperatura y precipita- 
ción sobre las áreas de influencia climática delimitadas de acuerdo a Gómez et al. (2006). Se realizó 
un balance de humedad bajo condiciones actuales y de cambio climático aplicando la metodología de 
Thornthwaite modificada versión III (Monterroso y Gómez, 2003). Actualmente se reporta 9,8% de la 
superficie del estado con tipos de vegetación con los que se asocia el cedro rojo, y las estimaciones 
con las características climáticas del escenario base muestran que el 30,4% del estado presenta algún 
grado de aptitud para el desarrollo de esta especie. Los resultados del modelo GFDL-R30 indican un 
incremento en la superficie con algún nivel de aptitud, con respecto al escenario base, de 3,1 % y 4,4% 
para el año 2020 y 2050, respectivamente. El aumento es diferencial dentro de las clases de aptitud. 
Para el modelo HadCM3 la superficie con algún nivel de aptitud disminuye 0,9% para el año 2020 y 
0,2% para el año 2050. Sin embargo, el nivel de Moderadamente apto pasa de 10,5% en el escenario 
base a 0% y 1,3% para el año 2020 y 2050, respectivamente, ubicándose prácticamente todas las 
áreas en el menor nivel de aptitud. 

PALABRAS CLAVE: 

Balance de humedad, cambio climático, cedro rojo, zonificación agroclimática. 


ABSTRACT 

The potential climatic areas for the distribution of red cedar ( Cedrela odorata L.) were obtained in 
the State of Hidalgo under actual conditions, with meteorological information of 1961-1990 as base line, 
and under climate change conditions wich were obtained with two models: GFDL-R30 and HadCM3 for 
two time scenarios (2020 and 2050). The adjustment rates in temperature and rainfall obtained for each 
model were applied on climatic influence areas delimitated according to Gómez et al. (2006). Asoil mois- 
ture balance in the state’s surface was calculated under current and climate change conditions using the 
Thornthwaite, modified versión III methodology (Monterroso and Gómez, 2003). The present percentage 
of the State surface with the vegetation types associated with red cedar is of 9,8% and the surface percen- 
tage estimated with some degree of suitability for this specie, using the climatic characteristics of the 
baseline scenario, is of 30,4%. The results of applied GFDL-R30 model shows an increment on the total 
surface with some degree of suitability with respect of the baseline scenario of 3,1% and 4,4% for the 

1 Universidad Autónoma Chapingo, Departamento de Suelos, Km. 38.5 carretera México-Texcoco. Chapingo. 
México, C.P. 56230. dgomez@correo. chapingo. mx 
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years 2020 and 2050, respectively, with a 
differential increment within the suitability 
classes. The surface estimated with some 
degree of suitability applying the HadCM3 
model shows a decrement of 0,9% for the 
year 2020 and 0,2% for 2050. However, 
the class of Modérate suitable, change 
from 10,5% on the baseline scenario to 
0% and 1,3% for the years 2020 and 
2050, respectively, with almost all the 
areas in the lowest level of suitability. 

KEY WORDS: 

Water balance, climate change, red 
cedar, agro climatic zoning. 

INTRODUCCIÓN 

La especie objeto de estudio en este 
trabajo, el cedro rojo ( Cedrela odorata L.), 
se localiza primordialmente en ecosis- 
temas húmedos, en específico en el 
bosque tropical perennifolio, el bosque 
tropical subcaducifolio, el bosque tropical 
caducifolio y el bosque mesófilo de 
montaña (conafor, 2000). Partiendo de 
esta información y de la distribución 
reportada en el inventario nacional 
forestal para estos tipos de vegetación 
(semarnat-unam, 2001), se estima que el 
cedro rojo en el estado de Hidalgo, se 
encuentra distribuido en áreas fragmen- 
tadas en la zona de barlovento que 
comprende parte de la planicie costera 
del Golfo norte, los flancos al este de la 
Sierra Madre Oriental hasta alturas inter- 
medias menores a los 2000 msnm. La 
superficie de los tipos de vegetación con 
los que se asocia el cedro rojo cubren 
2040,3 km 2 lo que equivale al 9,8% del 
territorio total del estado de Hidalgo. 

La distribución de las comunidades 
vegetales está determinada por diferentes 
factores ambientales como las caracterís- 
ticas del suelo y el clima; siendo este 
último el principal responsable (Zunino y 
Zullini, 2003; Berry et al., 2002). Cuando 
alguno de dichos factores cambia, en 
especial las condiciones climáticas, la 
estructura y composición de los ecosis- 


temas también pueden sufrir modifica- 
ciones. La respuesta por parte de las 
plantas y animales depende de la inten- 
sidad y velocidad del cambio en las condi- 
ciones climáticas, si el cambio es paulatino 
y con poca intensidad, es probable que 
gran parte de las especies logre evolu- 
cionar y sus descendientes se adapten a 
las nuevas condiciones ambientales. Sin 
embargo, si el cambio se produce con gran 
intensidad y en un periodo relativamente 
corto, un bajo porcentaje de la biodiver- 
sidad logra adaptarse a estos cambios, lo 
que produce la extinción de las especies 
con baja capacidad de amortiguamiento. 
Esta capacidad de adaptación de las espe- 
cies se ha venido presentando desde la 
aparición de la vida sobre el planeta Tierra 
(Pérez, 2005). 

La temperatura y la precipitación son 
dos de los elementos climáticos determi- 
nantes para el establecimiento de una 
especie y/o una comunidad vegetal en 
determinado sitio (Velázquez, 2002). Por 
su parte, Ortiz (1987) establece que la 
mayoría de los procesos fisiológicos que 
se realizan para el crecimiento y desarrollo 
de las plantas están fuertemente influen- 
ciados por la temperatura. En términos 
generales, ésta ejerce su influencia prin- 
cipal controlando la proporción de reac- 
ciones químicas involucradas en varios 
procesos de crecimiento dentro de la 
planta. La temperatura, además, afecta los 
mecanismos hormonales involucrados en 
la floración y fructificación de las plantas. 

El agua por su parte, resulta indispen- 
sable para la existencia y evolución de los 
ecosistemas vegetales y animales. Es un 
elemento asociado a condiciones climá- 
ticas y edafológicas, por lo que su cantidad 
en la atmósfera y en el suelo tiene una 
repercusión directa en los ecosistemas. El 
agua es la parte esencial del protoplasma 
vegetal, constituye un elemento importante 
de la fotosíntesis y un medio indispensable 
de suministro de elementos nutritivos. En 
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forma de vapor, el efecto principal del agua 
consiste en dar humedad al aire; un alto 
contenido de humedad del aire reduce la 
evaporación de los suelos así como la 
evapotranspiración de las plantas (Agui- 
lera y Martínez, 1996). 

Los cambios en el sistema climático 
global han ocurrido durante toda la 
historia del planeta. Dichas modifica- 
ciones se han presentado por causas 
naturales como por ejemplo: cambios en 
la órbita terrestre, alteraciones en la 
excentricidad del planeta, actividad volcá- 
nica intensa e impacto de meteoritos 
(Rivera, 1999). Sin embargo, desde hace 
varias décadas se ha incluido al hombre 
como una de las causas de cambio climá- 
tico. Tal cambio le es atribuido directa o 
indirectamente a las actividades humanas 
que alteran la composición de la atmós- 
fera, que agregadas a la variabilidad 
climática natural observada en periodos 
comparables de tiempo, han creado un 
fuerte consenso científico de que el clima 
global se verá alterado significativamente 
en este siglo, como resultado del 
aumento de la concentración de Gases 
de Efecto Invernadero (gei), tales como el 
dióxido de carbono (C0 2 ), el metano 
(CH 4 ), los óxidos de nitrógeno (N 2 0) y los 
clorofluorocarbonos (CFCs) (ipcc, 2001). 
La concentración de gases de efecto 
invernadero en la atmósfera continúa 
incrementándose, provocando en general 
un efecto de calentamiento. Aunque en 
algunas regiones hay también un incre- 
mento en la concentración de aerosoles, 
éstos tienen un efecto opuesto en el 
balance de la radiación, esto es, tienden a 
enfriar. Los cambios en la concentración 
de gei y aerosoles son los que causan 
cambios regionales y globales en la 
temperatura, precipitación y otras varia- 
bles climáticas. Esto implica cambios 
globales en la humedad del suelo, incre- 
mentos en el nivel del mar y la proyección 
de eventos más extremosos en cuanto a 
temperaturas elevadas, inundaciones y 
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sequías en una gran mayoría de locali- 
dades (Houghton et al., 1996). El ipcc 
(Intergovernmental Panel on Climate 
Change) establece que las concentra- 
ciones de C0 2 durante el periodo 1750- 
2000 aumentaron en un 31%, debido 
principalmente a la quema de combusti- 
bles fósiles y al cambio del uso del suelo 
(ipcc, 2002). 

A consecuencia de la ocurrencia de 
eventos extremos, centros de investiga- 
ción y universidades de todo el mundo se 
han dado a la tarea de crear herramientas 
para la simulación de los impactos que los 
cambios en el clima podrían ocasionar a 
los ecosistemas y por consecuencia a las 
actividades humanas. Para estimar las 
condiciones futuras del clima con base en 
las nuevas concentraciones de los gases 
de efecto invernadero se han desarro- 
llado diferentes modelos que simulan los 
cambios en las variables climáticas como 
por ejemplo la temperatura, precipitación, 
evaporación, etc. (Magaña et al., 2000). 

A partir de la aplicación de los 
modelos de cambio climático se obtienen 
las variables base para realizar simula- 
ciones a futuro del comportamiento de 
especies, de ecosistemas o bien de 
sistemas de producción agrícola, gana- 
dero o silvícola como es el caso de las 
investigaciones realizadas por Bakkenes 
et al. (2006); Berry et al. (2002); Yáñez et 
al. (1998); Olfert y Weiss (2006); Herr- 
mann et al. (2005), así como Hitz y Smith 
(2004). 

Como se observa, es importante 
profundizar en estudios sobre los 
impactos que un posible cambio en las 
condiciones globales climáticas tendrían 
sobre los ecosistemas del país y de una 
manera más regional, con motivo de 
apoyar los criterios de planeación que 
actualmente se llevan a cabo en 
programas nacionales tales como el 
Programa Nacional de Reforestación 


32 


Distribución del cedro rojo (Cedrela odorata L.) en el estado de Hidalgo 


(pronare), Programa de Desarrollo 
Forestal (prodefor) y el Programa de 
Plantaciones Forestales Maderables 
(prodeplan), entre otros. 

OBJETIVOS 

1 . Elaborar cartografía actualizada de las 
variables climáticas, temperatura y 
precipitación del estado de Hidalgo con 
datos de 1961 a 1990. 

2. Realizar un balance de humedad en el 
estado de Hidalgo bajo condiciones 
actuales y en escenarios de cambio 
climático con la aplicación de dos 
modelos, el norteamericano (gfdl-R3o) 
y el inglés (FladCM3), para evaluar los 
posibles cambios en la disponibilidad 
de humedad en el suelo que afecta la 
distribución del cedro rojo. 

3. Comparar la distribución actual del 
cedro rojo ( Cedrela odorata L.), su 
distribución potencial actual, así como 
la distribución potencial bajo escena- 
rios de cambio climático con la aplica- 
ción de los dos modelos señalados. 


METODOLOGÍA 
Área de estudio 

El estado de Hidalgo es una de las 31 
entidades federativas de la República 
Mexicana, se localiza en la zona centro 
del país, entre las coordenadas máximas: 
19° 36' 36” y 21° 24' 36” de latitud norte 
y 97° 57' y 99° 50' 24” de longitud oeste. 
Presenta un amplio rango altitudinal que 
va de los 100 hasta los 3400 metros 
sobre el nivel del mar. Con base en el 
inventario nacional forestal del año 2001 
(semarnat-unam, 2001), la agricultura es 
el principal uso de la tierra ocupando 47% 
de la superficie estatal, seguido por el 
pastizal inducido con 15%. Se presenta 
una gran diversidad de ecosistemas 


siendo los bosques templados los de 
mayor distribución con el 26% del terri- 
torio estatal. 

Delimitación de áreas de influencia 
climática 

Para la búsqueda de información climá- 
tica inicialmente se ubicaron las esta- 
ciones climatológicas que se encuen- 
tran dentro de los límites y en la peri- 
feria del estado de Hidalgo. Posterior- 
mente se utilizaron las normales clima- 
tológicas publicadas por el Servicio 
Meteorológico Nacional (smn, 2000) 
para la extracción de los datos prome- 
dios mensuales y anual de precipita- 
ción y temperatura para el periodo de 
1961-1990. 

Para el trazo de las isoyetas se 
utilizó la metodología propuesta por 
Gómez et al. (2006), en la cual se emplea 
como información base a espaciomapas 
y/o imágenes de satélite, modelos digi- 
tales de elevación, cartas topográficas en 
formato digital y los datos de precipitación 
de las estaciones meteorológicas. Utili- 
zando el software de sistemas de infor- 
mación geográfica ArcView v.3.2 (esri, 
1996) se sobrepusieron estos productos 
cartográficos y se trazaron las isoyetas 
tomando en cuenta la dinámica de circu- 
lación de los vientos, el efecto que tienen 
los sistemas orográficos en la dinámica 
de movimiento de los sistemas de vientos 
y su impacto en la precipitación, así como 
los diversos fenómenos meteorológicos 
que afectan al estado de Hidalgo como 
nortes y huracanes (Vidal, 2005). Se esta- 
blecieron analogías entre la cantidad de 
precipitación y la distribución de la vege- 
tación. Lo anterior para conformar el 
mapa de precipitación media anual y las 
precipitaciones mensuales. 

Para obtener las isotermas se dividió 
al estado en dos zonas determinadas por el 
efecto de las barreras orográficas en la 
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humedad dominante en el aire: barlovento 
y sotavento. Así, la vertiente que ve hacia el 
Golfo de México quedó dentro del área 
denominada barlovento y la vertiente inte- 
rior del estado en el sotavento. La línea 
divisoria se estableció en la isoyeta anual 
de 1 000 mm. Para cada zona se generaron 
modelos de regresión lineal simple para 
establecer el rango de variación de la 
temperatura con la altura, obteniéndose 
doce modelos mensuales y uno anual para 
cada área delimitada (Gómez et al., 2006). 
Aplicando los modelos obtenidos se 
trazaron las isolíneas de temperatura por 
cada rango altitudinal, conformando así los 
doce mapas de temperatura mensual y el 
anual en cada área. Se unieron posterior- 
mente, para hacer las correcciones en las 
áreas de confluencia cuando se presen- 
taron diferencias de temperatura. 

Finalmente, se obtuvieron las áreas 
de influencia climática (aic) al sobreponer, 
mediante el programa de sistemas de 
información geográfica ArcView v.3.2 
(esri, 1996) las isoyetas con las 
isotermas anuales. A cada una de las aic 
resultantes se le atribuyó su correspon- 
diente base de datos de temperatura y 
precipitación mensual. 

Aplicación de las razones de cambio 
estimadas con los Modelos de Cambio 
Climático 

Las AIC obtenidas se sobrepusieron con las 
zonas de cambio de temperatura y precipi- 
tación obtenidas por el Canadian Institute 
for Climate Studies (cíes, 2006) para el 
modelo de cambio climático gfdl-R30 
desarrollado por el Geophysical Fluid Dyna- 
mics Laboratory, de los Estados Unidos de 
Norteamérica y el modelo HadCM3 
desarrollado por el Hadley Centre for 
Climate Prediction and Research ubicado 
en el Reino Unido. Los modelos se 
evaluaron para los años 2020 y 2050. 
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Requerimientos climáticos que deter- 
minan la distribución del cedro rojo 
( Cedrela odorata L) 

La matriz de los requerimientos climáticos 
para la distribución de Cedrela odorata L. 
fue una compilación de información de 
diferentes fuentes oficiales, siendo la prin- 
cipal la ficha técnica del SIRE (Sistema 
de Información para la Reforestación) 
distribuida por la Comisión Nacional 
Forestal (conafor, 2000) así como de la 
base de datos de Especies Nativas 
Valiosas para la Reforestación (Vázquez 
et al., 1999) distribuida por la Comisión 
Nacional para el Conocimiento y Uso de 
la Biodiversidad (conabio) para esta 
especie. 

Distribución actual y potencial del 
cedro rojo 

La distribución actual del cedro rojo se 
estimó con base en los reportes de la 
conafor (2000). Estos indican que esta 
especie se localiza primordialmente en 
ecosistemas húmedos, en específico, en 
el bosque tropical perennifolio, el bosque 
tropical subcaducifolio, el bosque tropical 
caducifolio y el bosque mesófilo de 
montaña. 

La superficie en la que actualmente se 
distribuye se estimó a partir de la infor- 
mación del inventario nacional forestal 
del año 2001 (semarnat-unam, 2001) en 
donde se presentan los tipos de vege- 
tación con los que se asocia al cedro 
rojo. 

Para delimitar la distribución poten- 
cial de esta especie, se determinaron las 
condiciones climáticas en las que es 
posible que prospere considerando 
cuatro niveles de aptitud: 1) Apto, 2) 
Moderadamente Apto, 3) Marginalmente 
Apto y 4) No Apto. 
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Simulación de la distribución del 
cedro rojo con cambio climático 

Para delimitar las posibles zonas en 
donde se distribuirá el cedro rojo, se iden- 
tificaron los valores de temperatura y 
precipitación modificados por los dos 
modelos de cambio climático de las áreas 
de influencia que correspondieran con los 
rangos de distribución encontrados. 

Considerando la distribución actual y 
potencial de Cedrela odorata L. y bajo esce- 
narios de cambio climático, se obtuvieron 
seis escenarios de comportamiento del 
cedro rojo: 1) el actual o base usando la 
información reportada por el inventario 
nacional forestal (semarnat-unam, 2001); 2) 
el potencial de acuerdo a las condiciones 
actuales; 3) y 4) el potencial futuro de 
acuerdo al modelo norteamericano de 
cambio climático a los años 2020 y 2050 y 5) 
y 6) la posible distribución de acuerdo al 
modelo inglés de cambio climático para los 
años 2020 y 2050. 

Balance de humedad en el estado de 
Hidalgo 

Para evaluar los cambios en la distribución 
de Cedrela odorata L. por efecto del cambio 
climático se realizó un balance de humedad 
en el estado de Hidalgo para los escenarios 
descritos (2, 3, 4, 5 y 6) y así identificar los 
posibles cambios en la cantidad de 
humedad disponible en el suelo. El balance 
se realizó utilizando la metodología de 
Thornthwaite modificada versión III (Monte- 
rroso y Gómez, 2003). Esta metodología 
fue creada en la década de los 90 por 
países de la Unión Europea para contar con 
mayor información sobre el impacto del 
cambio climático en sus regiones. También 
fue incorporada por la Convención de las 
Naciones Unidas de Lucha Contra la Deser- 
tificación en Países Afectados por Sequía 
Grave y Desertificación, para delimitar las 
zonas áridas y semiáridas bajo un mismo 
criterio (Monterroso y Gómez, 2004). 


El sistema incluye la aplicación del 
concepto del balance hídrico. Este 
consiste en la relación entre el agua apor- 
tada por la precipitación y la requerida por 
la vegetación en sus funciones de evapo- 
ración y transpiración (Dunne y Leopold, 
1978). El sistema de Thornthwaite consi- 
dera un nuevo concepto en la clasifica- 
ción climática, como es la eficiencia de 
los parámetros meteorológicos de precipi- 
tación y temperatura. La evaluación del 
primero se determina a través de un 
balance de humedad con la finalidad de 
determinar qué tan seco o qué tan 
húmedo es el clima en relación con el 
crecimiento de las plantas. En la evalua- 
ción del segundo parámetro se estima 
qué tan caliente o qué tan frío es el clima 
para el crecimiento de las plantas. El 
sistema considera además la capacidad 
de almacenamiento de humedad del 
suelo, que es la cantidad máxima de agua 
disponible en las capas del suelo 
ocupadas por las raíces y que puede ser 
tomada por las mismas cuando la planta 
lo requiera para realizar sus funciones. En 
su estimación se considera la cantidad de 
humedad a partir de la situación del punto 
de marchitez permanente hasta capa- 
cidad de campo. En la Tabla 1 se 
presenta la descripción de los tipos de 
climas de la metodología de Thornthwaite 
modificada, resultantes del balance de 
humedad. 

RESULTADOS 

Delimitación de áreas de influencia 
climática. 

Para el trazo de las isotermas, se esti- 
maron inicialmente los modelos de regre- 
sión lineal simple que se presentan en la 
Tabla 2 en donde se observa que los 
valores del coeficiente de regresión (R 2 ) 
son altos, mayores de 0,9 para la zona de 
barlovento y mayores de 0,8 para sota- 
vento. 
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Tabla 1. Descripción de los tipos de climas usados en la metodología. Thornthwaite 

modificada versión III 


CLIMA 

PERIODO DE 

DÉFICIT DE 

DEMASÍAS DE 

LLUVIAS 

HUMEDAD 

AGUA 

Árido (E) 

< 3 meses 

Muy severo 

No existen 

Semiárido Seco (DI) 

3 meses 

Severo 

Mínimas 

Semiárido Moderado (D2) 

3 a 4 meses 

Severo 

Mínimas 

Semiárido Húmedo (D3) 

3 a 5 meses 

No tan severo 

Mínimas 

Subhúmedo Seco (C1) 

4 a 5 meses 

Reducido 

Moderadas 

Subhúmedo Húmedo (C2) 

4 a 6 meses 

Mínimo 

Mayores 

Húmedo (B) 

> 5 meses 

Mínimo 

Mayores 

Per Húmedo (A) 

> 6 meses 

No existe 

Mayores 


Fuente: Monterroso y Gómez, 2003 


Tabla 2. Modelos de regresión lineal simple para temperatura (T) y la altura (h) en la 
zona de barlovento y sotavento en el estado de Hidalgo 


MES 

BARLOVENTO 


SOTAVENTO 


MODELO 

R 2 

MODELO 

R 2 


ENERO 

T = -0,0036h + 18,86 

0,97 

T = -0,0049h + 22,90 

0,82 

FEBRERO 

T = -0,0038h + 20,12 

0,96 

T = -0,0053h + 24,78 

0,83 

MARZO 

T = -0,0041 h + 23,56 

0,94 

T = -0,006h + 29,07 

0,87 

ABRIL 

T = -0,0049h + 26,89 

0,97 

T = -0,0063h + 31 ,40 

0,87 

MAYO 

T = -0,0055h + 28,74 

0,97 

T = -0,0068h + 33,38 

0,88 

JUNIO 

T = -0,0063h + 29,18 

0,98 

T = -0,0067h + 32,21 

0,86 

JULIO 

T = -0,0062h + 28,38 

0,97 

T = -0,0067h + 31 ,33 

0,86 

AGOSTO 

T = -0,0063h + 28,46 

0,98 

T = -0,0068h + 31,72 

0,86 

SEPTIEMBRE 

T = -0,006h + 27,76 

0,98 

T = -0,0063h + 30,18 

0,86 

OCTUBRE 

T = -0,0056h + 25,45 

0,98 

T = -0,0058h + 27,80 

0,84 

NOVIEMBRE 

T = -0,0049h + 22,84 

0,97 

T = -0,0055h + 25,71 

0,85 

DICIEMBRE 

T = -0,0042h + 20,69 

0,95 

T = -0,0052h + 24,06 

0,85 

ANUAL 

T = -0,0051 h + 25,07 

0,97 

T = -0,006h + 28,72 

0,89 


En el mapa que contiene los rangos de 
altura se digitalizó la información obtenida 
de los rangos altitudinales de variación de 
la temperatura de los modelos de regre- 
sión lineal. Con esto se obtuvieron gráfi- 
camente las isotermas mensuales y 
anuales para todo el estado (Figura la). 

Las regiones donde se presentan las 
temperaturas más altas, son la zona de 


barlovento asociada a la planicie costera 
y zona baja de la Sierra Madre Oriental, 
así como partes bajas de las cañadas de 
los ríos que desembocan en el Golfo de 
México, con valores medios anuales de 
23 a 27°C. En las partes bajas de la 
barranca de Metztitlán y la parte baja del 
Valle del Mezquital en la zona de Ixmi- 
quilpan, los rangos de temperatura media 
anual son también altos y van de 17 a 
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23°C. Por otro lado, las áreas con mayor 
elevación como la caldera de Huichapan, 
las montañas de Ixmiquilpan-Zimapan y 
la región sureste del estado son las que 
presentan las temperaturas más bajas 
con valores medios anuales de 9 a 11°C. 

La Figura Ib muestra que la parte 
más húmeda del estado se presenta en 
las regiones noreste y este a lo largo de la 
zona de barlovento de la Sierra Madre 
Oriental con valores medios anuales de 
1900 a 2300 mm. La parte más seca se 
ubica en el centro y sur de la entidad que 
comprende las áreas del Valle del 
Mezquital y Pachuca, cuyos valores de 
precipitación media anual están entre los 
300 y 500 mm. En la mesa central, la 
precipitación vuelve a aumentar hacia los 
límites con el Estado de México, con 
valores de 500 a 700 mm anuales e 
incluso con algunas áreas con precipita- 
ción media anual de hasta 800 mm. La 
barranca de Metztitlán es otra zona con 
valores bajos de precipitación, con 
valores de 300 a 500 mm aumentando 


considerablemente la precipitación al 
ascender los vientos en las elevaciones 
montañosas de los municipios de Mineral 
del Chico y Mineral del Monte, alcan- 
zando valores de más de 1000 mm 
(Tinoco, 2005). 

Con la sobreposición de la carto- 
grafía resultante de las isoyetas e 
isotermas se obtuvo un mapa con 7981 
áreas de influencia climática. Cada aic se 
completó con su respectiva información 
de temperatura y precipitación mensual y 
anual, que posteriormente sirvió para 
ubicar las zonas de distribución potencial 
del cedro rojo. 

Aplicación de modelos de cambio 
climático 

Se sobrepuso el mapa de las aic a la 
cartografía de cambios de temperatura y 
precipitación definidos por el Canadian 
Institute for Climate Studies (cíes, 2006) 
para los modelos gfdl-R 30 y HadCM3, 
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Figura 1. Mapas de temperatura y precipitación media anual del estado de Hidalgo: 
a = temperatura media anual; b = precipitación media anual. 
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obteniéndose con esto las variables 
climáticas para los años 2020 y 2050 para 
el caso de los dos modelos utilizados. 

Al aplicar los cambios a la tempera- 
tura y precipitación, se observa que de los 
dos modelos usados, el gfdl-r30 es el 
que pronostica un aumento en las dos 
variables climáticas (Figura 2). El valor 
promedio de incremento en toda la 
entidad para la precipitación media anual 
aplicando este modelo para los años 
2020 y 2050 es de 100 mm y 120 mm, 
respectivamente, esto comparado con el 
escenario base que es de 1300 mm 
anuales, el mayor aumento para los dos 
periodos se presenta en el mes de julio. 
El mes de octubre, por su parte, es en 
donde se observa la mayor disminución 
de la lluvia bajo este modelo. En cuanto a 
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temperatura, los dos años simulados 
establecen un aumento en todos los 
meses, siendo el escenario para el 2050 
en el que mayor incremento se observa, 
con un promedio anual de 1,7°C para 
todo el estado de Hidalgo a partir del valor 
del escenario base que es de una tempe- 
ratura media anual de 17,9°C para todo el 
estado. 

El modelo HadCM3 pronostica un 
aumento mayor de temperatura en toda la 
entidad con 1,1 °C para el escenario del 
año 2020 y 2,5°C para el año 2050 
(Figura 2). En cuanto a la precipitación se 
refiere, el modelo inglés establece una 
merma en la cantidad de lluvia promedio 
anual para los dos periodos con relación 
al escenario base de 137,8 mm para el 
año 2020 y 131,6 mm para el 2050. 



Figura 2. Temperatura y precipitación promedio anual para las condiciones actuales y 

los escenarios de cambio climático 
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Distribución del cedro rojo (Cedrela odorata L.) en el estado de Hidalgo 


Balance de humedad en el estado de 
Hidalgo 

En la Figura 3 se muestra la proporción 
de los tipos climáticos en los cinco esce- 
narios evaluados. En el escenario actual 
(1961-1990), el tipo de clima que predo- 
mina es el Semiárido húmedo (D3) 
ocupando el 32,9% del territorio estatal, le 
sigue por cobertura el clima Húmedo (B) 
con el 19,6% y es seguido por el clima 
Per húmedo (A) con el 18,9%. Los climas 
con menor cobertura fueron el Semiárido 
moderado (D2), Subhúmedo seco (C1) y 
el Subhúmedo húmedo (C2) con una 
cobertura del 14,9%, 9,6% y 4,2%, 
respectivamente. 

Los periodos simulados de cambio 
climático (2020 y 2050) con el modelo 
gfdl-R30 no presentan mucha diferencia 
con respecto a los climas del escenario 
base. Para el año 2020, el clima con 


mayor superficie es también el Semiárido 
húmedo con 32,7% de la superficie total 
del estado, seguido por el clima Húmedo 
con 21,7% y el clima Per húmedo con 
13%. Los climas C2, C1 y D2 ocupan 
cerca del 10% cada uno de ellos. Para el 
año 2050 la proporción es muy parecida 
al anterior escenario, donde el clima 
Semiárido húmedo ocupa la mayor super- 
ficie con 30,2%, seguido del clima 
Húmedo con 22,1%, el Per húmedo, 
Subhúmedo húmedo, Subhúmedo seco y 
el Semiárido moderado ocupan cerca del 
10% cada uno. 

Los resultados al aplicar el modelo 
HadCM3 muestran una disminución de 
los climas Per húmedo (A) con respecto 
al escenario base y un incremento 
sustancial en el clima Semiárido mode- 
rado (D2). Para el año 2020 se pronos- 
tica que el clima Per húmedo ocupará el 
4,1% del estado y el clima Húmedo el 



A = Per húmedo £2 -Subhúmedo húmedo D3 = Semiárido húmedo DI = Semiáridoüeco 

B - Húmedo C1 4 Subhúmedo seco D2 - Semiárido moderado 


Figura 3. Distribución de los climas estimados con la metodología de Thornthwaite 
modificada versión III, bajo condiciones actuales y escenarios de cambio climático 
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20,4%. El clima Subhúmedo húmedo se 
pronostica en un 8,4% del territorio 
estatal, el Subhúmedo seco en el 12,6%, 
el Semiárido húmedo en el 26,9% siendo 
también el de más amplia distribución, 
seguido del Semiárido moderado con 
24,4%. En el escenario del año 2050, el 
clima A se pronostica en el 1,9% de la 
superficie estatal y el clima B se obtuvo 
para el 20,6%. Los climas C2 y C1 dismi- 
nuyen su cobertura en comparación con 
el periodo anterior al ubicarse en el 7,4% 
y 8,1% de la superficie, respectivamente. 
Por otro lado, los climas D3 y D2 
aumentan su superficie con respecto al 
escenario anterior, estimándose que 
ocuparán el 30,4% y 27,7%, respectiva- 
mente. 

En las simulaciones bajo el modelo 
inglés resalta la aparición del tipo de clima 
Semiárido seco (DI), ya que en los otros 
escenarios no se registra. Para el año 
2020 este modelo establece que dicho 
clima se distribuirá en el 3,1% del estado 
y para el año 2050 en el 3,9% de la super- 
ficie estatal. 
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Distribución actual de Cedrela odorata L. 

En la Tabla 3 se muestra la compilación 
de los requerimientos climáticos del cedro 
rojo que se obtuvieron de las fuentes de 
información consultadas (conafor, 2000; 
Vázquez et al., 1999). 

La distribución actual del cedro rojo 
se presenta en la Figura 4a, misma que 
se obtuvo de la información del inven- 
tario nacional forestal (semarnat-unam, 
2001) relacionándola con la distribución 
de los tipos de vegetación con los que, 
de acuerdo a la conafor (2000), se 
asocia esta especie. Estos tipos de 
vegetación fueron el bosque tropical 
perennifolio, el bosque tropical subcadu- 
cifolio, el bosque tropical caducifolio y el 
bosque mesófilo de montaña. La infor- 
mación obtenida muestra que el cedro 
rojo actualmente se encuentra distri- 
buido en áreas fragmentadas a lo largo 
de la vertiente del Golfo de México, 
cubriendo una superficie de 2040,3 km 2 
lo que equivale al 9,8% del territorio total 
del estado de Hidalgo. 


Tabla 3. Características climáticas y nivel de aptitud en las que crece el cedro rojo 

( Cedrela odorata L.) 


VARIABLE 


NIVEL DE APTITUD* 


1 

2 

3 

4 

Precipitación media 
Anual (mm) 

2500-4000 

4000-5000 

2000-2500 

>5000 

1200-2000 

<1200 

Duración del periodo 
Seco (meses) 

0-3 

3-5 

5-6 

> 6 

Temperatura media 
Anual (°C) 

20-26 

26-28 

18-20 

28-30 

16-18 

> 30 
< 16 

Temperatura mínima 
Promedio (°C) 

16-22 

22-24 

14-16 

24-26 

11-14 

> 26 
< 11 

Temperatura máxima 
Promedio (°C) 

24-30 

30-34 

22-24 

34-38 

20-22 

> 38 
< 20 


Donde: 1 = Apto, 2 = Moderadamente Apto, 3 = Marginalmente Apto y 4 = No Apto 
*Adaptado por los autores a partir de: CONAFOR (2000) y Vázquez et al. (1999). 
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Distribución del cedro rojo (Cedrela odorata L.) en el estado de Hidalgo 



Figura 4. a= Distribución actual de los tipos de vegetación con los que se asocia 
Cedrela odorata L. y b= niveles de aptitud de distribución potencial actual 


Distribución potencial actual de 
Cedrela odorata L. y bajo escenarios 
de cambio climático 

Aplicando los requerimientos climáticos 
del cedro rojo a la información actual 
de las áreas de influencia climática, se 
determinó la distribución potencial 
actual (Figura 4b), en donde se aprecia 
que el área con potencial climático para 
el desarrollo de la especie se 
encuentra ubicada en la zona de barlo- 
vento. 

En la tabla 4 se presentan las superficies 
cubiertas por nivel de aptitud a partir de las 
variables climáticas para los cinco escena- 
rios de acuerdo a los requerimientos de la 
especie expresados en la Tabla 3. El factor 
más restrictivo en todos los escenarios fue 
la precipitación media anual, no presentán- 
dose ninguna área con el nivel más alto de 
aptitud (Apto) para todos los escenarios 
estudiados. En los escenarios del modelo 
GFDL-R30 para este factor, en general 
aumenta la superficie con algún nivel de 


aptitud, presentándose un incremento con 
respecto al potencial actual de 472,84 y 
628,64 km 2 para los escenarios de 2020 y 
2050, respectivamente. Dentro de los 
niveles de aptitud, hay un incremento en la 
superficie para el nivel Moderadamente 
apto y una disminución para el nivel Margi- 
nalmente apto. En cambio para el modelo 
HadCM3, hay una disminución de la super- 
ficie con cierto grado de aptitud por este 
factor, con respecto al potencial actual, de 
406,78 y 367,28 km 2 para los escenarios de 
2020 y 2050, respectivamente. La mayor 
parte de la superficie se ubica en la clase 
Marginalmente apto, y la disminución en la 
clase Moderadamente apto tiene valores de 
cero para el escenario del 2020 y de tan 
sólo 271,39 km 2 para el año 2050. 

Para las variables asociadas a 
temperatura, en ambos modelos se 
presenta un aumento en la superficie con 
cierto grado de aptitud. Para los escena- 
rios del GFDL-R30, el incremento en la 
superficie con aptitud con respecto al 
potencial actual del escenario base fue de 
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Tabla 4. Superficie estimada por nivel de aptitud de Cedrela odorata L. de acuerdo a las 
variables climáticas para los cinco escenarios evaluados 


Escenario 


Superficie por Nivel de Aptitud (km 2 ) 

Apto 

Moderadamente 

Apto 

Marginalmente 

Apto 

No Apto 

Precipitación 

Potencial Actual 

0,00 

2 565,24 

4 202,83 

14 061,25 

GFDL-R30 al 2020 

0,00 

3 959,92 

3 115,88 

13 588,40 

GFDL-R30 al 2050 

0,00 

4 000,53 

3 231,06 

13 432,61 

HadCM3 al 2020 

0,00 

0,00 

6 196,18 

14 468,03 

HadCM3 al 2050 

0,00 

271,39 

5 964,32 

14 428,50 

Temperatura 

Potencial Actual 

4 573,10 

2 314,82 

4 334,73 

9 606,67 

GFDL-R30 al 2020 

5 081,62 

4 192,26 

6 470,19 

4 920,14 

GFDL-R30 al 2050 

1 495,63 

3 036,59 

13 438,01 

2 694,09 

HadCM3 al 2020 

5 081,62 

4 192,26 

6 470,19 

4 920,14 

HadCM3 al 2050 

7 456,29 

7 822,09 

4 576,06 

811,76 

Periodo seco 

Potencial Actual 

7 313,99 

1 457,16 

9 993,17 

2 064,99 

GFDL-R30 al 2020 

6 690,21 

8 532,11 

3 391,65 

2 050,24 

GFDL-R30 al 2050 

7 092,55 

5 838,55 

5 044,28 

2 688,83 

HadCM3 al 2020 

6 494,89 

2 331,32 

5 536,62 

6 301,38 

HadCM3 al 2050 

6 307,68 

5 646,68 

4 034,62 

4 675,23 


4686,53 y 6912,58 km 2 para el año 2020 y 
2050, respectivamente. Para el modelo 
HadCM3, el incremento de la superficie 
con aptitud fue de 4686,53 y 8794,94 km 2 
para el año 2020 y 2050, respectivamente. 

Para la variable duración del periodo 
seco que se estimó a través del balance de 
humedad entre la precipitación y la 
demanda evapotranspirativa, según el 
modelo gfdl-R30 no se tiene gran variación 
de la superficie con cierto grado de aptitud 
con respecto al escenario base. Hubo un 
ligero incremento de 14,59 km 2 para el año 
2020 y un decremento de 623,84 km 2 para 
el 2050. En cambio, para el modelo 
HadCM3, para ambos escenarios se esti- 
maron decrementos significativos de 
4236,39 y 2610,24 km 2 para el año 2020 y 
2050, respectivamente. 


La variable con mayor restricción 
determinó los niveles de aptitud, 
mismos que se presentan resumidos 
en la tabla 5. La superficie potencial 
estimada bajo el escenario actual con 
categoría Moderadamente Apta se 
ubica en la parte noreste del estado, 
sobre las laderas intermedias de los 
flancos orientales de la Sierra Madre 
Oriental así como una pequeña porción 
en la zona sureste (Figura 4b), 
cubriendo 2169,80 km 2 , que corres- 
ponden al 10,5% del total de la super- 
ficie estatal. El nivel de aptitud Margi- 
nalmente Apto ocupa una superficie de 
4102,73 km 2 , que corresponde al 
19,86% del total de la superficie. El 
69,65% de la superficie estatal corres- 
ponde al nivel No apto. 
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Tabla 5. Superficie estimada por nivel de aptitud de Cedrela odorata L. para los cinco 

escenarios evaluados 


Superficie por Nivel de Aptitud (km 2 ) 

escenario 

Apto 

Moderadamente 

Apto 

Marginalmente 

Apto 

No Apto 

Potencial Actual 

0,00 

2 169,80 

4 102,73 

14 391,58 

GFDL-R30 al 2020 

0,00 

3 708,86 

3 195,83 

13 759,52 

GFDL-R30 al 2050 

0,00 

2 304,20 

4 879,91 

13 480,09 

HadCM3 al 2020 

0,00 

0,00 

6 091,38 

14 572,83 

HadCM3 al 2050 

0,00 

270,50 

5 962,97 

14 430,74 


En la figura 5 se muestran los resultados 
del comportamiento del cedro rojo apli- 
cando los dos modelos de cambio climá- 
tico. Para el modelo gfdl-R 30 , a nivel 
general se observa una disminución del 
nivel No Apto para ambos periodos (2020 
y 2050) siendo de 66,59% y 65,23%, 
respectivamente, en comparación del 
69,65% que se tiene para el periodo poten- 
cial actual (1961-1990). En el nivel de 
Moderadamente apto, en ambos escena- 
rios de cambio climático se tuvo un 


aumento, pasando de 10,5% en el esce- 
nario base a 17,95% para el año 2020 y a 
11,15% en el 2050. El Nivel Ligeramente 
apto disminuye para el escenario del 2020 
con respecto al potencial actual pasando 
de 19,86% de la superficie a 15,47%. Sin 
embargo, para el año 2050 aumenta a un 
23,62% del total de la superficie. La distri- 
bución espacial de los diferentes niveles 
de aptitud con este modelo se muestra en 
la Figura 6a y 6b para el año 2020 y 2050, 
respectivamente. 



BASE (1 961-1990) 


GFDL-R3Ü 2020 


GFDL-R30 2050 

Escenario 


HADCM3 2020 


HADCM3 2050 


MA = Moderadamente Apto 


LA = Marginalmente Apto 


NA = No Apto 


Figura 5. Cambios en los niveles de aptitud de Cedrela odorata L. bajo condiciones 

actuales y con escenarios de cambio climático 
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Figura 6. Distribución potencial de Cedrela odorata L. bajo escenarios de cambio 
climático estimados con el modelo GFDL-R30: (a) 2020; (b) 2050 


El comportamiento del nivel de 
aptitud con el modelo FladCM3 (Figura 6), 
muestra un pequeño aumento en la 
superficie catalogada como No apta, con 
un porcentaje de la superficie estatal de 
70,52% para el año 2020 y de 69,83% en 
el 2050. En la clase de Moderadamente 
apto, se da una gran disminución, 
pasando de 10,5% para el escenario 


potencial actual a cero en el año 2020 y 
1,31% en el 2050. El Nivel Marginalmente 
apto, aumenta de 19,86% para el esce- 
nario base a 29,48% y 28,86% para los 
años 2020 y 2050, respectivamente. La 
distribución espacial de los diferentes 
niveles de aptitud con este modelo se 
muestra en la Figura 7a y 7b para el año 
2020 y 2050, respectivamente. 



Figura 7. Distribución potencial de Cedrela odorata L. bajo escenarios de cambio 
climático estimados con el modelo HadCM3: (a) 2020; (b) 2050 
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DISCUSIÓN 

Delimitación de áreas de influencia 
climática 

Los modelos de regresión lineal para las 
dos zonas presentan variaciones en el 
rango de incremento de la temperatura al 
aumentar la altura. La zona de barlo- 
vento, para todos los meses y el modelo 
anual, tiene rangos de variación altotér- 
mica menores que la de sotavento. Esto 
se asocia al mayor contenido de vapor de 
agua del aire en esta zona, que al 
ascender sobre la pendiente, se presenta 
la condensación y la consecuente libera- 
ción de calor latente al pasar el agua de 
vapor a líquido, mismo que calienta el aire 
y reduce el rango de enfriamiento 
(Arhens, 2003). Los coeficientes de regre- 
sión mayores para la zona de barlovento 
se asocian con el mayor rango de altura 
en que están localizadas las estaciones 
meteorológicas que se incluyen en los 
modelos (Gómez et al., 2006). 

El trazo de las isoyetas mensuales 
muestra que la distribución de la precipi- 
tación en todo el estado es primordial- 
mente en verano como resultado del 
efecto de los vientos alisios provenientes 
del Golfo de México, así como también de 
los huracanes que se forman en el Caribe 
y siguen su curso por el Golfo entre los 
meses de Julio a Octubre. Durante el 
invierno las lluvias son influenciadas por 
los frentes fríos llamados “nortes”, 
mismos que están asociados al compo- 
nente invernal del Monzón de Norteamé- 
rica, los cuales también aportan lluvia de 
forma considerable, principalmente en la 
zona de barlovento, aunque en menor 
cantidad que los sistemas de verano 
(Vidal, 2005). 

Las zonas secas comprenden las 
áreas donde los vientos que pasan la 
barrera montañosa de la Sierra Madre 
Oriental al descender se comprimen, lo 


que ocasiona sombra de lluvia, y por lo 
tanto disminuye la precipitación. La preci- 
pitación aumenta de nuevo conforme las 
masas de vientos se retiran de las áreas 
de compresión y al desplazarse sobre 
valles altos donde la convección local es 
un elemento que aumenta el ascenso del 
aire y la formación de nubes (Mosiño y 
García, 1974). El aumento considerable 
de la precipitación en el macizo monta- 
ñoso de Mineral del Chico y Mineral del 
Monte, se asocia con el ascenso forzado 
del aire sobre las laderas orientales de 
dichos sistemas montañosos y con la 
formación de nubes al enfriarse el aire por 
debajo del punto de condensación del 
vapor de agua y la formación de nubes, 
sobre todo cuando el gradiente de 
presión es considerable entre la zona del 
Golfo de México y el altiplano central que 
establece el flujo de los vientos alisios 
(Tinoco, 2005). 

Aplicación de modelos de cambio 
climático. Ambos modelos pronostican un 
incremento de la temperatura, siendo 
mayor dicho aumento para el escenario 
del año 2050. El incremento de tempera- 
tura está ampliamente explicado por la 
acumulación de Gases de Efecto Inverna- 
dero (GEI) que, al ir aumentando su 
concentración en la atmósfera, causan 
una mayor retransmisión a la superficie 
de la radiación de onda larga irradiada 
por la superficie de la tierra, lo que 
calienta más el aire y provoca un incre- 
mento en la temperatura. Para el esce- 
nario más tardío se preve mayor acumu- 
lación de estos gases y la magnificación 
del fenómeno de calentamiento (IPCC, 
2001 ). 

El modelo GFDL-R30 pronostica un 
incremento menor de temperatura que el 
modelo HadCM3, con 0,95°C para el 
2020 y 1,7°C para el 2050 comparados 
con 1 ,1 °C y 2,5°C para estos años que se 
obtienen con el modelo inglés. Esta dife- 
rencia se asocia con las series de tiempo 
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con las que se calibraron dichos modelos. 
El modelo norteamericano considera una 
cantidad significativa de estaciones de 
zonas con alta concentración de aero- 
soles, que son comunes en la costa este 
de dicho país, lo que reduce el calenta- 
miento. En cambio, el modelo inglés 
presenta datos más globales donde estos 
elementos tienen menor peso (McGregor 
y Nieuwolt, 1998). 

Las diferencias en precipitación esti- 
madas en los escenarios de los modelos, 
indican que el modelo inglés reporta una 
disminución que se acentúa ligeramente 
al aumentar los años considerados. Esto 
se explica por un desplazamiento 
pequeño hacia el sur de la zona de 
corrientes descendentes y el menor 
avance de los centros de alta presión 
hacia el océano en la temporada caliente 
del año, lo que disminuirá la precipitación 
en las zonas aledañas al movimiento de 
estos sistemas. En cambio el modelo 
norteamericano pronostica una mayor 
evaporación sin que esta zona se vea 
afectada considerablemente por el movi- 
miento de los sistemas anteriores, lo que 
trae consigo un incremento en la precipi- 
tación (McGregor y Nieuwolt, 1998). 

Balance de humedad en el estado de 
Hidalgo 

Ambos modelos reportan un incremento 
en el déficit de humedad, con el desplaza- 
miento a tipos climáticos más secos. Este 
incremento es mayor para los resultados 
del modelo inglés, ya que en sus escena- 
rios se estiman mayores incrementos en la 
temperatura, lo que aumenta la demanda 
evapotranspirativa, que aunado a la dismi- 
nución de la precipitación, el balance de 
humedad presenta mayor déficit. En las 
estimaciones con el modelo norteameri- 
cano, aunque se preve que aumente la 
precipitación, esta no lo hace en la misma 
proporción que la demanda evapotranspi- 
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rativa, y aunque el cambio en el balance 
de humedad es menos drástico que lo 
reportado en el otro modelo, también está 
presente un ligero aumento en el déficit y 
disminución de las demasías (Monterroso 
y Gómez, 2003). 

Distribución actual de Cedrela odorata L. 

Se observa que C. odorata L. es una 
especie característica de los ecosistemas 
definidos por los grupos climáticos del 
trópico húmedo y trópico subhúmedo. 
Dentro de los cuales se llega a encontrar 
tanto en las partes bajas de los bosques 
mesófilos de montaña, como en las partes 
bajas de los bosques de pino y en los 
bosques de pino-encino, además de los 
bosques tropicales perennifolio, subperen- 
nifolio, caducifolio y subcaducifolio. 
Presenta una distribución amplia en 
cuanto a precipitación media anual. Sin 
embargo en zonas con precipitaciones 
notablemente menores a 2500 mm, la 
especie no se desarrolla favorablemente y 
presenta fustes cortos y frecuentemente 
torcidos (Vázquez et al., 1999). 

Distribución potencial actual de 
Cedrela odorata L. y bajo escenarios 
de cambio climático 

La condición climática más limitativa es la 
precipitación media anual, ya que en el 
estado de Hidalgo no rebasa los 2500 
mm. Esto implica que no se encontrarán 
áreas con el nivel de aptitud mayor 
(Apta). Sin embargo, con los datos del 
escenario base se estima que aproxima- 
damente el 30% de la superficie estatal 
tenga un cierto grado de aptitud para esta 
especie. La ubicación de estas áreas se 
presenta en la zona de barlovento, sobre 
parte de la planicie costera y las laderas 
orientales de la Sierra Madre Oriental así 
como en las cañadas que dan al Golfo de 
México. Los resultados de la aplicación 
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del modelo gfdl-r30 indican un ligero 
incremento de un poco menos de 5% en 
la superficie con cierto grado de aptitud 
para esta especie, y se asocia con el 
incremento de la altura en la que puede 
prosperar el cedro al tenerse tempera- 
turas más altas y el rango de precipitación 
y duración del periodo seco favorables 
dentro de estos niveles de aptitud. 

Los escenarios del modelo HadCM3 
reportan una disminución drástica en los 
niveles de aptitud, concentrándose toda 
la superficie en prácticamente el nivel 
más bajo (Marginalmemnte apto); aunque 
la proporción de la superficie no apta 
permanece casi sin cambio respecto al 
escenario base. Sin embargo, en las 
zonas con menor altitud, mismas que les 
corresponden las temperaturas más altas 
con el incremento estimado, la evapo- 
transpiración potencial se incrementa 
exponencialmente, determinando que en 
el balance de humedad la duración del 
periodo seco se incremente, pasando 
estas áreas a ser no aptas para la especie 
(Figura 7a y 7b). En los escenarios resul- 
tantes de este modelo también se tiene 
un desplazamiento a mayor altura de la 
cota de áreas con cierto grado de adapta- 
bilidad en la zona de barlovento. 

CONCLUSIONES 

Aunque actualmente se reporta que el 
9,8% de la superficie del estado presenta 
los tipos de vegetación con los que se 
asocia Cedrela odorata L., las estima- 
ciones a partir de las características 
climáticas del escenario base (1961- 
1990), muestran que el 30,4% del estado 
presenta algún grado de aptitud para el 
desarrollo de esta especie. Las áreas con 
aptitud se ubican fundamentalmente en la 
zona de barlovento, que comprende la 
porción de la planicie costera del Golfo, el 
flanco este de las laderas de la Sierra 
Madre Oriental y las cañadas que desem- 
bocan a la planicie costera. 


Los escenarios de cambio climático 
estimados con los modelos GFDL-R30 y 
el HadCM3 muestran una variabilidad en 
la distribución del cedro rojo. Ambos 
modelos pronostican un aumento en la 
temperatura, lo que incrementa la altura 
(sobre el nivel del mar) del área con 
aptitud en las laderas orientales de la 
Sierra Madre. En cambio, el modelo 
norteamericano pronostica un incremento 
menor de temperatura y un aumento 
ligero en la precipitación, lo que aumenta 
el área con algún nivel de aptitud con 
respecto al escenario base actual, con un 
3,1% para el año 2020 y de 4,4% para el 
año 2050. Sin embargo, este aumento es 
diferencial dentro de las clases de aptitud, 
para el primer periodo la clase Moderada- 
mente apta aumenta de 10,5% en el 
escenario base a 18,0% para el año 
2020, pero para el escenario de 2050 
disminuye a 11,2% que, aunque mayor 
que el escenario base, el incremento es 
reducido. 

Bajo el modelo HadCM3 la distribu- 
ción del cedro rojo se ve afectada. Esta 
se desplaza de la zona baja de la planicie 
costera a rangos de mayor altitud sobre 
las laderas del flanco este de la Sierra 
Madre Oriental. La superficie con nivel de 
aptitud disminuye 0,9% para el año 2020 
y 0,2% para el año 2050. Sin embargo, el 
nivel de Moderadamente apto pasa de 
10,5% en el escenario base a 0% para el 
año 2020 y a 1,3% para el año 2050, 
ubicándose prácticamente todas las 
áreas en el menor nivel de aptitud. 

Los resultados diferentes de los dos 
modelos de cambio climático resaltan lo 
complejo de estos estudios, en donde 
existe un rango amplio de incertidumbre 
que aún se tiene ante cambios en los 
componentes climáticos globales y regio- 
nales. Al ser dos escenarios que señalan 
niveles muy amplios de comportamiento 
futuro, éstos deben tomarse como un 
primer acercamiento en la evaluación de 
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los posibles impactos de un cambio climá- 
tico en la especie y su ecosistema. Por 
ello, es menester que se realicen estudios 
más detallados donde se evalúe local- 
mente a la especie. Sin embargo, los 
resultados aquí presentados señalan que 
independientemente del escenario que se 
tome, se observan posibles cambios en el 
comportamiento de la especie, por lo que 
deben incorporarse acciones inmediatas 
dentro de los programas regionales de 
manejo y conservación de la especie. 

Las regiones que deben ser sujetas 
de monitoreo son aquellas donde el cedro 
rojo podría dejar de tener las condiciones 
favorables para su desarrollo, ya que 
estas son las zonas más vulnerables para 
la especie al mostrar potenciales de 
cambio. 

En este contexto, se requieren estu- 
dios futuros a escala regional y local 
donde se involucre a campesinos, toma- 
dores de decisiones, dueños de la tierra y 
gobiernos en sus distintos niveles, para 
que en conjunto se desarrollen acciones y 
planes de adaptación ante el cambio 
climático. 
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ARTÍCULO DE INVESTIGACIÓN 


Reconstrucción del historial de incendios 
y estructura forestal en bosques mixtos 
de pino-encino en la Sierra Madre Oriental 

Reconstruction of fire history and forest structure in mixed 
pine-oak forest in the Sierra Madre Oriental 

Marco Aurelio González Tagle 1 , Luitgard Schwendenmann 2 , 
Javier Jiménez Pérez 3 y Wibke Himmelsbach 4 

RESUMEN 

Los incendios forestales en México han tenido gran relevancia después de la temporada de 
incendios en 1998 debido a su contribución directa a la deforestación, el cambio en la estructura y 
composición del bosque. Además recientemente se considera su participación directa en el incre- 
mento de bióxido de carbono en la atmósfera. Sin embargo, es muy escaso el conocimiento que se 
posee acerca de los efectos de los incendios sobre la estructura, composición y dinámica de los 
bosques de la Sierra Madre Oriental (SMO). En este trabajo, se presenta una cronosecuencia de 
incendios forestales (1 34 años) determinada a través del análisis de anillos de crecimiento en cohortes 
post-incendio con diferentes edades. Los objetivos del estudio fueron: (1) establecer y determinar el 
rango histórico de la ocurrencia de incendios forestales, (2) evaluar parámetros dendrométricos para 
conocer la estructura presente y así establecer sitios de referencia para programas de restauración 
ecológica y (3) incrementar nuestro conocimiento en los procesos de disturbio en los bosques mixtos 
de pino-encino en la SMO. El estudio revela que existieron diferencias significativas entre las cohortes 
post incendios en términos de riqueza de especies, composición y estructura del rodal. La estructura 
forestal fue menos diversa en rodales jóvenes comparados con los intermedios, maduros y sobrema- 
duros. Una alta diversidad de especies arbóreas fue generalmente encontrada en rodales jóvenes, 
mientras que en los rodales intermedios y maduros mostraron la diversidad más baja. 

PALABRAS CLAVE: 

Análisis dendrocronológico, cronología de incendios, estructura forestal, incendios forestales, restau- 
ración, Sierra Madre Oriental. 


ABSTRACT 

Forest fires have had a remarkable relevance since the extraordinary fire season in 1998, due 
to their direct contribution to deforestation, changes in forest structure, species composition, and 
recently to their impact on the increment of carbón dioxide in the atmosphere. Nevertheless, there is 
a lack of detailed information about fire impact on forest structure, species composition, and forest 
dynamics in the Sierra Madre Oriental (SMO). In order to explain possible changes in those factors, 
this study has examined a chronosequence (134 years), which was determined by analysing tree ring 
samples in post-fire cohorts of different ages. The objectives of this study were to: (1 ) determine forest 
fire historical range of occurrence, (2) quantify present forest structure and tree species composition 
along a fire chronosequence in order to determínate reference areas to assist restoration programs, 
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and (3) increase the understanding of 
disturbance patterns in mixed pine-oak 
forest in this región. The study revealed 
significant differences between post-fire 
cohorts in terms of woody species rich- 
ness, stand composition, and structure. 

High tree diversity was generally found in 
young stands, while the intermedíate and 
major ones showed the lowest diversity. 

In oíd stands, tree diversity was incre- 
ased again. 

KEY WORDS: 

Dendrochronological analysis, master 
fire chronology, forest structure, forest 
fire, restoration, Sierra Madre Oriental. 

INTRODUCCIÓN 

El fuego al igual que otros disturbios (e. g. 
viento, agua) es un factor que se encuentra 
presente de manera natural en muchos 
ecosistemas y paisajes a lo largo del mundo 
(Whelan, 1995). Sin embargo, su aparición 
y empleo dentro de los ecosistemas fores- 
tales es vista de manera controversial en 
áreas de manejo y conservación biológica 
(Pyne et al., 1 996). Esto ha llevado a que se 
realicen investigaciones científicas en 
ecología de fuegos en algunos países, 
especialmente en los Estados Unidos, 
Canadá y Australia. Estos países han 
desarrollado estrategias para manejar y 
suprimir incendios forestales, pero también 
han aprendido a aceptar e integrar la 
ocurrencia de estos eventos naturales en 
programas de restauración, conservación y 
manejo (Wouters, 1993; Bergeron et al., 
2003; Stephens y Ruth, 2005). 

Para el caso de México los incendios 
forestales son un fenómeno ampliamente 
distribuido en los ecosistemas, principal- 
mente en áreas naturales protegidas, 
tales como parques nacionales, reservas 
de la biosfera entre otros (Jardel et al., 
2003; Rodríguez y Fulé, 2003). De 
acuerdo con las estadísticas oficiales en 
el año de 1998 cerca de 9000 incendios 
devastaron aproximadamente 220 000 ha 
de bosques. Para el año 2003 se consu- 
mieron cerca de 100 000 ha de áreas 


boscosas y hacia el mes de octubre de 
2006 se han registrado 8 657 incendios 
forestales afectando una superficie de 
alrededor de 240 000 ha (conafor, 
2006). 

En la Sierra Madre Oriental (smo), 
especialmente en los ecosistemas de 
bosques mixtos de pino-encino los incen- 
dios forestales tienen diversos efectos 
sobre la vegetación como consecuencia 
de la gran complejidad que presentan las 
comunidades vegetales y de la interfe- 
rencia ejercida por las acciones antropo- 
génicas como son las actividades agrí- 
colas, silvícolas y de pastoreo (conafor, 
2006). Aunque también debido a las dife- 
rentes respuestas de la vegetación al tipo 
de intensidad de cada fuego, la estación 
del año en la que tuvo lugar y la 
frecuencia del evento (Agee, 1993; 
Whelan, 1995). 

Al conjunto de estos elementos 
descritos anteriormente y sus combina- 
ciones se conoce como el régimen de 
incendios de un determinado lugar (Agee, 
1993; Whelan, 1995). En este sentido el 
uso del conocimiento histórico provee a 
largo plazo perspectivas para entender 
patrones y procesos ecológicos dentro de 
las comunidades vegetales. El empleo de 
este tipo de información se ha incremen- 
tado entre los científicos y manejadores 
de recursos para determinar el rango 
histórico y variabilidad de los procesos 
ecológicos para su evaluación y la toma 
de decisiones dentro de programas de 
restauración ecológica (Kaufmann et al., 
1994; Landres et al., 1999; González, 
2005b). 

Actualmente, diferentes técnicas son 
empleadas para identificar, datar y 
describir la incidencia histórica del fuego 
en una región determinada (Agee, 1993; 
González, 2005b) y para una gran 
cantidad de ecosistemas en diferentes 
regiones del mundo (Swetnam et al., 
2001; Johnson, 1992). Estás técnicas son 
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muy variadas e incluyen diversos 
elementos del paisaje como son: (1) 
técnicas dendrocronológicas que incluyen 
el muestreo de cicatrices en los fustes de 
los árboles quemados y de manera parti- 
cular en sus anillos de crecimiento (e.g. 
Dieterich y Swetnam, 1984; Fulé y 
Covington 1994; González, 2005a; 
Stephens et al., 2003); (2) muestreo de 
sedimentos de polen y fragmentos de 
carbón en lagos y pantanos (e.g. Cope y 
Chaloner, 1985); y (3) empleando docu- 
mentación oral y escrita (e.g. Lorimer, 
1980). Sin embargo, estas técnicas aún 
no han sido empleadas en la smo para 
determinar un régimen histórico del fuego 
y su impacto en ecosistemas forestales 
(González, 2005a; Rodríguez y Fulé, 
2003). 

El presente estudio tiene como obje- 
tivo el determinar el historial de incendios 
en una fracción de la smo para: (1) esta- 
blecer y determinar el rango histórico de 
la ocurrencia de incendios forestales, (2) 


evaluar parámetros dendrométricos para 
conocer la estructura presente y así esta- 
blecer sitios de referencia para 
programas de restauración ecológica y 
(3) incrementar nuestro conocimiento en 
los procesos de disturbio en los bosques 
mixtos de pino-encino en la smo. 

METODOLOGÍA 
Área de estudio 

En el presente trabajo se hace referencia 
al Parque Ecológico Chipinque (pech), el 
cual es una fracción del Parque Nacional 
Cumbres de Monterrey localizado dentro 
en la smo (Fig. 1). El pech tiene una 
superficie de 1624 ha (25°34' - 25°38' N, 
100°18' - 100°24' W) y presenta una 
topografía abrupta con una altitud que 
varía entre 650 y 1800 msnm. Los tipos 
de vegetación presentes son matorral 
submontano, bosque de encino y bosque 
de pino-encino (Rzedowski, 1978). 



Figura 1. Localización del área de estudio (pech). 
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Colecta de datos y análisis 

La revisión intensiva de literatura especia- 
lizada no aporta mayores antecedentes 
respecto a la ocurrencia y efectos del 
fuego al interior del pech. De esta manera 
se realizaron recorridos en campo para 
detectar evidencia visible acumulada por 
el fuego de acuerdo con las técnicas 
descritas por Agee (1993, Fig. 2) en 
árboles vivos y tocones para las especies 
representativas en el pech ( Pinus teocote, 
Pin us pseudostrobus y Quercus spp.). 
De esta manera, la combinación de las 
cicatrices en los árboles sobrevivientes, 
junto con el recuento de los anillos de 
crecimiento, los cuales cubren parcial- 
mente la herida permite datar incendios 
con una resolución anual (Agee, 1993; 
Dieterich y Swetnam, 1984). 

Producto de ello, se lograron esta- 
blecer cinco rodales (cohortes) post- 
incendio (Tabla 1) basados en la estruc- 
tura de edades del renuevo forestal 
(Johnson, 1992). 


El muestreo se llevó a cabo dentro 
de las cohortes identificadas en el pech, 
empleando una red de puntos estable- 
cidos sistemáticamente para determinar 
simultáneamente el historial de incendios 
y la estructura presente del bosque. El 
espaciamiento de la red de puntos fue de 
300m x 300m en cohortes de gran exten- 
sión y de lOOm x lOOm en aquellas de 
menor extensión (Fulé y Covington, 
1996). 

Un total de 23 puntos de muestreo 
se establecieron en las cinco cohortes 
post incendio identificadas (Tabla 1). 
Posteriormente, a lo largo de los tran- 
sectos determinados por la red de puntos 
y dentro de un radio de 15 m se marcaron 
todos los árboles que presentaron alguna 
evidencia de incendios para proceder con 
la extracción de virutas con el taladro de 
Pressler (Fig. 3). 

Para establecer un historial de 
incendios lo más preciso posible, se obtu- 
vieron en cada cohorte post incendio 


corteza 


restos de cicatriz del 
carbón primer incendio 


cicatriz del cicatriz del 

segundo año tercer año 



Figura 2. Ejemplo de la evidencia dejada por incendios forestales en el fuste de 
coniferas (Modificado de Morrison y Swanson, 1990) 
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Tabla 1. Características de los cohortes distribuidos en el pech. 


Cohortes 

Unidades 

de 

muestreo 

Año 

incendio 

tsf* 

Pendiente 

n 

Elevación 

(m) 

Virutas 

extraídas 

PECH98 

5 

1998 

4 

23-36 

1220 

0 

PECH84 

5 

1984 

18 

30-36 

1300 

11 

PECH72 

6 

1972 

30 

22-23 

1380 

15 

PECH40 

4 

1940 

62 

27 

1300 

8 

PECH868 

3 

1868 

134 

28 

1195 

30 


*tsf= tiempo sin la ocurrencia de incendio 


virutas de árboles en pie con al menos 
una marca visible en el fuste provocada 
por la presencia de fuego (Amo y Sneck, 
1977). Cada viruta fue procesada de 
acuerdo al método estándar de pulido 
para así incrementar la visibilidad de los 
anillos bajo el microscopio (Agee, 1993; 
Amo y Sneck, 1977). Las cicatrices ( fire 
scars) fueron identificadas de acuerdo a 
las anomalías características en los 
anillos de crecimiento (Dieterich y 
Swetnam, 1984). El año preciso de cada 
incendio fue determinado utilizando el 


conteo de anillos con las virutas colec- 
tadas (Strokes y Smiley, 1997). 

Para determinar la composición y 
estructura actual de los cohortes post- 
incendio, se establecieron parcelas tempo- 
rales de 0,1 ha donde todas las especies 
leñosas con un diámetro > 5 cm fueron 
registradas. El diámetro fue medido a la 
altura de pecho (dbh) y el coeficiente de 
variación de la variable diámetro fue 
empleada para medir la complejidad 
estructural de las cohortes junto con la 


corteza 



cicatriz 



t 


Figura 3. Formas de extracción de las virutas de árboles con anomalías ( fire scars) en 
los anillos de crecimiento (Modificado de Agee, 1993) 
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riqueza de especies (S), la cual está repre- 
sentada por el total de especies leñosas en 
cada parcela temporal de muestreo (Caroi 
et al., 2004). El índice de importancia (IVI) 
(i.e. el papel estructural) de cada especie 
fue calculado usando la densidad relativa, 
abundancia relativa y la frecuencia relativa 
(Magurran, 2004). Estadísticas básicas 
fueron usadas para comparar los prome- 
dios y estimar la variabilidad (Fowler et al., 
1998). Todos los valores concernientes a 
la estructura del rodal están expresados en 
± S.D. Para determinar si existen diferen- 
cias significativas entre la riqueza de espe- 
cies fue empleada la prueba no paramé- 
trica de Kruskall-Wallis acompañada de la 
prueba de comparación de medias de 
Dunn (Fowler et al., 1998). 

RESULTADOS 

Cohortes post incendio 

Mediante el análisis de los anillos de 
crecimiento obtenidos de las virutas 
colectadas y la estructura de edades 
dentro de las unidades de muestreo, se 
determinó el año de ocurrencia de 
disturbio (fuego) y al mismo tiempo corro- 
boró la fecha de origen de la cohorte. De 
esta manera se obtuvieron cinco cohortes 
post incendio las cuales en conjunto 
forman una cronosecuencia de 134 años 


donde la existencia de incendios se ve 
reflejada en la estructura presente del 
bosque (Tabla 1). Aunque la intensidad 
entre los incendios forestales fue muy 
variable y no fue determinada en este 
trabajo, el fuego fue responsable de 
destruir la vegetación en la mayoría de 
las cohortes dentro del pech. Dichos 
incendios fueron registrados y datados 
dentro del periodo que comprende entre 
los años de 1860-2002 en el área del 
pech (Fig. 4). 

Dentro de la cronología se observó 
que el primer evento se registró en el año 
de 1868, seguido por un amplio intervalo 
libre de incendios de 72 años. A partir del 
año 1940 se observan cuatro eventos 
donde el intervalo libre de incendios se 
reduce primeramente de 32 a 14 y 12 
años, respectivamente (Fig. 4). 

Riqueza de especies forestales 

A lo largo de los sitios, la diversidad de 
especies leñosas varió entre las cohortes 
de acuerdo al tiempo sin la ocurrencia de 
incendio. Para caracterizar los sitios a una 
escala más particular en cada cohorte se 
dividieron de acuerdo a su exposición, 
norte y sur. Los sitios con exposición norte 
mostraron la diversidad más alta de espe- 
cies leñosas que aquellas localizadas en la 
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Figura 4. Cronología de incendios para el pech la cual comprende el periodo de 

1868-2002 
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exposición sur (Fig. 5). La mayor riqueza 
de especies en la exposición norte fue la 
que se presentó en la cohorte más vieja 
(10,7±0,6 spp/1000 m 2 ), mientras que la 
cohorte de 40 años mostró la riqueza de 
especies más baja (3,7±1,2 spp/1000 m 2 ). 
La prueba no-paramétrica de Kruskal- 
Wallis mostró diferencias significativas 
(P<0.05) entre la media de la riqueza de 
especies de los sitios con exposición norte. 
Los sitios de la cohorte pech868 fueron 
significativamente diferentes de los sitios 
de la cohorte pech40 (PO.05) usando la 
prueba de Dunn. Entre los sitios con expo- 
sición sur, no se detectaron diferencias 
significativas al utilizar la prueba de 
Kruskal-Wallis. 


Estructura forestal 

La densidad media de plantas leñosas 
(dap ± 5 cm) en los rodales de exposición 
norte osciló desde 3400 ± 701 /ha (pech98) 
hasta 1066±120 (PECH868) individuos por 
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hectárea. Dentro de los rodales con expo- 
sición sur se registró una densidad menor 
que oscila entre 1265 ± 248/ha (pech98) 
hasta 850/ha (pech4o) (Tabla 2). 

El diámetro promedio (dbh) de los 
sitios con exposición norte se incrementa 
a medida que se incrementa el tiempo sin 
la ocurrencia de incendios, oscilando 
desde 3,1± 1,3 cm (pech98) hasta 22,7 ± 
10,6 cm (pech40). Sin embargo el 
diámetro promedio disminuye en los sitios 
maduros (PECH868). No obstante, los 
rodales maduros son más complejos que 
los rodales jóvenes e intermedios, debido 
a que el diámetro dentro del pech868 
presentó la mayor variabilidad, expresado 
por el coeficiente de variación (Tabla 2). 
Contrariamente en la exposición sur, los 
rodales jóvenes fueron estructuralmente 
más complejos que los rodales interme- 
dios y maduros. 

Estructuralmente, la especie más 
importante a lo largo de las cohortes fue 


Tabla 2. Parámetros dendrométricos de los cohortes dentro del pech. 
Se representan los valores medios ± desviación estándar (sd). 


Parámetros 

Exposición 

PECH98 

PECH84 

PECH72 

PECH40 

PECH868 

Densidad (ha-1) 

N 

34007701 

16207780 

11757247 

7537169 

10667120 


S 

12657248 

11257417 

14177140 

850 

— 

Diámetro (cm) 

N 

3,171,3 

7,772,5 

8,9?3,7n 

22,7710,6 

14,179,8 


S 

3,872,4 

7,071,9 

9,073,2 

22,0710,1 

— 

Altura (m) 

N 

2,170,7 

6,571,6 

7,071,9 

13,175,5 

8,179,8 


S 

2,771,2 

5,371,3 

6,771,9 

8,973,6 

— 

Área basal 







(m2 ha-1) 

N 

3,470,6 

8,374,1 

8,570,7 

37,375,3 

28,372,9 


S 

2,270,3 

4,772,2 

10,274,2 

38,970,5 

— 

CV1 

N 

0,46 

0,31 

0,41 

0,46 

0,69 


S 

0,63 

0,24 

0,36 

0,45 

— 


1 CV= Coeficiente de variación (diámetro) 
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Q. rysophylla, la cual presentó el índice 
de importancia más alto. P. teocote fue 
dominante en los rodales maduros y 
codominante en los rodales jóvenes e 
intermedios. P. pseudostrobus apareció 
como dominante en los rodales sobrema- 
duros (Tabla 3). La especie más impor- 
tante en las primeras etapas fue Q. 
rysophylla y Q. canbyi ya que los 
rebrotes (dbh< 5cm) eran más frecuentes 
en los rodales jóvenes. 

DISCUSIÓN 

El presente estudio reveló la estructura 
forestal encontrada en las cohortes post 
incendio en una fracción smo. Existieron 
importantes diferencias entre los cohortes 
post incendio en términos de riqueza de 
especies, índice de importancia, coefi- 
ciente de variación y estructura. 


Régimen de incendios en el pech 

De acuerdo a los registros de las virutas 
provenientes de los árboles incendiados, 
la aparición del fuego ha sido parte inte- 
gral de los bosques de pino-encino en la 
fracción de la smo (pech) durante los 
últimos 140 años. Con esto, los datos 
resultantes de la aplicación del método de 
la historia del fuego indican algunas rela- 
ciones coherentes entre los patrones 
espaciales y temporales de los incendios 
y la variabilidad climática. De esta 
manera, los registros en el marco crono- 
lógico (Fig. 4) mostraron la ausencia de 
incendios para el período comprendido 
entre 1860 y 1940. Sin embargo, este 
periodo libre de incendios coincidió con la 
información histórica documentada por 
Stephens et al. (2003), en la cual señalan 
que se presentó una disminución en la 
ocurrencia de incendios en el suroeste de 


Tabla 3. índice de Valor de Importancia (IVI) para las especies arbóreas. 


Especies 

PECH98 

PECH84 

PECH72 

PECH40 

PECH868 


NE 

SE 

N 

S 

N 

S 

NE 

s 

N 

Arbutus xalapensis 

5,2 

8,7 

3,8 

3,1 

4,0 

2,9 

— 

— 

3,3 

Ceanothus coeruleus 

6,1 

4,4 

— 

— 

— 

— 

— 

— 

— 

Cercis canadensis 

1,7 

— 

1,7 

— 

— 

— 

— 

5,3 

4,0 

Pinus pseudostrobus 

— 

— 

6,9 

4,6 

8,7 

3,3 

10,7 

0,0 

30,4 

Pinus teocote 

— 

— 

22,8 

11,4 

14,6 

23,4 

53,4 

42,7 

6,9 

Quercus canbyi 

15,4 

21,3 

6,0 

40,8 

6,3 

18,4 

14,2 

23,9 

18,7 

Quercus rysophylla 

47,2 

35,8 

39,9 

18,5 

39,9 

33,2 

18,7 

8,7 

19,9 

Quercus virginiana 

10,8 

26,2 

11,7 

2,9 

12,6 

— 

3,1 

6,0 

3,1 

Quercus laeta 

1,4 

3,7 

— 

14,7 

— 

12,2 

— 

8,4 

— 

Otras especies 

12,2 

— 

7,1 

4,0 

13,7 

6,5 

— 

5,1 

13,7 

Total 

100,0 

100,0 

100,0 

100,0 

100,0 

100,0 

100,0 

100,0 

100,0 


Otras especies= spp. 1, spp. 2, Juglans mollis, Juniperus flaccida, Prosopis glandulosa, Quercus laceyi, Q. poly- 
morpha.- No se registro 
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los Estados Unidos y el norte de México 
a inicios del presente siglo (1899/1900). 
Por otra parte Kitzberger et al. (2001), 
demostraron que el decremento en la 
ocurrencia de incendios forestales en 
esta época estuvo relacionado con 
cambios en la frecuencia y amplitud de El 
Niño-Oscilación del Sur (enso). 

A partir del año 1940 los intervalos 
libres de fuego varían desde 32 hasta 12 
años en las áreas del pech. Tal reducción 
en los intervalos podría ser resultado de 
la relación causal entre el aumento en la 
frecuencia de incendios y la expansión 
de la infraestructura para la explotación 
de madera y la iniciación de campañas 
para detección, prevención y supresión 
de incendios forestales en México (Fulé y 
Covington, 1994). 

Para los disturbios de 1984 y 1998, 
no se cuenta con información contun- 


dente si fueron iniciados de manera 
natural o antropogénica. Sin embargo, 
estos incendios se encuentran asociados 
con los años de 1983 y 1998 fecha en las 
cuales se presentaron índices positivos 
del fenómeno de El Niño-Oscilación del 
Sur donde intensas sequías y tempera- 
turas extremas afectaron el norte de 
México (Heyerdahl y Alvarado, 2003). 


Riqueza de especies leñosas 

Mientras que la riqueza de especies 
declina a medida que se incrementa el 
paso del tiempo en la cronología hasta 
llegar a escasa en los rodales maduros, 
en los rodales sobremaduros (pech868) la 
riqueza de especies aumenta nueva- 
mente (Fig. 5). El decremento en la 
riqueza de especies arbóreas pudo ser 
debido en primera instancia por la elimi- 
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Cohortes post-incendio 

Figura 5. Riqueza de especies leñosas a lo largo de las diferentes cohortes post 
incendio. Las barras de error representan sd. 
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nación de algunas especies pioneras, las 
cuales fueron suprimidas por especies de 
rápido crecimiento y de aquellas especies 
que poseen la habilidad de rebrotar 
(Miller, 2000). De tal manera que a 
medida que se incrementó la edad del 
rodal, los niveles de luz fueron menores al 
interior del sotobosque debido probable- 
mente por el cierre gradual del dosel 
como se ha presentado en otros estudios 
(e. g. Caroi et al., 2004). Tal efecto 
también fue documentado por Gilí et al. 
(1999) como una característica de 
bosques presentes en zonas áridas y 
semiáridas. 

En general, los sitios con exposición 
norte presentaron una mayor densidad de 
árboles y arbustos, lo cual promueve una 
complejidad estructural más diversa, 
comparada con los sitios de exposición 
sur. Los resultados de la investigación 
también fueron consistentes con las 
observaciones en un estudio similar por 
Guo (2001) en el chaparral californiano, 
donde los sitios con exposición norte 
mostraron una diversidad más alta que en 
los de exposición sur. 

La riqueza de especies es una 
medida de la diversidad, la cual es muy 
importante para la funcionalidad, estabi- 
lidad e integridad del ecosistema. Aparte 
de las razones ecológicas, mantener la 
diversidad en ecosistemas forestales 
posee un especial interés económico, 
científico y educacional (Magurran, 2004). 

Estructura forestal 

Se encontró que a lo largo de la cronose- 
cuencia la variable tsf tiene diferentes 
impactos en la estructura y composición de 
bosque, no obstante de las diferencias 
creadas por las condiciones geográficas 
como elevación, exposición y pendiente. El 
área basal total para cada cohorte se 
incrementó a medida que aumenta la 


variable tsf (Tabla 2), como es de espe- 
rarse la dinámica de crecimiento de 
bosques secundarios (Barnes etal., 1988). 

La mayoría de los cambios en la 
vegetación arbórea a lo largo de la crono- 
secuencia fue causada por el incremento 
en el índice de valor de importancia de la 
especie P. pseudostrobus, P. teocote y Q. 
canbyi (Tabla 3) y el decremento significa- 
tivo en las especies de Q. rysophylla y Q. 
virginiana. La tendencia para estas dos 
especies de encino fue un decremento a 
partir de la cohorte PECH72 aproximada- 
mente 30 años después del disturbio (Tabla 
3). Esta tendencia probablemente repre- 
senta el período del cierre de copas, la 
disminución de la densidad e incrementos 
en área basal relativamente lentos, 
factores que afectan otras propiedades de 
los ecosistemas (Barnes et al., 1 988). 

CONCLUSIONES 

Alrededor del mundo existen numerosas 
evidencias de que los incendios forestales 
desempeñan un papel importante y posi- 
tivo en ciertos ecosistemas (Agee, 1993; 
Pyne et al., 1996). Sin embargo, en 
México se relacionan con la pérdida de 
hábitat, diversidad y con la devastación 
de la belleza escénica. No obstante, el 
presente trabajo demuestra mediante la 
reconstrucción histórica de los incendios 
forestales combinada con la evaluación 
ecológica y silvícola de las áreas en consi- 
deración, que los incendios forestales no 
necesariamente afectan la biodiversidad 
sino que contribuyen a la misma y además 
promueven la heterogeneidad del ecosis- 
tema. De ahí, que la reconstrucción histó- 
rica de los incendios forestales adquiera 
una gran importancia, puesto que ayuda a 
la generación de conocimientos especí- 
ficos sobre el disturbio y sus efectos 
ecológicos en las condiciones físico- 
geográficas, ecológicas y sociales de una 
región en particular. 
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Cambios en la estructura y biodiver- 
sidad fueron registrados en las diferentes 
etapas sucesionales del bosque. Por lo 
anterior, en el caso particular del presente 
estudio, los resultados refuerzan la idea 
de que los incendios forestales han sido 
parte importante e integral para mantener 
la diversidad estructural y de especies en 
los bosques de pino-encino de la región. 

El monitoreo de las áreas incendiadas 
facilitaron suficiente información para inte- 
grar y establecer áreas o sitios referencia 
para futuros programas de restauración de 
ecosistemas. Sin embargo, la nueva polí- 
tica establecida para la supresión sistemá- 
tica de incendios en México podría modi- 
ficar e incluso aumentar el riesgo de que se 
presenten eventos más severos tanto en 
intensidad como en tiempo. De tal manera 
que se podría contar con un menor número 
de sitios de referencia y áreas para deter- 
minar la evidencia histórica de la ocurrencia 
natural de los incendios forestales. 

Por último, la acción del fuego sobre 
el ambiente ha experimentado cambios 
drásticos en su concepción misma. Su 
papel ecológico en el mantenimiento y 
funcionamiento de ecosistemas ha sido de 
cierta manera ensombrecido por los 
grandes desastres. Pero la investigación 
ha revelado que para los ecosistemas 
donde naturalmente ocurren incendios 
forestales, el fuego desempeña un papel 
importante para mantener una diversidad 
estructural y en especies de los mismos 
ecosistemas. Un incendio forestal no es 
necesariamente un desastre, solamente se 
podría considerar como tal cuando se 
enfrenta con los intereses económicos y 
sociales del bosque. Es por eso que las 
organizaciones y agencias encargadas de 
los recursos naturales deben crear 
acciones encaminadas a promover polí- 
ticas de educación y difusión a distintos 
niveles socioculturales, capacitar a la 
población rural sobre el uso responsable 
del fuego y sobre todo fomentar la investi- 
gación en materia de incendios forestales. 
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El tamaño sí importa: 

los frutos grandes de Juniperus deppeana 
Steud. (sabino) son más susceptibles a 
depredación por insectos 

Size matters: the Juniperus deppeana Steud. (sabino) 
big fruits are more sensitive to insect depredation 

Armando J. Martínez 1 >* Pilar Sainos 1 , Ernesto Lezama Delgado 1 

y Guillermo Angeles-Álvarez 2 

RESUMEN 

El siguiente trabajo documenta la relación entre el tamaño de los frutos (gálbulos) de Juniperus 
deppeana Steud. y el daño ocasionado por insectos depredadores de semillas en un fragmento de 
bosque de sabinos, característico del estado de Tlaxcala, México. En el estudio se seleccionaron 
árboles femeninos del sabino para cuantificar la maduración de los gálbulos y se clasificó la disponi- 
bilidad de éstos en tres categorías de tamaño: pequeños (5 a 9 mm), medianos (10 a 12 mm) y 
grandes (13 a 15 mm), cuantificando al mismo tiempo el número y la densidad de gálbulos intactos y 
aquellos que habían sido dañados por insectos, comparaciones que se realizaron durante cinco mués- 
treos. La producción de gálbulos en J. deppeana fue continua durante las cinco fechas de colecta. 

Entre 21 a 46 de los individuos presentaron frutos en 50-70% del árbol (por estimación visual). 
También, las categorías de 5 a 30% de abundancia de gálbulos en las copas la presentaron entre 19 
a 31 de los individuos femeninos de sabinos durante las cinco fechas de cuantificaciones porcen- 
tuales. Los sabinos hembras produjeron frutos continuamente de agosto de 2004 a agosto de 2005. 
Aunque hubo siempre una amplia disponibilidad de gálbulos grandes, pequeños y medianos, fueron 
los de mayor tamaño los más dañados por insectos, ya que hasta 30% de los frutos grandes presen- 
taban orificios de salida de insectos depredadores de semillas. 

PALABRAS CLAVE: 

Biología reproductiva, gálbulos, insectos semillívoros, Juniperus deppeana , sabino. 


ABSTRACT 

The relationship between fruit size of Juniperus deppeana and damage caused by seed 
predator insects, in a forest fragment of Juniperus deppeana (sabino), typical of the State of Tlaxcala, 
México was analyzed. Female trees were selected to quantify fruit maturation. Their availability were 
classified into three categories: small (5 to 9 mm), médium (10-12 mm), and large (13 to 15 mm). At 
the same time, the amount and density of intact fruits and those presenting insect damage were quan- 
tified on five different samplings. Between 21-46 of the individuáis presented fruits in 50 to 70% of the 
tree (by visual estimation). Also, categories of 5 to 30% of galbule abundance in the crowns were 
presented in 19 to 31 of the female individuáis during the five different dates when the per cent quan- 
tifications were made. Female individuáis produced fruits continuously from august 2004 to august 
2005. Although there was always a wide availability of large, médium, and small fruits, the first ones 
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were the most damaged, since up to 
30% of the fruits ¡n that category 
presented exit holes caused by seed- 
predator insects. 

KEY WORDS: 

Reproductive biology, fruit productivity, 
seed-predator insects, Juniperus 
deppeana , sabino. 

INTRODUCCIÓN 

Los frutos son una fuente importante de 
alimento en diversos ambientes, para 
distintos grupos de animales. Entre los 
consumidores de frutos más evidentes se 
encuentran los que se alimentan del 
mesocarpo y son los que en su mayoría 
dispersan las semillas (Willson, 1992), 
pero en esta interacción también podrían 
estar actuando insectos depredadores 
(Janzen 1969, 1971a). Es así como 
muchas de las interacciones antago- 
nistas, como la de los insectos hacia los 
frutos, tienen un fuerte impacto sobre el 
desarrollo de estos últimos (Herrera, 
1982; Janzen, 1977) y sobre la dispersión 
de semillas. Procesos que han sido 
objeto de amplio interés ecológico porque 
pueden reducir la progenie de distintas 
especies vegetales (Sallabanks y 
Courtney, 1992). En este sentido, algunas 
poblaciones de insectos pueden llegar a 
tener un gran impacto como depreda- 
dores de la pulpa y semillas de los frutos, 
lo cual repercute en el desarrollo de inter- 
acciones mutualistas fruto-dispersor, 
debido a que un fruto, cuando está 
dañado, es menos atractivo para los 
vertebrados frugívoros (Mazur y 
Courtney, 1984; Sallabanks y Courtney, 
1992; Valburg, 1992ab). 

En este contexto, el nivel de la inten- 
sidad del daño en los frutos y el consumo 
de semillas predispersión están definidos 
como una interacción antagonista que 
tiene amplia importancia en los paráme- 
tros poblacionales, ya que reduce el 
reclutamiento de nuevos individuos 


(Jordano et al., 2004). Estos dos sistemas 
de interacción tienen un fuerte vínculo 
dentro del marco reproductivo de la 
planta. Es así como el efecto del balance 
entre el consumidor de frutos y semillas 
vs. la labor de los dispersores, es el resul- 
tado de una evidente dualidad que 
presentan los frutos como fuente de 
alimento para distintos taxa con papeles 
funcionales distintos. Además, evolutiva- 
mente, el fruto de las plantas es el punto 
de reunión de las presiones selectivas por 
parte de parásitos y dispersores, por lo 
que sus rasgos y características como 
estructura alimenticia representan costos 
y beneficios directos al sistema reproduc- 
tivo de la planta (Herrera 1984,1986). 

El sabino es un claro ejemplo de 
planta del grupo de las gimnospermas que 
produce gálbulos, los cuales funcional- 
mente se despliegan sobre las copas de 
las plantas como frutos carnosos, homolo- 
gando a las angiospermas. Los frutos de 
algunas especies de Juniperus europeas y 
del sur de los Estados Unidos de América 
son consumidos por aves y mamíferos 
frugívoros (Livingston, 1972; Schupp et al., 
1997) aunque también son consumidos 
por insectos que dañan la pulpa o semillas 
antes de su dispersión (Traveset y Sans, 
1994; García, 1998; Guido y Roques, 
1 996). Por ejemplo, en J. communis la tasa 
de daño causado a sus semillas por depre- 
dadores como Megastigmus bipunctatus, 
no se relacionó con características del 
fruto; en cambio, el mayor daño a los frutos 
causado por Carulaspis juniperi, se rela- 
cionó con gálbulos pequeños, de peso 
ligero y baja proporción de pulpa:semillas 
(García, 1998). 

El sabino es una especie dioica que 
ocupa el segundo lugar en distribución en 
México, pues habita desde ambientes 
templados hasta semiáridos (Martínez, 
1963). El estado de Tlaxcala tiene una 
extensión aproximada de 16288 hectáreas 
de superficie forestal de bosques de J. 
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deppeana (semarnap-uach, 1999; cona- 
fora, 2003). A pesar de esto, existe muy 
poca información acerca de los aspectos 
biológicos y ecológicos de esta especie 
(Reyes, 2003). Por lo tanto, los gálbulos 
frescos de Juniperus representan un 
sistema de estudio idóneo para escudriñar 
algunos aspectos básicos sobre interac- 
ción ecológica planta-insecto y para inves- 
tigar el efecto del impacto de los insectos 
consumidores de semillas en relación al 
tamaño de los gálbulos. 

Existen algunas investigaciones en 
donde se menciona el efecto de los 
insectos sobre distintas estructuras del 
sabino, entre los que están el estudio de 
Vázquez-Yanes et al. (2002), quienes 
señalan que J. deppeana es resistente al 
ataque por termitas. En cambio Cibrián et 
al. (1995), registraron que Frankliniella 
molesta (Thripidae) daña los estróbilos 
tanto masculino como femenino. Salazar- 
Figueroa (2001), determinó reciente- 
mente que las larvas de un insecto de la 
familia Gelechiidae destruyen hasta un 
50% de los conos de los sabinos. 

OBJETIVO 

Examinar la producción de gálbulos de J. 
deppeana durante un período anual, 
además de evaluar la relación entre cate- 
gorías de tamaño y la cantidad de frutos 
intactos y dañados por insectos consumi- 
dores de semillas en una población de 
sabinos del estado de Tlaxcala. 

MATERIALES Y MÉTODOS 

Descripción de la especie 

El género Juniperus se ubica dentro de la 
familia Cupressaceae. Son árboles 
perennes con hojas escuamiformes 
opuestas o verticiladas. Están reportadas 
aproximadamente 50 especies para 
América, Asia, Europa y África (Farjon et 


al., 1993). J. deppeana es una especie 
dioica y su fruto es una gálbula de color 
rojizo-canela hasta rojizo-moreno con un 
diámetro de 8 a 15 mm, el cual no se abre 
aunque llegue a la madurez. Es posible 
que exista variación en el tamaño de los 
frutos entre poblaciones. Estas estruc- 
turas están colocadas sobre pedúnculos 
cortos, escamosos, macrocárpicos, polis- 
permos (Martínez, 1948). Tienen sabor 
agridulce y son de color violáceo o rojizo, 
con una capa de cera. El mesocarpio es 
fibroso y fuertemente adherido a las semi- 
llas, las cuales varían de una a cuatro por 
gálbula. Las hembras de los insectos 
consumidores de semillas ponen general- 
mente un huevo por semilla, después 
pupan las larvas y los adultos emergen 
abandonando los frutos a través de un 
agujero de salida típicamente circular, al 
igual que en otras especies de Juniperus 
(García, 1998; Vikberg, 1966). Las semi- 
llas son dispersadas por algunos mamí- 
feros (Lezama, 2007) y supuestamente 
pasan por el tracto digestivo de esos 
mamíferos sin sufrir daño aparente 
(Sainos, 2007). El sabino florece de 
febrero a marzo (Vázquez-Yanes et al., 
2002) sus gálbulos maduran entre agosto 
y octubre del segundo año de desarrollo 
(Salazar-Figueroa, 2001) y una vez 
maduros, los gálbulos ya no aumentan de 
tamaño. En México hay registradas 14 
especies de Juniperus (Oldfield et al., 
1998; Rzedowski, 1986), distribuidas en 
una amplia extensión territorial (Farjon et 
al., 1993; Zanoni y Adams, 1979; 
Rzedowski, 1986) y donde J. deppeana 
destaca, puesto que crece en Durango, 
Chihuahua, Jalisco, Zacatecas, Veracruz, 
Hidalgo, Puebla y Tlaxcala, formando el 
bosque escuamifolio. Varía desde 
pequeños matorrales hasta rodales de 15 
m de altura, establecidos en suelos 
pobres y calizos. Los bosque de sabinos 
en Tlaxcala son sistemas vegetales fuer- 
temente perturbados y que han sido 
desplazados por la agricultura y plantíos 
de Eucalyptus globulus. Es común encon- 
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trar a J. deppeana asociada con especies 
vegetales del matorral xerófito, bosques 
de encino y especies no nativas (inegi, 
1998; Ortega, 2000). 

Área de estudio 

El área de estudio fue la zona boscosa 
anexa a la Facultad de Biología de la 
Universidad Autónoma de Tlaxcala, la 
cual se ubica a 19° 1976” de latitud norte 
y 98°1676” de longitud oeste, a una 
altitud de 2300 m. El clima es templado 
subhúmedo con un régimen de lluvias de 
mayo a septiembre, donde la precipita- 
ción promedio anual es de 500 a 700 mm 
y con una incidencia máxima de lluvia 
entre junio y agosto, la cual está acompa- 
ñada de fuertes ráfagas de viento. El 
periodo más seco corresponde a los 
meses de diciembre y enero, con un valor 
promedio de 7 mm de lluvia. El período 
más caluroso es de marzo a mayo y el 
más frío es de diciembre a enero, con una 
temperatura de 11 a 12°C. La tempera- 
tura promedio mínima anual registrada es 
de 6,1°C y la máxima de 24,6°C, en tanto 
que la temperatura anual promedio es de 
13°C (spp, 1987). Los bosques de 
sabinos en Tlaxcala son sistemas vege- 
tales fuertemente perturbados y han sido 
desplazados por la agricultura y plantíos 
de Eucalyptus globulus. Es común encon- 
trar a J. deppeana asociada a especies 
vegetales del matorral xerófilo, bosques 
de encino y especies no nativas (INEGI, 
1998; Ortega, 2000). 

Distribución y tamaño de los sabinos 

En una hectárea se ubicaron, mediante 
coordenadas cartesianas (x, y), todos los 
árboles de sabino a los cuales se les 
midió la altura, diámetro basal y cobertura 
de la copa. Esta última variable fue esti- 
mada utilizando la longitud del eje mayor 
y menor de la proyección de la copa, utili- 
zando la fórmula de la elipse ( C= 0.25 


p .DM* dm). Los datos fueron analizados 
utilizando sólo los árboles seleccionados 
al azar para el seguimiento del estudio de 
producción de frutos mediante regre- 
siones lineales (Zar, 1999), para 
comparar las relaciones de tamaño entre 
pares de variables de respuesta 
(diámetro basal vs. altura, por ejemplo). 

Categorías de gálbulos 

Después de marcar los árboles de J. 
deppeana, estos fueron categorizados 
funcionalmente como árboles machos, 
hembras y juveniles. Posteriormente, con 
sólo las plantas hembras y productoras 
de gálbulos en distinta etapa de 
desarrollo se realizaron estimaciones 
visuales de la producción de frutos, utili- 
zando la siguiente escala: categoría 1 (de 
5 a 30%), 2 (de 30 a 50%), 3 (de 50 a 70 
%) y 4 (de 70 a 100%), en septiembre de 
2004, y marzo y septiembre de 2005. El 
criterio de la categorización fue a partir 
del volumen de la copa del árbol que 
contenía gálbulos. Posteriormente con 
estos datos se calcularon los porcentajes 
de árboles en dichas categorías y se 
compararon entre las tres fechas de 
registro para el mismo número de árboles 
en cada fecha (n=117). 

Producción de gálbulos en ramas 

Se realizó una selección aleatoria de 30 
plantas de los 117 árboles hembras 
etiquetados, en los cuales se marcaron al 
azar dos ramas de una longitud de 30 cm 
de largo desde el área apical hacia la 
base del dosel inferior de la copa. Esto 
permitió homogeneizar el tipo de rama 
para tener menor variabilidad entre 
árboles y así realizar el seguimiento de la 
maduración de los gálbulos. Para esto, en 
cada rama se contó el número de frutos 
producidos, separándolos por su tamaño 
en tres categorías: pequeños (de 5 a 9 
mm), medianos (10 a 12 mm) y grandes 
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(12 a 15 mm). Se efectuaron cinco regis- 
tros de conteo durante un año. Esta 
variable de respuesta fue analizada utili- 
zando análisis de varianza de medidas 
repetidas (andevaiti) con Modelos Line- 
ales Generalizados (glm, Statistica Stat- 
Soft, Inc., 6.1 1984-2003) empleando el 
siguiente modelo de factores fijos y 
anidados: 

número de frutos = tamaño + árbol j rama j + Fecha 
* árbo1 [rama] + error . 

donde: tamaño = categoría de gálbulo y 
fecha la medida repetida con las respec- 
tivas interacciones, y entre corchetes se 
ubica el factor anidado aleatoriamente. 
Además, la variable de respuesta conteo 
de gálbulos fue transformada a rangos 
para no violar el supuesto de normalidad 
y homogeneidad de varianzas (Conover e 
Imán, 1981). 

Densidad de gálbulos 

En la hectárea donde se realizó el censo, 
otras 30 plantas fueron seleccionadas al 
azar. En un punto al azar bajo el dosel de 
la copa de cada árbol hembra se colocó 
una red de malla de forma circular, cuya 
área abarcaba una superficie de 0,25 m 2 , 
para colectar los frutos desprendidos de 
la planta madre y cuantificar la disponibi- 
lidad de gálbulos intactos y los que 
presentaban daños por insectos. Estas 
dos variables de respuesta que fueron 
gálbulos intactos y dañados, fueron anali- 
zadas utilizando análisis de varianza no 
paramétricos de medidas repetidas de 
Friedman (Zar, 1999). 

Estimación de la producción de 
gálbulos 

Posteriormente, a partir de la estimación de 
la cobertura de la copa de los árboles 
hembras y con los datos de número de 
gálbulos colectados en las trampas de cada 
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árbol se obtuvo la densidad total por árbol, 
cálculo que se realizó a partir de la densidad 
de frutos colectados en la redes de 0.25 m 2 , 
extrapolándola para el área total de cober- 
tura de la copa de cada árbol. Así, se obtuvo 
una medición aproximada de la producción 
de frutos intactos y dañados, variables de 
respuesta que usamos de manera indepen- 
diente para un análisis de regresión lineal 
cuando los gálbulos no tienen daño, y poli- 
nomial de segundo orden para frutos 
dañados (Zar, 1 999). 

RESULTADOS 

La altura de los árboles de J. deppeana 
en los que se colocaron las redes para 
colectar los gálbulos fue en promedio de 
6 m ± 0,18 EE, los cuales tienen una 
cobertura promedio de 27 m 2 ± 2,26 EE y 
de diámetro de la base del tronco 
promedio de 30 cm ± 0,9 EE. Asimismo, a 
partir de las regresiones lineales se deter- 
minó que la altura de los árboles utili- 
zados para el seguimiento fenológico no 
se relacionó con el diámetro de la base 
del tronco (F 1 2 s = 4, r 2 - 0,08 P> 0,06), ni 
tampoco fue significativa la relación entre 
el diámetro basal vs. cobertura de la copa 
(Fi 28 = 2, r 2 - 0,02, P- 0,2), de igual 
forma tampoco hubo asociación entre la 
altura vs. cobertura (F 1 28 = 4,2, r 2 = 0,09, 
P> 0,05). Esto sugiere por un lado que no 
hay relación entre variables de tamaño de 
los árboles seleccionados, además de ser 
homogéneos al menos en la cobertura de 
la copa. 

Categorías de gálbulos 

De los 328 árboles de J. deppeana 
marcados dentro de la hectárea de 
estudio el 35,6% eran individuos feme- 
ninos y el 49,4% masculinos; el restante 
15% de las plantas registradas eran juve- 
niles. Estas proporciones se mantuvieron 
con cierta regularidad durante el año de 
monitoreo de la producción de frutos. 
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Durante septiembre de 2004 se 
registró que 45 de ios árboles hembras se 
situaron dentro de la categoría de 50 a 70% 
de abundancia de gálbulos sobre la copa 
de los árboles, en tanto, que 35 de los 
sabinos registraron una presencia de frutos 
dentro de la categoría de 30 a 50%. En 
cambio, la categoría de 5 a 30% sólo se 
presentó en 19 de los árboles y los 18 
árboles restantes de sabino hembras estu- 
vieron en la categoría de 70 a 100%. Estos 
conteos de árboles con abundancia de 
gálbulos se mantuvieron muy similares 
durante la medición de marzo de 2005. 
Aunque el número de árboles en las dife- 
rentes categorías varió en relación a la 
reducción o aumento de plantas en las 
escalas porcentuales con respecto al 2004. 
En cambio durante septiembre del mismo 
año, la categoría de 30 a 50% fue la que 
presentó el mayor número de árboles 
(n=39), seguida de la categoría de 5 a 30% 
en la cual se registraron (n=31 ). Empero, la 
categoría 50 a 70% presentó una sustan- 
cial reducción de árboles hembras con tan 


sólo 21 plantas y el restante de los sabinos 
(n=26) correspondió a la categoría de 70 a 
100% (Fig. 1). 

Producción de gálbulos en ramas 

La comparación del número de gálbulos 
producidos en ramas de J. deppeana en 
tres diferentes categorías de tamaño, 
indica que hay diferencias significativas 
entre los diferentes tamaños de gálbulos 
(F 2120 P< 0,0001), además de 

haber un efecto significativo de la interac- 
ción tamaño x fecha (F 823 4 =23,7, 
P< 0,0001). Asimismo, el húmero de 
gálbulos chicos fue poco abundante 
durante las cinco fechas de registro. Al 
contrario, los frutos medianos fueron los 
más abundantes en agosto de 2004 con 
respecto a la categoría de frutos grandes. 
No obstante, en diciembre del mismo año 
el promedio fue mayor para los gálbulos 
de la categoría grande, aunque en marzo 
y junio de 2005 nuevamente el número 
promedio de gálbulos fue mayor para la 
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Figura. 1. Número de árboles correspondiente a cuatro diferentes estimaciones de 
producción de gálbulos en valor porcentual, en tres distintas fechas, donde el número 

total de plantas hembras fue de 117. 
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categoría de tamaño mediano. Esto se 
debe a que en marzo, junio y agosto del 
2005 los frutos grandes se desprendieron 
de las ramas por el nivel de madurez. 
Aunque, durante la última fecha no hubo 
diferencias significativas entre estas dos 
categorías de tamaño, pero sí difirieron 
de los frutos pequeños (Fig. 2). Esto 
indica que el desarrollo posiblemente sea 
por pulsos entre categorías de tamaño 
hasta el momento en que se desprenden 
de la planta madre. 

Fenología de gálbulos 

El número de gálbulos obtenidos en las 
redes colocadas bajo la copa de los 
árboles hembras de J. deppeana varía a 
lo largo de las cinco fechas de registro. 
Así, los frutos sin daño (viables) de la 
categoría de tamaño pequeño no difieren 
significativamente entre las tres primeras 
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fechas de registro de agosto 2004 a 
marzo de 2005. Donde las dos últimas 
fechas tuvieron un mayor número de 
gálbulos al realizarse las comparaciones 
múltiples de Tukey (Fig. 3a). En contraste, 
la cantidad de frutos con daño para la 
misma categoría de tamaño mostró menor 
abundancia en las trampas durante todo el 
periodo de registro y fueron los dos 
últimos registros los que presentaron la 
mayor cantidad de gálbulos perforados 
(Fig. 3b). La categoría de tamaño de 
gálbulos medianos con y sin daño 
muestra el mismo patrón en la abundancia 
entre las fechas de registro (Fig. 3cd). 
Asimismo, el número de frutos dañados 
por insectos en la categoría de tamaño 
mediano fue significativamente menor en 
relación con la categoría mediana de junio 
y agosto (Fig. 3d). Este patrón se repite 
para la categoría de frutos grandes (Fig. 
3ef). No obstante, los valores de 
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Figura 2. Promedio de gálbulos para las tres categorías de tamaño de frutos, durante 
los cinco periodos de muestreo de agosto de 04 a agosto de 05. Las líneas verticales 

indican errores estándar. 
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Figura 3. Número de gálbulos con daño y sin daño durante las cinco fechas de registro 
dentro del fragmento de bosque de sabinos anexo a la Facultad de Biología 
Agropecuaria, ab) frutos de la categoría chicos, cd) medianos y ef) grandes. Las líneas 
horizontales, cajas y líneas verticales indican la mediana, cuartil e intervalos de valores 
y los círculos oscuros muestran los valores extremos. 


medianas para gálbulos dañados dentro 
de esta categoría fueron más altos en 
relación con la categoría de tamaño de 
frutos medianos y chicos (Fig. 3bd). 

Las densidades de gálbulos viables en 
las categorías de tamaño chico y 
mediano fueron las que presentaron los 
valores de medianas más altas durante 


junio y agosto de 2005. Además dentro 
de cada categoría de tamaño, hay dife- 
rencias significativas entre fechas (Fig. 
4) y el comportamiento de las medianas 
obtenidas en la categoría de frutos 
grandes también presenta diferencias 
significativas entre las medidas repe- 
tidas (Fig. 4e). No obstante, esta misma 
categoría de tamaño para la densidad 
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Figura 4. Densidad de gálbulos por (m2) con y sin daño durante cinco fechas de 
registro, ab) frutos de la categoría chicos, cd) medianos y cf) grandes. Las líneas 
horizontales, cajas y líneas verticales indican la mediana, cuartil e intervalos de valores 

y los círculos oscuros son valores extremos. 


de gálbulos dañados por insectos fue la 
que presentó mayor densidad al compa- 
rarla con las categorías de tamaño de 
gálbulos chicos y medianos (Fig. 4bd). 


Estimación de número de gálbulos 

La estimación de la producción de 
gálbulos de J. deppeana es mayor para el 
caso de frutos viables, en donde la rela- 
ción alométrica entre la cobertura de copa 
de los árboles está positivamente relacio- 


nada y el ajuste lineal muestra una alta 
densidad de gálbulos por metro cuadrado 
(y= -1 ,14e 4 +3285,8x, F 28 =15, 

P<0,0006, r 2 - 0,35). Asimismo, con el 
mismo valor de cobertura de las plantas, 
se relacionó la densidad de gálbulos 
perforados por insectos y en este caso el 
modelo fue un ajuste polinomial de 
segundo grado (y= -3,59e 4 + 4872,8-57,1 
x, F . | 28 =6, P= 0,007, ^=0,30). Por lo 
tanto, se obtuvieron diferencias en la esti- 
mación de la densidad de gálbulos en las 
dos condiciones de frutos (viables vs. 
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c 

Figura 5. Densidad estimada de gálbulos por (m 2 ) y su relación con la cobertura de 
copa de los árboles de sabino Juniperus deppeana. Los círculos oscuros representan 
gálbulos sin daño y los vacíos los gálbulos dañados. 


dañados, Fig. 5). Además, utilizando una 
prueba de t para datos no independientes 
se obtuvieron diferencias significativas 
(fi 29 = 3,4, P<0,001), donde los gálbulos 
perforados tienen en promedio 47213 ± 
5053,6 EE y los intactos 80493,8 ± 
12659,8 EE. Esto en general se repre- 
senta como distintos valores porcen- 
tuales, ya que el 63% correspondió a los 
gálbulos sin daño, mientras que los perfo- 
rados por insectos presentan un 37%. 

DISCUSIÓN 

Los sabinos presentaron producción de 
gálbulos todo el año, aunque de forma 
variada y con distintas cantidades en las 
categorías de tamaño de los frutos, 
siendo proporcionalmente la categoría de 
abundancia de 30 a 70% la más frecuente 
de observar en los árboles de J. 
deppeana. Por lo que, todos los árboles 
hembras en forma continua están produ- 


ciendo frutos de distintas categorías de 
maduración a lo largo del año. 

Además, durante los muéstreos no 
se registraron plantas sin gálbulos en la 
copa y no determinó una cierta época en 
donde los gálbulos maduros fueran más 
abundantes. Esto difiere en parte de la 
idea de que esta especie sólo madura los 
frutos durante los meses de agosto a 
octubre (Salazar-Figueroa, 2001; Vázquez- 
Yanes et al., 2002). No obstante, se 
requieren más estudios a largo plazo 
sobre la fenología del sabino, ya que la 
producción de conos puede estar relacio- 
nada con las condiciones ambientales en 
donde crecen las plantas (García et al., 
2002). Esto lo han demostrado Ortiz et al. 
(2002) con J. communis en seis locali- 
dades del sureste de España en donde la 
producción de conos decrece conforme 
aumenta la altitud. Asimismo, para el 
sabino J. deppeana, al igual que para 
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muchas gimnospermas, la producción de 
frutos puede variar de un año a otro 
(Powell, 1977, Jordano 1991; Arista et al., 
1997) y ello puede maximizar o reducir el 
esfuerzo reproductivo de las plantas para 
producir nuevamente frutos (Stephenson, 
1981). En este sentido, Ortiz et al. (2002) 
obtuvieron datos de los valores promedio 
de producción de conos en ramas de 
plantas hembras y determinaron que 
había diferencias entre años y dentro de 
un gradiente de altitud, registrando, en 
promedio, de 14 hasta 268 conos por 
rama. 

Con nuestros resultados, a pesar de 
no ser comparables con los de J. 
communis, se determinó que las ramas 
revisadas produjeron en promedio de 
cinco hasta 40 gálbulos. Diferencias que 
pueden deberse a un efecto del tamaño 
de la rama, dato omitido por Ortiz et al. 
(2002), ya que no mencionan la longitud 
de las ramas utilizadas en J. communis. 

Por otro lado, la gran cantidad de 
gálbulos colectados en las redes para los 
dos últimos registros fue sobresaliente. 
Esto puede deberse, por un lado, al 
efecto de las fuertes lluvias torrenciales 
acompañadas con ráfagas de viento 
durante junio y agosto, por lo que la 
acción climática puede facilitar el 
desprendimiento de los frutos en cierta 
forma, aún cuando no hayan completado 
su desarrollo, pues es frecuente el 
desprendimiento de gálbulos que no han 
completado su madurez, como son los 
frutos pequeños. 

En el caso de efecto materno, éste 
se puede deber también a un aborto 
selectivo de los embriones, como ocurre 
en muchas gimnospermas (Owens et al., 
1991), lo cual está relacionado con la 
disponibilidad y origen del polen cuando 
la densidad de individuos es baja. De ahí 
que la producción de frutos sea menor 
para Abies pinsapo (Arista y Talavera, 


1996). Además del efecto materno, el 
daño que causan los depredadores de 
frutos y semillas puede mermar a la 
población. 

En este último sentido, el destino y 
características de los gálbulos varían en 
cuanto a la cantidad y tamaño, ya que son 
los frutos de categoría chicos y medianos 
los menos depredados por insectos 
consumidores de semillas entre las 
distintas fechas de registro. Aunque, sería 
necesario realizar el seguimiento de 
frutos pequeños, medianos y grandes 
para determinar el momento del daño. A 
pesar de ello, esta información concuerda 
con lo reportado para algunas especies 
de insectos hembras de Megastigmus 
spermotrophus que parasitan óvulos 
fértiles y los cuales prefieren usar frutos 
de tamaño grande para ovipositar, de tal 
forma que garantizan las suficientes 
reservas alimenticias para las larvas (von 
Aderkas et al., 2005). 

Así, en otras especies de gimnos- 
permas el daño por insectos del género 
Megastigmus está más asociado a 
tamaños de cono más pequeños (Rappa- 
port y Roques, 1991; Rappaport et al., 
1993). Lo contrario ocurre en J. deppeana 
en donde son los gálbulos grandes los 
que tuvieron agujeros de salida de 
insectos consumidores de semillas. Esto 
se relaciona en cierto modo con las 
características de tamaño de los frutos 
(Góttsche, 1977) pues en su estudio 
experimental, documentó que cuando 
ofrecía gálbulos de enebro a hembras de 
M. bipunctatus, éstas visitaban frutos de 
tres a cuatro años de desarrollo, mientras 
que las de dos años de edad pasaban 
inadvertidas. Posiblemente esto se debía 
a que eran de menor tamaño. Siendo así 
los gálbulos de tres años los más selec- 
cionados para la oviposición, quizá 
porque esta etapa de desarrollo es la más 
adecuada para el crecimiento larval. Sin 
embargo, este patrón de daño por 
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insectos no es similar al de otras gimnos- 
permas, ya que por ejemplo, los conos de 
Cupressus sempervirens que son 
atacados por M. wachtli no difieren en 
cuanto a su tamaño con relación a los 
conos que no sufren daño (Guido et al., 
1995). 

En este sentido, los gálbulos de 
tamaño pequeño durante las cinco fechas 
de colecta fueron los que tuvieron menos 
orificios de salidas de insectos. Esto 
indica la preferencia por oviposición de 
los insectos en frutos que garanticen la 
supervivencia de las larvas hasta la fase 
de adultos. Sin embargo, la densidad de 
huevos por fruto no necesariamente está 
relacionada con el tamaño del fruto, tal y 
como se ha reportado en el caso de Acro- 
comia aculeata, debido a que el insecto 
Speciomerus revoili no selecciona los 
frutos grandes para ovipositar (Ramos, 
2001 ). 

En J. deppeana aún no se ha identi- 
ficado a la especie de insecto que daña 
las semillas. Pero la abundancia de este 
tipo de gálbulos no es baja, ya que de 
acuerdo con nuestras estimaciones, 
puede llegar a superar el 30%. Propor- 
ción que varía al compararla con otras 
especies de gimnospermas, como es el 
caso de conos de pino, a los que Eupi- 
thecia spermaphaga puede llegar a 
infestar hasta un 20% de los conos, 
aunque sólo destruye un máximo de 10% 
de las semillas en los conos infestados 
(Hedlin et al., 1980). 

Lo mismo ha sido reportado para 
Tsuga canadensis, la cual sufre daño en 
sus conos por E. mutana y M. hoffmeyeri 
de hasta un 21% por la primera especie y 
de 9 a 40% por la segunda, además de 
que también demuestran que los porcen- 
tajes de conos dañados son mayores que 
el de las semillas que son depredadas 
(Turgeon et al., 2004). En este caso, a 
pesar de no haber hecho estas estima- 


ciones, los gálbulos que tienen orificio de 
salida de insecto poseen dos de tres 
semillas que han sufrido el daño, por lo 
que el porcentaje de semillas dañadas 
también puede variar en este sentido 
(Sainos y Martínez, datos no pub.). 

En el caso de J. communis, García 
(1998) determinó una alta correlación 
entre la tasa de frutos depredados por 
avispas de Megastigmus y el total de 
semillas dañadas por planta. En este 
sentido, el orificio de salida del insecto en 
los gálbulos también significa muy posi- 
blemente que la selección de las hembras 
al ovipositar sobre los frutos depende de 
las características y viabilidad de las 
semillas, tal como lo reportan Herrera 
(1984) y Knight (1987). En el caso 
contrario, si la oviposición se realiza 
sobre gálbulos con semillas inviables, las 
larvas potencialmente mueren. Pero en J. 
monosperma las larvas de insectos que 
se encuentran en frutos que sufrieron 
abscición tienen más alta probabilidad de 
morir que aquellas que se encuentran en 
los gálbulos que aún están unidos al 
árbol. Este proceso de desprendimiento 
de frutos lo interpretan Fernández y 
Whitham (1989) como una estrategia de 
la planta para reducir los costos hacia 
gálbulos dañados y destinados a morir y 
además de que la abscición puede 
reducir el número de depredadores de las 
frutas y semillas. 

Bajo este contexto, una estrategia 
de evasión a la depredación, tanto de 
frutos como de semillas, la describen 
Fuentes y Shupp (1998), quienes docu- 
mentan que J. osteosperma produce 
frutos sin semillas como consecuencia del 
fenómeno de partenocarpia y esto reduce 
el efecto del ave depredadora de frutos y 
semillas Parus inornatus, proceso que 
también se observa en J. communis al 
reducir la depredación por Apodemus 
sylvaticus (García et al., 2000). En este 
sentido, la saciación es otra vía por la 
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cual puede reducirse el consumo de 
semillas y daño de frutos. Era de esperar 
que a más cobertura de copa de árbol, 
fuera mayor la producción de frutos 
sanos. 

Sin embargo, la densidad de 
gálbulos dañados no sobrepasó los 
100,000 frutos, aún cuando la cobertura 
de la copa aumentó y en cierta forma se 
esperaba un patrón parecido al de los 
frutos no dañados. Por lo tanto, estos 
resultados sugieren un posible efecto de 
saciedad, ya que a mayor número de 
gálbulos los depredadores se sacian y así 
el resto de los frutos escapan al potencial 
daño. 

Se desconoce si los frutos que 
permanecen sujetos a las plantas tienen o 
no perforaciones de salida de insectos, 
por lo cual no se tiene el punto de compa- 
ración. Aunque estos datos sugieren que 
si los gálbulos que se colectaron en las 
trampas son los que sufren abscición a 
consecuencia del daño por depredadores 
de semillas, entonces la categoría de 
frutos chicos con daño deberían ser más 
abundantes que los de tamaño mediano y 
grande, hecho que no ocurre en J. 
deppeana. 

Esto contrasta con lo que sugieren 
Janzen (1969, 1971b) y Queller (1985), 
en el sentido de que la depredación de 
semillas por insectos conlleva a una 
abscición selectiva de frutos jóvenes. 
Para este caso, los datos demuestran 
que la abscición de gálbulos es más 
intensa en los frutos de tamaño chico que 
no presentaban orificio de salida de 
insectos, lo cual contrasta con lo docu- 
mentado por los anteriores trabajos. 

Es así como el desprendimiento de 
gálbulos bajo la copa de J. deppeana 
contribuye a que los frutos estén disponi- 
bles para mamíferos que consumen e 
integran estas estructuras en su dieta. 


Por ejemplo, en el trabajo de Lezama 
(2007), se documenta que el conejo 
montés Sylvilagus cunicularius consume 
gálbulos debido a que en sus excretas 
están presentes hasta en un 10% las 
semillas de J. deppeana. Asimismo, 
Sainos (2007) también documenta que el 
cacomixtle Bassariscus astutus también 
consume los frutos. Por lo tanto, el hecho 
de que los frutos se desprendan de la 
planta madre promueve que estos dos 
mamíferos integren en su dieta a los 
gálbulos, pues durante los muéstreos en 
el sitio de estudio no se observaron aves 
frugívoras sobre el sabino, como se 
reporta para otras especies de Juniperus 
de Europa (Fuentes y Schupp, 1998). Por 
ello, el desprendimiento de frutos facilita 
que estén disponibles bajo la copa de los 
árboles, aunque una suposición es que 
conejos y cacomixtles consumen sólo 
aquellos gálbulos que no presentan perfo- 
raciones de salida de los insectos 
(Sainos, 2007). 

En conclusión, la producción de 
gálbulos en J. deppeana fue continua 
durante las cinco fechas de colecta y el 
número de árboles de sabino con catego- 
rías porcentuales de mayor abundancia de 
frutos (70 a 100%) se mantuvo entre 18 a 
26 plantas. Sin embargo, las categorías de 
30 a 50% y de 50 a 70% de frutos fueron 
las que presentaron el mayor número de 
plantas con gálbulos. Además, las catego- 
rías de 5 a 30% de frutos gálbulos no fue 
menor de 17 árboles, por lo que los Juni- 
perus hembras siempre mantuvieron 
producción de frutos. Asimismo, hay una 
amplia distinción entre la cantidad de 
gálbulos con y sin daño por insectos, lo 
cual está asociado con el tamaño o estado 
de desarrollo de los frutos. Finalmente, se 
requieren más investigaciones sobre el 
sistema reproductivo de J. deppeana rela- 
cionadas con aspectos de esfuerzo repro- 
ductivo y su asociación con procesos 
ecológicos como dispersión de semillas y 
depredación. 
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Características anatómicas y propiedades 
mecánicas de la madera de Oecopetalum 
mexicanum Greenm. & C.H. Thomps. (Icacinaceae) 
de la sierra de Misantla, Veracruz, México 

Mechanical properties and anatomical characteristics 
of Oecopetalum mexicanum Greenm. & C.H. Thomps. 
(Icacinaceae) wood of the la sierra de Misantla, Veracruz, México 

Maite Lascurain 2 , Guillermo Angeles-Álvarez 1 , 

Fernando Ortega Escalona 1 , Víctor Rubén Ordóñez Candelaria 2 , 

Mima Ambrosio 4 y Sergio Avendaño 3 

RESUMEN 

Se describe por primera vez la anatomía de la madera de Oecopetalum mexicanum (cachichín) 
Greenm. & C.H. Thomps., la cual es una especie que crece en los bosques y ecosistemas de café de la 
sierra de Misantla, Veracruz, México. Con el fin de conocer sus posibles usos, se determinaron las propie- 
dades mecánicas en flexión, del tipo: esfuerzo en el límite de proporcionalidad, trabajo unitario en el límite 
de proporcionalidad, módulo de elasticidad, módulo de ruptura y trabajo unitario a carga máxima; en 
compresión paralela a la fibra, como lo fue el esfuerzo máximo, esfuerzo en el límite de proporcionalidad 
y módulo de elasticidad. Todo ello, en su madera en condición verde. Se investigó también la densidad 
relativa básica de cada muestra, así como el diámetro, la distribución y densidad de los vasos, la abun- 
dancia y distribución del parénquima axial y radial. A partir de macerados, se identificó el tamaño de los 
elementos de vaso, tamaño de las fibras y contenidos celulares. La descripción anatómica se basó en la 
IAWA List of microscopio features for hardwood Identification (IAWA 1989). Los resultados obtenidos a la 
fecha indican que la madera del cachichín tiene resistencia similar a la del pino y, en algunas propiedades, 
incluso lo supera. Por lo tanto, es apta para la construcción, siempre y cuando se utilice todo el tronco y 
solo cuando sea posible aserrarla, se usarán piezas prismáticas aserradas. 

PALABRAS CLAVE: 

Anatomía de la madera, Oecopetalum mexicanum , Propiedades mecánicas en flexión, Sierra de 
Misantla, Veracruz, México. 


ABSTRACT 

The wood anatomy of Oecopetalum mexicanum Greenm. & C. H. Thomps. is described for 
the first time. This species grows in the forest and coffee ecosystem of the sierra de Misantla, Vera- 
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cruz, México. In order to recommend 
possible uses, mechanical properties 
essays of different types were carried 
out. These were: bending, like stress at 
proportional limit, unit work to propor- 
tional limit, modulus of elasticity, 
modulus of rupture, and unit work to 
máximum load in bending; paral leí 
compression to grain, like máximum 
strength, stress at proportional limit and 
modulus of elasticity. All these essays 
were performed in green wood condi- 
tion. For each sample, specific gravity 
was obtained. Vessel diameter, density 
and distribution, and estimation of 
abundance and distribution of axial and 
radial parenchyma were measured from 
permanent slides. Fiber and vessel 
elements length, as well as cell 
contents, were determined from mace- 
rated wood samples. Anatomical 
descriptions were made according to 
the IAWA List of microscopic features 
for hardwood identification (IAWA 
1989). The results so far obtained 
suggest that el cachichín wood is as 
resistant as pine wood and, for some 
mechanical properties, this is better 
than pine wood. Therefore, it is apt as 
building material, provided the whole 
bole is used; only if dimensions allowed, 
it could be saw for making prismatic 
pieces. 

KEY WORDS: 

Wood anatomy, Oecopetalum mexi- 
canum, Mechanical properties in 
bending, Sierra de Misantla, Veracruz, 
México. 


INTRODUCCIÓN 

Para hacer más eficiente el uso de nues- 
tros recursos naturales, es preciso contar 
con la información tecnológica que facilite 
la optimización de los procesos de su 
transformación. En el caso de Oecope- 
talum mexicanum Greenm. & C.H. 
Thomps. de la familia Icacinaceae, cono- 
cido con el nombre común de “cachichín”, 
su aprovechamiento como árbol de 
sombra de cafetales es muy conocido, al 
igual que sus frutos que son comestibles y 
representan una fuente de ingresos relati- 
vamente importante en la sierra de 
Misantla. Esta sierra se localiza en el 
estado de Veracruz, México y en una de 


sus localidades se realizó la colecta de 
material para realizar el trabajo de labora- 
torio, mismo que conforma la parte prin- 
cipal del presente estudio. La madera de 
estos árboles se utiliza para fabricación de 
muebles rústicos, cuando son derribados 
porque ya no producen suficientes frutos, o 
porque se requiere abrir claros en el 
cafetal o cambiar este cultivo por otra acti- 
vidad más rentable. No obstante que el 
uso de los árboles de cachichín es muy 
antiguo en la región Misanteca, hasta el 
momento se desconocen las propiedades 
anatómicas y tecnológicas de su madera. 
Con el fin de contribuir al conocimiento de 
esta madera, que permita su aprovecha- 
miento más integral, se presentan algunas 
de sus propiedades físicas y su descrip- 
ción anatómica. Con base en esto, se 
sugieren algunos usos adecuados a sus 
propiedades mecánicas. 

Botánica 

La familia Icacinaceae está conformada por 
aproximadamente 60 géneros y 400 espe- 
cies; la mayoría de éstas tienen distribución 
tropical pero algunas crecen en zonas 
templadas. En el estado de Veracruz 
(México) se registran las siguientes cuatro 
especies: Calatola mollis, C. laevigata, 
Mappia racemosa y Oecopetalum mexi- 
canum (Gutiérrez, 1994). 

Oecopetalum mexicanum fue 
descrita por Greenman y Thompson en 
1914 basándose en una colecta de 
Purpus de 1912, proveniente de Misantla, 
Veracruz. Esta especie aparte de crecer 
en Veracruz también se reporta para 
Chiapas y Guatemala (Gutiérrez, 1994). 
Se conoce con los nombres comunes de 
cacaté y cacaté de septiembre (Simo- 
jovel, Chiapas), cachichín (Misantla, 
Veracruz), jamacuquiaca (lengua zoque, 
Tapalapa, Chiapas) (Martínez, 1987). 
Oecopetalum greenmanii Standl. & 
Steyerm. se distribuye en México 


Madera y Bosques 13(2), 2007:83-95 

(Chiapas y Tabasco), Guatemala y Costa 
Rica. Actualmente no se cuenta con la 
información de usos o los nombres 
comunes que se le aplican a esta 
especie. 

De acuerdo con Gutiérrez (1994), 
O. mexicanum es un árbol o arbusto de 
2-25 m de altura, las ramas son glabras 
o pubescentes. Las hojas son simples, 
pecioladas, haz verde obscuro, envés 
verde claro, elíptico-lanceolado, de 15- 
25 cm de largo, 7-10 cm de ancho, 
glabras, margen entero, ápice apiculado. 
La inflorescencia es de tipo terminal o 
axilar, la flor posee cinco pétalos 
blancos, oblongo-lanceolados con 
anteras amarillas; el fruto es una drupa 
verde y café (al madurar), globosa, 
rugosa, glabra de 2 a 3 cm de largo, 1 a 
2 cm de ancho, glabra. La semilla es de 
color blanco. 


Área de estudio 

La sierra de Misantla también llamada 
de Chiconquiaco, está localizada en la 
parte central del estado de Veracruz, en 
la vertiente del Golfo de México (19° 52’ 
y 21° 02’ latitud norte y 96° 46’ y 97° 59’ 
longitud oeste). Está formada por un 
complejo de cañadas y barrancas con 
importantes cerros como el Algodón, 
Maquilquilitlán, San Miguel y Santa Rita. 
La vegetación y el clima de la sierra 
están conformados por tres unidades 
fisiográficas: pinares que corresponden 
al clima templado húmedo C(fm)B(c)g 
(Gutiérrez, 1993; Gómez-Pompa, 1966) 
y se encuentran a 2100-2500 msnm; el 
bosque caducifolio o mesófilo tiene un 
clima semicálido húmedo (A)C(fm) a(c) 
(Gutiérrez, 1993, Gómez-Pompa, 1966) 
y se desarrolla a 1100-2100 msnm. El 
tercer ecosistema es la selva alta 
subperennifolia y le corresponde un 
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clima cálido húmedo Af (m)(c) con 
lluvias todo el año y una altitud de 
400-1100 msnm (Gutiérrez, 1993, 
Gómez-Pompa, 1966). 

En la región de Misantla, el cachi- 
chin es más frecuente de observar 
desde los 400 a 1100 msnm, prospera 
en suelos húmedos, laderas pronun- 
ciadas y cerca de caídas de agua. 
Puede encontrarse como especie 
primaria en ecotonos de bosque caduci- 
folio o mesófilo y selva alta subperenni- 
folia. Algunas de las especies que 
caracterizan estos ecotonos son: Beils- 
chmiedia anay (escalón), Ulmus mexi- 
cana (olmo), Calatola mollis (calatola), 
Persea schiedeana (pagua), Pseudol- 
media oxiphilaria (tepetomate), Ficus 
cooki (en Chiapas, en la Lacandona le 
dicen “mutut” o “motut”, “chumiz” y 
“Jitzicui”; como a muchos otros Ficus, 
también se le conoce como “higo” o 
“higuera”) y Nectandra sp. (laurel). 
También se le localiza en la selva alta 
subperennifolia en la región misanteca. 
Las especies más comunes de este tipo 
de vegetación son: Brosimum alicastrum 
(ojite o ramón), Cedrela odorata (cedro), 
Guarea grandifolia (pozote), Asplundia 
I abela (cola de pescado), Alchornea lati- 
folia (calabacino), Beilschmiedia anay 
(escalón) y B. mexicana (escalanique), 
entre otras (Gutiérrez, 1993; Gómez- 
Pompa, 1966). Otra área donde el 
cachichín se desarrolla es en el sistema 
agroforestal de café, donde es utilizado 
como árbol de sombra. 

La colecta de las muestras de 
madera se realizó en la localidad Roca 
de Oro, municipio de Yecuatla, Veracruz. 
Sus coordenadas son 19° 52’ y 21° 02’ 
latitud norte y 96° 46’ y 97° 59’ longitud 
oeste, su altitud promedio es de 1000 
msnm y se ubica aproximadamente a 10 
kilómetros al sur de la ciudad de 
Misantla (Fig. 1). 
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Figura 1. 

Localización del sitio de colecta en el Municipio de Yecuatla, Veracruz. 


OBJETIVOS 

Determinar las propiedades mecánicas 
en flexión (Esfuerzo en el Límite de 
Proporcionalidad elp, Módulo de Elasti- 
cidad moe, Módulo de Ruptura, mor) y en 
compresión paralela a la fibra (elp, moe, 
mor) de la madera del cachichín ( Oeco- 
petalum mexicanum Greenm. & c.h. 
Thomps.) en condición verde. 

Identificar la densidad relativa básica de 
esa madera. 

Describir la anatomía de la madera y 
proporcionar sugerencias sobre su uso. 


MATERIALES Y MÉTODOS 

Pruebas mecánicas 

Las pruebas para determinar las propie- 
dades mecánicas de la madera, tanto en 
flexión estática como en compresión para- 
lela a la fibra, se llevaron a cabo utilizando la 
norma de astm DI 43-94 (astm, 2005a) de 
acuerdo con los incisos 6, 8, 9 y 21. Para la 
determinación de la densidad relativa básica 
de la madera se usó la norma astm D2395- 
94 (astm, 2005b) siguiendo el método “B”. 
En este trabajo sólo se realizaron los 
ensayos con material en verde y aplicando 
los métodos secundarios establecidos en la 


norma astm. 
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El equipo utilizado en las pruebas 
fue: 

• Máquina Universal de pruebas de 250 
kN. Con las celdas Interface 1220-50k y 
la 121 0-2k 

• Transductores de deformación Transtek 
0242 para medir deflexiones en flexión y 
Omega LD400-2.5 para las deforma- 
ciones en compresión paralela. 

Obtención de probetas 

Las probetas se obtuvieron de la madera 
colectada en dos parcelas con café, 
ubicadas en la localidad de Roca de Oro, 
municipio de Yecuatla, Veracruz, en donde 
se derribaron tres árboles de 11 años de 
edad y con un diámetro de 15 cm cada 
uno. Se obtuvo una troza de 125 cm de 
longitud por árbol. La tabla 1 detalla, para 
cada uno de los tres árboles muestreados, 
el número que se les ha asignado en las 
bitácoras y registros del LiNCE-Laboratorio 
de pruebas de productos maderables del 
Instituto de Ecología, A.C. También se 
anota la altura de los árboles, largo del 
fuste comercial y la categoría de dificultad 
que presenta el aserrado de su madera. 

Estudio anatómico 

Las muestras para este estudio se obtu- 
vieron también de los tres árboles colec- 
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tados. A cada uno se le cortó una troza de 
más de 125 cm de alto. Después se 
ajustó su longitud a esta medida con una 
motosierra y a partir de los 125 cm se les 
cortó una rodaja de 5 cm de espesor. De 
cada rodaja se obtuvieron dos cubos de 
1,5 cm por lado, perfectamente orien- 
tados de acuerdo con los planos de la 
madera. De ellos se obtuvieron cortes 
transversales, radiales y tangenciales de 
15 pm de espesor, usando un micrótomo 
de deslizamiento. 

Los cortes se tiñeron diferencial- 
mente con safranina “O” acuosa al 1% y 
verde rápido al 0,1% en alcohol de 96°. 
Se montaron en resina sintética disuelta 
en xilol. Se obtuvieron fotomicrografías 
digitales con un microscopio marca 
Nikon, modelo Eclipse E600. 

Los cubos empleados en la obten- 
ción de cortes histológicos se aprove- 
charon para obtener imágenes digitales al 
microscopio electrónico de barrido (meb). 
Para ello, con una navaja desechable se 
redujo su volumen un 60% y después se 
deshidrataron en una estufa a 60 °C, 
durante 72 horas. Una vez secos, se 
montaron en barriles de bronce de 6 mm 
de diámetro y se cubrieron con oro- 
paladio. Después en ellos se hicieron 
observaciones con un meb marca jeol 
modelo jsm-5600, utilizando un voltaje de 
aceleración de 20 KeV y un haz de elec- 
trones de 28 pm. 


Tabla 1 . Datos de la colecta de Oecopetalum mexicanum 


No. de 
árbol 1 

Altura 

Total 

cm 

Largo fuste 
comercial 
cm 

Aserrado 

192 

340 

130 

fácil 

193 

530 

190 

fácil 

194 

700 

200 

fácil 


1 Número consecutivo dado a las probetas en el LINCE-LP 
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Características anatómicas y propiedades mecánicas de ia madera de Oe copeta! um... 


Para obtener material disociado se 
obtuvieron astillas de madera de aproxima- 
damente 3 mm de espesor. Aproximada- 
mente diez de estas astillas se sumergieron 
en una solución de volúmenes iguales de 
peróxido de hidrógeno y ácido acético 
glacial en un tubo de ensayo, a 50 °C 
durante 24 horas (Johansen, 1940). Poste- 
riormente, los elementos celulares ya diso- 
ciados (fibras, elementos de vaso y parén- 
quima) se lavaron y tiñeron con Pardo de 
Bismarck al 2% en solución acuosa 
(Johansen, 1940). Después de teñido el 
material, se lavó abundantemente con agua. 
Las células disociadas se montaron en gela- 
tina glicerinada. En estas preparaciones se 
midieron las longitudes de fibras y la de 
elementos de vaso, así como el diámetro de 
sus punteaduras. Otros caracteres obte- 
nidos del material disociado fueron: diámetro 
total y espesor de pared de fibras, tipos de 
punteaduras de vaso a vaso y de vaso a 
radio. Las descripciones de los caracteres 
microscópicos de la madera se realizaron de 
acuerdo con la lista de iawa (1989). 


RESULTADOS 

Pruebas mecánicas 

La tabla 2 muestra los resultados de las 
pruebas de flexión estática en pequeñas 
probetas. En los casos donde no se 
presentan los valores de incertidumbre de 
los resultados, se debe a que el equipo de 
medición no cuenta con información para 
obtener estos valores. 

La madera 

Madera con albura y duramen de similar 
color, rosa hue 7.5YR 8/4 o amarillo rojizo 
hue 7.5YR 8/6 a 7/6. Carece de olor y 
sabor. Tiene veteado suave dado princi- 
palmente por los límites de las zonas de 
crecimiento y los radios. Su textura es 
mediana, el lustre alto, el hilo es recto y 
su dureza media. La figura 2 (A-F) 
muestra los planos característicos de la 
madera de cachichín. 


Tabla 2. Resultados de pruebas de flexión estática en pequeñas probetas. 

Norma astm DI 43 incisos 8. 


Condición: 

Especie: 


Verde 

Oecopetalum mexicanum 


Número (s) de árbol 192, 193 y 194 


RESULTADOS (SI) 

Clave 

Esfuerzo 
en el límite 
de 

proporcio- 

nalidad 

Trabajo 
unitario en el 
límite de 
proporcio- 
nalidad 

Módulo 

de 

elasticidad 

Módulo de 
ruptura 

Trabajo 
unitario a 
carga 
máxima 

Contenido 

de 

humedad 

Densidad 

básica 

relativa 


MPa 

N-mm/mm3 

GPa 

MPa 

N-mm/mm3 

% 

pa/vv 

FEV-192-A1 

32,5 

0,0079 

7,6 

60,4 

0,173 

113 

0,556 

FEV-1 92-A2 

32,0 

0,0073 

8,1 

60,9 

0,183 

117 

0,545 

FEV-193-A1 

30,7 

0,0056 

9,3 

65,0 

0,141 

107 

0,563 

FEV-1 93-A2 

32,1 

0,0070 

8,3 

58,4 

0,130 

107 

0,561 

FEV-1 93-A4 

37,0 

0,0085 

8,8 

66,1 

0,110 

105 

0,571 

FEV-1 94-A1 

33,7 

0,0067 

9,9 

65,3 

0,101 

101 

0,576 

FEV-1 94-A2 

40,3 

0,0090 

10,0 

69,9 

0,145 

102 

0,569 

Promedio 

34,1 

0,0074 

8,9 

63,7 

0,140 

107 

0,563 

d.e. 

3,4 

0,0011 

0,9 

4,0 

0,030 

6 

0,010 

c.v. 

10 

15 

10 

6 

22 

5 

2 

n 

7 

7 

7 

7 

7 

7 

7 

Incertidumbre 
de los 

± 0,683 

.. 

± 0,234 

± 0,937 

.. 

0.008 

0.0006 

resultados 
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Figura 2 

Características de la madera de Oecopetalum mexicanum (cachichín). 

A. Plano transversal, mostrando aparentes anillos de crecimiento. B. Detalle de A. 
C. Plano tangencial. D. Plano radial. Las flechas muestran los radios parenquimáticos. E. Detalle de un radio 
(flecha); F. Hilo recto de la madera. Escalas: A= 2mm; B = 20mm; C = lOmm; D = lOmm; E = 20mm; F = 5mm. 


Descripción anatómica 

La tabla 3 presenta en resumen las carac- 
terísticas anatómicas sobresalientes de la 
madera del cachichín. 


Porosidad 

Madera con porosidad difusa y zonas de 
crecimiento demarcadas por un cambio 
en la abundancia de poros, notándose 
donde una zona con muchos poros 
colinda con otra de menos y el ancho de 


ambas zonas es muy variable. Las zonas 
con menor cantidad de poros resaltan en 
oscuro porque en ellas además es muy 
escaso el parénquima axial (Fig. 3A). Los 
poros son de contorno anguloso a 
redondo u oval (Fig. 3B) y se disponen de 
la siguiente manera: solitarios (55%), 
múltiples radiales de 2 a 5 (5%), múltiples 
oblicuos de 2 a 3 (21 %), múltiples tangen- 
ciales de 2 a 4 (15%) y los agregados de 
3 a 7 (4%). En cuanto a abundancia, 
existen en promedio 28 poros por mm 2 y 
su diámetro tangencial es de 63 im en 
promedio. 
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Tabla 3. Resultados de pruebas de compresión paralela a la fibra en pequeñas 

probetas. 

Norma astm DI 43 incisos 8. Condición verde 


RESULTADOS (SI) 


Esfuerzo 


RESULTADOS (MKS) 


Esfuerzo 


Clave 

en el 
límite de 
proporcio 
nalidad 

Módulo de 
elasticidad 

Esfuerzo 

máximo 

Contenido 

de 

humedad 

Densidad 

básica 

relativa 

en el 
límite de 
proporcio 
-nalidad 

Módulo de 
elasticidad 

Esfuerzo 

máximo 


MPa 

GPa 

MPa 

% 

pa/vv 

kg/cm 2 

Mg/cm 2 

kg/cm 2 

CPV-192A-01 

- 

7,8 

24,5 

109 

0,531 

- 

79 

249 

CPV-192A-02 

16,5 

8,5 

23,6 

96 

0,543 

169 

87 

241 

CPV-193A-01 

16,2 

9,1 

24,7 

91 

0,569 

165 

92 

252 

CPV-193A-02 

15,7 

8,2 

21,6 

98 

0,567 

160 

83 

220 

CPV-194A-01 

15,0 

10,0 

25,9 

79 

0,553 

153 

102 

264 

CPV-194A-02 

16,2 

11,3 

24,1 

87 

0,560 

165 

116 

246 





Promedio 

15,9 

9,2 

24,1 

93 

0,554 

162 

93 

245 

Desv. 

0,6 

1,3 

1,4 

10 

0,015 

6 

14 

14 

c.v. 

4 

15 

6 

11 

3 

4 

15 

6 

n 

5 

6 

6 

6 

6 

5 

6 

6 


No se reportan valores en el 
caso de compresión paralela. 
Incertidumbre Prueba no acreditada. Sin 
délos calibrar transductor de 

resultados deformación 


± 0.008 


± 0.0006 


No se reportan valores en el caso 
de compresión paralela. Prueba 
no acreditada. Sin calibrar 
transductor de deformación 


Elementos de vaso 

Los elementos de vaso presentan una 
longitud de 2020 pm en promedio. La 
mayoría tienen en sus extremos cola o 
lígula (Fig. 4 A-D) y placa de perforación 
escalariforme con 20 a 55 barras 
delgadas (en promedio 36), muchas de 
ellas se anastomosan (Fig. 4D, F, G). Las 
punteaduras de vaso a vaso son de dos 
tamaños: las pequeñas son areoladas y 
ornamentadas, su disposición es irre- 
gular, aunque tienden a presentarse 
como opuestas o alternas (Fig. 4C, E); las 
grandes son alargadas simples o escala- 
riformes y su diámetro promedio es 
mayor a 10 im (Fig. 4A). Este último tipo 
de punteaduras prevalece en los campos 
de cruce o áreas de la pared de los vasos 
que hacen contacto con los radios, donde 
también puede haber algunas puntea- 
duras pequeñas, areoladas (Fig. 4E). 


Parénquima axial 

Es apotraqueal difuso y abundante (Fig. 
3B). En el plano tangencial las series de 
parénquima axial son de 2 a 15 células, 
con 8 en promedio. 

Parénquima radial 

Los radios son uniseriados (40%) y multise- 
riados (60%) de 2 a 1 0 células de ancho (Fig. 
3C,D). Los primeros son de 2 a 43 células de 
alto, con un promedio 15. Los segundos 
pueden o no tener porción uniseriada pero su 
área multiseriada está compuesta de una 
combinación aleatoria de células procum- 
bentes, cuadradas y erectas. Las cuadradas 
pueden mostrarse como células marginales 
cuando no hay porción uniseriada (Fig. 4C, 
F) o como células envolventes, formando 
una vaina incompleta a completa que deli- 
mita la porción multiseriada. 
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Figura 3 

Características anatómicas de la madera de Oecopetalum mexicanum. 

A. Corte transversal. La porosidad es difusa. Las flechas delimitan una zona de crecimiento. B. Detalle de A, 
mostrando los vasos (v). La flecha apunta a una fibra de pared gruesa. Radios (r) con abundante almidón. C. Corte 
tangencial. Fibras (f). Radios (r) con células alargadas en los extremos (flecha). D. Vista tangencial al MEB. La 
flecha señala los granos de almidón en los radios. Fibras (f) flanqueando a los radios. E. Fibra traqueida con punte- 
aduras simples, alargadas (flecha). F. Corte radial. Ambos extremos formados por varias células erectas (flechas). 
G. Vista radial al MEB. Células erectas (flecha), con algunos granos de almidón. H. Fibra libriforme, con puntea- 
duras simples (flecha). 

Escalas: A =800im; B = lOOim; C = 80im; D, E = 20im; F= 701m; G = 100 im; H = lOim. 
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Tabla 4. Características anatómicas microscópicas de Oecopetalum mexicanum. 



CARÁCTER 

X 

MODA 

MÍN. 

MÁX. 

DE 

VASOS 

No./mnV 

28 

25 

13 

45 

5.05 

* Diámetro tangencial 

63 

66 

30 

90 

12.03 

No. de barras 

35 

35 

18 

55 

6.54 

* Longitud 

2020 

1800 

1296 

3348 

342.05 

JIMA RADIAL 

Uniseriados 

No./mm tangencial 

13 

12 

8 

17 

1.96 

* Altura 

942 

1000 

120 

3000 

537.40 

* Ancho 

21 

18 

12 

30 

4.59 

No. de células 

15 

9 

2 

43 

8.76 

PARÉNQl 

Poliseriados 

* Altura 

2774 

2000 

320 

7200 

1378.70 

* Ancho 

155 

120 

80 

480 

56.95 

No. de series 

6 

5 

2 

10 

1.77 

FIBRAS 

* Longitud 

2705 

2880 

1680 

3888 

369.98 

* Diámetro 

41 

40 

24 

60 

5.56 

* Grosor de pared 

13 

12 

8 

24 

2.20 


* Valores en pm, X = Media, DE = Desviación Estándar 


Los radios uniseriados tienen en 
promedio 942 pm de alto y 21 pm de 
ancho; los multiseriados tienen en 
promedio 2774 pm de alto y 155 pm de 
ancho. Tomados en conjunto, los radios 
multiseriados y uniseriados se cuentan en 
promedio 13/mm. 

Fibras 

Las fibrotraqueidas tienen punteaduras 
areoladas, con areolas pequeñas. Son 
largas, en promedio de 2705 pm, 41 pm 
de diámetro y 13 pm de grosor de pared 
(Fig. 3E). Las fibras libriformes tienen 
punteaduras simples (Fig. 3H). 


Contenidos 

Se presentan profusamente granos de 
almidón en los parénquimas axial y radial 
(Fig. 3D-E, 4G). En las células del parén- 
quima radial, hay algunos cristales cúbicos, 
romboidales o poliédricos, además abun- 
dancia de taninos (Fig. 3B,C). 

DISCUSIÓN 

La anatomía de la madera del cachichín 
muestra ciertas peculiaridades, que 
pueden estar incidiendo en las propie- 
dades mecánicas. Las fibras son largas 
(iawa, 1989), con una longitud promedio 
de más de 2 mil pm y un grosor de pared 
de 13 pm. Esta cualidad le confiere a la 




Figura 4. Elementos de vaso. A. Sin cola o lígula. B. Con colas o lígulas en ambos 
extremos. C. Terminación sin cola. D. Terminación con cola. E. Punteaduras de vaso a 
vaso (flecha izquierda); de vaso o radio (flecha derecha). F. Vista radial al meb. 
Elementos de vaso (v) con plastinas escalariformes (flechas). G. Platina escalariforme 

con granos de almidón (flechas). 


Escalas: A, B: 120 jam; C, D: 60 |jm; F: 200 ¡im; G: 20¡jm 


madera resistencia mecánica a la 
compresión, flexión y tensión (Bárcenas- 
Pazos et al., 2005). Como las fibras 
predominan en la composición anatómica 
de la madera de esta especie, es de 
esperarse que sean éstas las que deter- 
minen las propiedades mecánicas de la 
madera en sentido paralelo a la fibra. Por 
otro lado, los radios parenquimáticos de 
esta madera son altos y gruesos (iawa, 


1989), coalescentes, abundantes. Los 
radios influyen negativamente en la resis- 
tencia mecánica a la compresión perpen- 
dicular al grano y hacen, además, que la 
madera se contraiga desigualmente 
durante el secado. Según nuestra expe- 
riencia (datos no publicados) la contrac- 
ción en cubos de 2 cm de lado, causa una 
deformación notable en la madera, en 
sentido tangencial. Por lo tanto, se debe 
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diseñar una estrategia de secado apro- 
piado, para evitar la deformación de las 
piezas fabricadas con esta madera. Por 
otro lado, el color y lustre atractivo de esta 
madera, así como su textura fina, la 
hacen un buen candidato para fabricar 
chapa y muebles. 

Se puede decir que la madera de 
cachichín tiene una resistencia mecánica 
similar a la de los pinos. Por ejemplo, 
comparándola con tres especies de pino 
provenientes del estado de Oaxaca, las 
propiedades en flexión del cachichín son 
similares a las de Pinus pseudostrobus 
(Ordóñez et al., 1989), salvo en el módulo 
de elasticidad, que tiene un valor más 
bajo. Sus propiedades mecánicas en 
compresión paralela son más parecidas 
al Pinus patula (Ordóñez et al., 1989). Se 
debe tener en cuenta que la madera es 
un material con alta variabilidad en sus 
propiedades, por lo que se puede 
aceptar, con las reservas del caso, que la 
madera de cachichín es apta como mate- 
rial de construcción, siempre y cuando se 
utilicen los rollizos y sólo cuando sea 
posible aserrarla se usarán piezas pris- 
máticas aserradas con dimensiones 
adecuadas para estructuras. Otra posibi- 
lidad es su uso para la elaboración de 
muebles o utensilios pequeños, por su 
color y veta atractivos. 

Cabe mencionar que los valores de 
resistencia obtenidos en el presente 
estudio fueron obtenidos en condición 
verde, por lo que se debe tener mucho 
cuidado al interpretarlos. Es bien sabido 
que la resistencia de la madera aumenta 
exponencialmente a medida que pierde 
humedad (Zink-Sharp, 2003). Sin 
embargo, cada especie se comporta de 
diferente manera, y la relación contenido 
de humedad-resistencia no se puede 
aplicar de manera generalizada. 


Como valores indicativos para la 
relación entre las propiedades mecánicas 
en condición verde y condición seca, se 
han obtenido datos de la madera de 20 
especies tropicales de la selva Lacan- 
dona (Bárcenas, 1995). En ese estudio se 
obtuvieron valores de 0.68 para el mor en 
condición verde contra 0.83 para condi- 
ción seca; para compresión paralela al 
grano se obtuvieron valores de 0.55 en 
condición verde y de 0.86 en condición 
seca. Esto coincide con la aserción de 
Zink-Sharp (2003) de que la compresión 
perpendicular al grano es la propiedad 
mecánica más influenciada por el conte- 
nido de humedad, seguida por la compre- 
sión paralela al grano. Además, debido al 
carácter anisotrópico del material, los 
cambios en contenido de humedad de la 
madera afectan a cada una de las propie- 
dades mecánicas de la madera en dife- 
rente manera. 
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Monitoreo de ecosistemas con estrategias 
geoestadísticas, una aplicación de gran 
escala en Jalisco, México 

Ecosystem monitoring through the use of geostatistic 
strategies, a large scale application in Jalisco, México 

José Germán Flores Garnica* 1 , Martín Alfonso Mendoza B .* 2 

y Celedonio Aguirre Bravo * 3 

RESUMEN 

El proyecto de Inventario y Monitoreo de los Recursos Naturales de Jalisco, México (IMRENAT), 
es un proceso innovador que propone una aplicación práctica de las herramientas de modelado esta- 
dístico geoespacial para reportar creíblemente la condición del territorio. Aprovechando las intercone- 
xiones que existen entre atributos del terreno, la vegetación, las actividades humanas y lo detectable 
con sensores remotos multiespectrales, se pueden construir modelos estadísticamente válidos que 
reflejan las relaciones espaciales y las tendencias dinámicas de los atributos del territorio. Este tipo 
de modelos permiten generar productos de geoinformación de mayor utilidad que la información 
geográfica de naturaleza cartográfica, y que la información pictórica proveniente de interpretar 
imágenes de sensores remotos. Estos productos geoespaciales tienen la capacidad de describir 
distintos niveles de resolución y de escala. 
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ABSTRACT 

The Inventory and Monitoring of Natural Resources of Jalisco, México (IMRENAT), is an inno- 
vative process suggesting a practical application for geospatial statistical modeling to credibly report 
the condition of land resources. Taking advantage of interconnections among terrain attributes, vege- 
taron, human activities and multispectral remóte sensing information, it is possible to develop statisti- 
cally valid models depicting land resources spatial attributes and dynamic trends. These types of 
models enable us to create geoinformation products of greater practical utility than regular cartogra- 
phic products, or geographic information derived from imagery interpretation. Geospatial products 
have the capability to describe different scales and resolution levels. 
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INTRODUCCIÓN 

Informar sobre el estado actual y los 
cambios en el ambiente, así como sobre 
las condiciones del mundo (p.ej. fao, 
2006), es un asunto de tal relevancia 
actual, que se ha creado de hecho una 
nueva rama de actividad conocida como 
monitoreo. El avance del conocimiento en 
esta rama ocurre no únicamente por su 
importancia, también es un vigoroso 
proceso impulsado por la dificultad 
técnica y humana que conlleva. 

Las principales metas buscadas en 
el monitoreo generalmente se relacionan 
con la credibilidad de la información resul- 
tante sobre la dinámica de los atributos 
del terreno, y en su momento, la posibi- 
lidad de sustentar técnicamente la calidad 
de esa información (Busch, 2005; 
Schreuder et al., 2003). 

Gran parte del impulso a la innova- 
ción y a la adopción de ideas y tecnología 
de monitoreo es motivado por la defi- 
ciente satisfacción de las necesidades y 
expectativas de los usuarios que han 
recurrido a los medios geoinformáticos 
tradicionales disponibles. En particular se 
ha cuestionado el papel de la información 
que proveen las agencias informativas de 
los gobiernos nacionales (Longhorn, 
2005). De esta información comúnmente 
se ha criticado la omisión de evidencias 
sobre la validez y procesos de validación 
de la geoinformación. 

Entre las tendencias fuertes que 
empiezan a mostrar alientos, se tiene la 
posibilidad de que la tecnología de 
cómputo actual ya permite interoperabi- 
lidad de una multiplicidad de bases de 
datos distribuidas en localidades y juris- 
dicciones distintas. Otra tendencia impor- 
tante en el monitoreo es la que pretende 
cambiar paradigmas geográficos, y en 
particular, se plantea un nuevo posiciona- 
miento mucho más modesto en las priori- 
dades de la sociedad moderna para el 


arte de la cartografía como modo de 
representación de la realidad geográfica 
(Brodersen, 2005). Esta nueva tendencia 
busca presentaciones de la geoinforma- 
ción que sean mejor adaptadas al 
contexto cultural y a las circunstancias del 
usuario, y que no impongan limitaciones o 
prejuicios producto de la tradición del arte 
cartográfico, o producto de la forma en 
que se definen y se representan las cate- 
gorías que son el medio para representar 
a los geodatos. 

Una de estas nuevas formas de 
generar productos de geoinformación es 
la que resulta de aplicaciones de geoesta- 
dística y modelado geoespacial 
(Schreuder et al., 2003). Este tipo de apli- 
caciones ha sido reportado en literatura 
(Kalas et al., 2003; Reich et al., 2004), en 
escalas de tipo experimental y ámbito 
temático limitado a uno o unos pocos 
asuntos. Hay dos casos reportados 
donde se tienen en marcha proyectos de 
gran envergadura destinados a reportar la 
condición de los ecosistemas en amplias 
zonas territoriales, una de ellas es Colo- 
rado Digital (Pellicane y Burns, 2005), y el 
otro, es el proyecto de Inventario y Moni- 
toreo de los Recursos Naturales de 
Jalisco (imrenat) en Jalisco, México 
(Alonso et al., 2006). 

Es intención de este trabajo dar 
evidencia de estos avances mediante el 
uso de un ejemplo práctico, que en este 
caso es el proyecto de imrenat. En el 
2001, el estado de Jalisco decidió 
proceder con urgencia a llenar demandas 
de geoinformación para las que eran insu- 
ficientes las cifras de instituciones nacio- 
nales como las de inegi y semarnat e 
internacionales como la fao. De ello 
surgió la necesidad de crear un instru- 
mento para proveer tal geoinformación 
bajo las directrices anotadas arriba. 
imrenat ha reportado ya un buen número 
de resultados que muestran las cuali- 
dades y capacidades de productos de 
geoinformación y geosíntesis que siguen 
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la pauta que marca la frontera del conoci- 
miento en esta materia. Un propósito 
colateral de este trabajo es difundir 
algunos de estos resultados en un 
formato que exhibe las cualidades y 
potencialidades de este diseño de moni- 
toreo de ecosistemas. 

ESTRATEGIA TEÓRICA DE IMRENAT 

En la figura 1 se muestra esquemática- 
mente la lógica que soporta metodológi- 
camente el trabajo de imrenat. Las impli- 
caciones lógicas de esta imagen aluden a 
que los elementos de carácter más 
estable, como la topografía, geología y 
suelos, son explicación eficaz de los 
procesos y cualidades superiores, como 
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sería el clima, que a su vez puede ser 
usado, junto con las capas anteriores 
para explicar la vegetación y fenómenos 
circunstanciales como el fuego. Todas 
estas capas temáticas, una vez acumu- 
lada suficiente porción de secuencias en 
tiempo de ellas, son base explicativa de 
los fenómenos antropogénicos como el 
uso del suelo y las variables geográficas 
de origen cultural, histórico, jurisdiccional 
y administrativo. Los elementos tangibles 
y los interpretados en esa imagen 
requieren de un marco abstracto de refe- 
rencia, representado aquí por la imagen 
central de un globo dentro de una retícula 
de coordenadas, que permita hacer las 
operaciones lógicas y matemáticas para 
unir los distintos tipos de información. 





X = Sistémas físico objeto de monitoreo (i=l,...n). 

Y = Aspectos de X objeto de monitoreo (i=l...n). 

Z = Datum geográfico de X y Y (e.g., latitud, longitud, altitud, 
tiempo). 


Figura 1. Esquema conceptual del proyecto imrenat 
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El aspecto relevante aquí es el fuerte 
apoyo que avances recientes en materia 
de estadística geoespacial (Schreuder et 
al., 2003), permiten manejar y potenciar la 
enorme masa de datos que implica un 
censo del territorio, implícito en el barrido 
de los sensores remotos multiespectrales. 
A propósito de esto, cabe mencionar que 
en imrenat se ha usado Landsat 7 TM, 
2002, 2003 y 2004, serie que se irá 
complementando con el tiempo con otras 
imágenes para 2005 en adelante. 

Con este gran poder estadístico se 
puede modelar el comportamiento de 
variables de la “realidad a nivel del suelo”, 
y lograr estimadas no sesgadas, de alta 
calidad y poder estadístico. Esta “realidad 
a nivel del suelo” es, en principio, 
producto de mediciones de una red de 
sitios permanentes, puesta en campo con 
un diseño que permite el estudio de rela- 
ciones espaciales. Para este último fin, la 
muestra se asigna en proporción a la 
variabilidad que los datos de campo 
exhiben respecto a las variables usadas 
en estratificación. Lamentablemente, los 
datos de los sitios en el terreno son sólo 
unas pocas cifras de eventos específicos, 
pero con los modelos temáticos apoyados 
en la información satelital y procedi- 
mientos de varianza forzada (tipo boots- 
trapping, jacknifing. Schreuder et al., 
2004), se da lugar a representaciones 
modeladas de estas capas temáticas con 
estimadores para el total de pixeles 
dentro del área de estudio. Estas repre- 
sentan no tanto a la situación del pixel en 
el momento de la medición, sino a la 
población de valores teóricamente plausi- 
bles en ese pixel para esa variable. 

Los procedimientos concretos para 
el modelado aludido se hacen por vía de 
árboles binarios de clasificación, para las 
variables binomiales y categóricas, y 
árboles de regresión y modelos de super- 
ficie de tendencia (Reich, 2005). Estos 
procedimientos se validan después 
mediante simulación Montecarlo que 


alimenta un proceso recursivo de 10 
ciclos de validación cruzada (Efron y Tibs- 
hirani, 1993), además de pruebas esta- 
dísticas de bondad de ajuste propias de 
modelos espaciales (Agterberg, 1984). 
Finalmente, los procesos y patrones de 
tendencia en escala corta (resolución de 
10 m, Reich et al., 2005b), se estudian 
modelando la dispersión de errores de los 
modelos de escala gruesa (resolución de 
30 m), mediante modelos de regresión y 
modelos de clasificación binaria, reite- 
rando hasta que no se puedan agregar 
variables con poder explicatorio estadísti- 
camente significativas, al modelo que 
predice los errores. 

Cabe aclarar que a lo largo de todo 
este proceso se mantiene una constante 
vigilancia de parte del analista para detectar 
elementos que requieran acción correctiva 
en alguna parte del proceso de modelado, y 
en especial en la porción de poda de los 
árboles binarios y árboles de regresión. 
Para auxiliar en esta tarea al analista, se 
desarrolló un procedimiento mecanizado 
que genera imágenes cartográficas de 
etapas intermedias del modelado, dentro de 
las cuales es posible examinar el modelo 
para visualmente detectar incongruencias o 
aspectos que puedan mejorar. Este proceso 
ha sido codificado en un programa llamado 
geosis (Solano et al., 2005), escrito en C+ 
en un entorno Linux (Red Hat Inc., 2006), 
para operar en estaciones de trabajo de 
gran capacidad, si bien ya se ha migrado a 
un formato que opera dentro de plataformas 
Windows (Microsoft Corp., 2006). geosis 
mecaniza el proceso de extracción de datos 
espectrales y de elevación, así como la 
compilación y edición de datos de campo y 
de fuentes secundarias, y prepara las 
matrices que alimentan el modelado esta- 
dístico. Una vez que se tienen preparadas 
estas matrices de datos, se usa un paquete 
estadístico (S-plus 7 Insightful Corp., 2006) 
para ajustar los modelos. A la salida del 
proceso estadístico, geosis recupera los 
reportes de resultados, resuelve los 


Madera y Bosques 13(2), 2007:97-104 

modelos estadísticos, y con sus predic- 
ciones despliega las representaciones 
cartográficas de los modelos, para conside- 
ración del analista, quien repite el proceso 
con modificaciones las veces que su criterio 
considere necesario para obtener un 
producto aceptable. Finalmente, se corren 
las librerías que fueron diseñadas en S-plus 
para fines de validación cruzada y mode- 
lado de errores. Si los parámetros de ajuste 
del modelo son satisfactorios, geosis toma 
el modelo obtenido, lo aplica a la malla de 
pixeles del total del área de estudio, genera 
predicciones pixel por pixel, y produce una 
imagen visual que aparenta ser un mapa, 
pero que en realidad es la superficie 
predicha por el modelo para la variable 
modelada. Esta imagen es estudiada 
visualmente para detectar posibles incon- 
gruencias, errores y necesidades de modifi- 
cación del modelo, que en caso de aparecer 
obligan a repetir el proceso para enmendar. 

De entre los varios resultados dispo- 
nibles al momento se tienen 24 capas 
temáticas de la zona central de Jalisco 
(Reich et al., 2005b), que fue usada como 
estudio piloto antes de emprender el esta- 
blecimiento de la red estatal de sitios 
permanentes de monitoreo. También se 
tiene para todo el estado de Jalisco, 
modelos de elevación (provistos por 
cortesía de inegi y nasa), zonas climá- 
ticas, y vegetación. El principal resultado 
ha sido un primer reporte de variables 
ambientales (Alonso et al., 2006), que 
contiene una variedad de temas mode- 
lados que enfatizan la situación de los 
recursos forestales estatales. Este 
énfasis en recursos forestales fue 
pensado para agregarse a la serie de 
previos inventarios forestales que han 
sido realizados a nivel federal y estatal. 

RESULTADOS Y DISCUSIÓN 

Una presentación amplia de resultados 
de este trabajo se presenta en Alonso et 
al. (2006). De este trabajo y a manera de 
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ilustrar las cualidades de los modelos ya 
disponibles, se extrajo el caso de la 
distribución de las coberturas del suelo 
(Figura 2). La imagen de este modelo se 
pone a la par de la vegetación como la 
reporta inegi, que si bien se trata este 
último de un estudio ya viejo (datos de 
2000), y a escala pequeña (1:250 000), 
sirve para mostrar de forma gruesa y 
algo exagerada dos atributos básicos de 
los productos de imrenat: la apariencia 
general del modelo de imrenat es más 
proclive a gradientes, a mostrar cierta 
apariencia granulosa en la distribución 
espacial de los atributos modelados, y a 
evitar la apariencia de fronteras nítidas 
que es producto del uso de polígonos 
fundados en interpretación de la imagen 
por expertos intérpretes, tal como 
despliega el caso de inegi. El otro 
aspecto, quizá el más importante de 
resaltar, es que al margen de las diferen- 
cias tecnológicas, de escala y obsoles- 
cencia, los patrones generales de distri- 
bución espacial son casi idénticos, lo 
que es un elemento empírico que 
sugiere implícitamente cierta confiabi- 
lidad de resultados imrenat. Por 
supuesto, la parte objetiva de esta 
segunda observación es la que resulta 
de los análisis estadísticos de los resul- 
tados del modelo. En este caso, el 
modelo mostró una confiabilidad general 
de 82,9% (Reich et al., 2005a). 

Admitiendo que los recursos fores- 
tales son motivo de especial atención, 
también se presenta como ejemplo de 
resultados meset la Figura 3, que 
despliega en tonalidades crecientes de 
gris, la distribución geográfica de la espe- 
sura de las masas arboladas de Jalisco. 
La espesura o densidad del arbolado en 
este caso está representada por el área 
basal (m 2 /ha). De nuevo la cualidad rele- 
vante en este otro resultado es la fina 
presentación de gradientes y clara distri- 
bución de la variabilidad de este indi- 
cador. 
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COBERTURAS DE SUELO 

I I BOSQUE DE PINO 
I I BOSQUE DE PINO - ENCINO 
I I BOSQUE DE ENCINO 

I I SELVA BAJA 

BOSQUE DE ABIES O CEDRO 
I I AGRICULTURA (O PASTIZAL) 

I I PASTIZALES 

I I HUMEDALES 

» — *1 CUERPOS DE AGUA 

ASENTAMIENTOS HUMANOS (INEGI) 
I 1 SUELO DESNUDO 



IMRENAT 


Figura 2. Cubierta del suelo modelada en IMRENAT, en relación a una presentación 
vectorial (inegi, serie 2, escala original 1:250 000). 

Datos INEGI: http://mapserver. inegi. gob.mx/map/datos_basicos/vegetacion/?c=556 


CONCLUSION 

imrenat es un esquema de monitoreo 
para el seguimiento espaciotemporal de 
las variables relacionadas con la susten- 
tabilidad de ecosistemas terrestres que 
ha mostrado, con sus resultados 
parciales disponibles en 2006, un poten- 
cial atractivo para representar fenómenos 
que tienen patrones geográficos de 
comportamiento modelables geoestadísti- 
camente. 

Una cualidad lateral en el diseño de 
imrenat es su capacidad de incorporar 
datos de fuentes diversas, y potencial 
interoperabilidad con otros esfuerzos y 
proyectos similares. 
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Figura 3. Distribución modelada del área basal del arbolado de Jalisco (m 2 /ha). 
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